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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. ‘The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
typed using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespectively of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zeits. f. Astrophys., 15, 239, 1938. 
A. S. Eddington, Internal Constitution of the Stars, Cambridge, p. 182, Table 24, 1926. 


3. Notation.—Authors should conform closely to the recommendations of Commission 
3 of the International Astronomical Union (Trans. J.A.U., Vol. VI, p. 345, 1938). Council 
has decided to adopt the I.A.U. 4-letter abbreviations for constellations where contraction 
is desirable (Vol. IV, p. 221, 1932). 


4. Diagrams.—These should be drawn about twice the size required in print 
and prepared for direct photographic reproduction except for the lettering, which 
should be inserted in pencil. Legends should be given in the manuscript 
indicating where in the text the figure should appear. Blocks are retained by the 


Society for 10 years; unless the author requires them before the end of this period they 
are then destroyed. 


5. Tables—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts——When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 
8. When submitting papers authors are requested to indicate whether they will be 


willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1950 OCTOBER 13 
Professor W. M. Smart, President, in the Chair 


The President referred to the death on 1950 September 21 of Professor E. A. 
Milne, and read the following tribute paid to him by the Council at its meeting 
that afternoon (the Fellows standing) :— 


“In the sudden death at the early age of fifty-four of Professor Edward 
Arthur Milne, the Royal Astronomical Society has suffered a grievous loss, 
one shared not only with astronomers the world over but with the whole world 
of science. Milne’s early work on radiation problems connected with the outer 
atmospheres of the Sun and stars won in 1935 the award of our Gold Medal, 
the highest honour that the Society can bestow. In his later years Milne devoted 
his keen intellectual faculties to a persistent and courageous attack on the 
problem of the ultimate nature of the physical universe. His book “‘ Kinematic 
Relativity ” will provide food for keen discussions for many years to come. 

‘“‘ As President, as Member of Council over many years, Professor Milne 
served the Society faithfully and loyally, never sparing himself, despite failing 
health, when asked to undertake any duty for the Society. His loss will be 
severely felt by our Fellows and the Council knows that they will wish to 
associate. themselves with an expression of deep sympathy with his family 
in their tragic bereavement.” 


The election by the Council of the following Fellows was duly confirmed :— 

Norman Robert Blackmore, 12 Victoria Street, Middle Brighton, Melbourne, 
Australia (proposed by E. B. Walton); 

Dennis Jack Fulcher, 19 Gloucester Road, Ipswich, Suffolk (proposed by 
E. H. Collinson); 

John Barrows Irwin, Goethe Link Observatory, Indiana University, 
Bloomington, Indiana, U.S.A. (proposed by W. A. Hiltner); and 


Eric Smith, 6 Fiennes Road, Herstmonceux, Sussex (proposed by W. A. 
Scott). 


Seventy-seven presents were announced as having been received since the 
last meeting, including :— 


British Astronomical Association, Index to the Association’s Fournal for 
Volumes 1-50 (presented by the Council of the Association); and 


Yusuke Hagihara, Foundations of Celestial Mechanics (presented by the 
author). 
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THE ELECTRICAL PHOTOMETRY OF STARS AND NEBULAE 


George Darwin Lecture, delivered by Professor Foel Stebbins on 
1950 October 13 


My feeling of pleasure on being with you today is mixed with a sense of 
responsibility, not to mention humility. Dr George E. Hale once said that the 
meetings of the Royal Astronomical Society are the most stimulating astronomical 
meetings in the world ; and I heard Professor Willem de Sitter say in this room 
that the Royal Astronomical Society acts as astronomical conscience, astronomers 
come from everywhere to confess what they have been doing. In the present 
case my task is lightened by the fact that some months ago our President told 
you all about my work ; he did such a thorough job that I suggested that I sign my 
name to his paper and call it the George Darwin Lecture. My one criticism of his 
address is that he was forced to limit it to the successes; he could not give you a 
true picture of the troubles and failures of the work in photometry. 

I never had the privilege of meeting Sir George Darwin, but I knew very well 
his disciple, Ernest W. Brown. Brown once told me that the training he received 
from Darwin consisted simply in going to his study and talking with him about 
Brown’s problems. Likewise it was an education to talk with, or rather listen to, 
Ernest Brown. I met Sir James Jeans several times, the last when he was presiding 
here over a meeting of this Society. My first collaborator in photometry, F. C. 
Brown, had been a student of Jeans at Princeton ; and so in a way I might claim to 
be a descendant of both Darwin and Jeans, though not in direct line. 

The electrical photometry of stars involves the technique of experimental 
physics at the end of a telescope, and since all my work in this field has been in 
collaboration with others, it is only fair at the outset to acknowledge my debt to 
those without whose aid I should notbe here. There was F. C. Brown, as mentioned, 
who started me off for a few years with the selenium photometer; the late Jacob 
Kunz, who for twenty-five years produced photoelectric cells of the best quality 
with his own hands, and cooperated in other ways until his death a dozen years 
ago; C. M. Huffer, who has divided the observing time at the telescope for many 
years; A. E. Whitford, who first applied the thermionic amplifier, and continued 
with notable improvements and applications; and Gerald E. Kron, a former 
student with whom I am working, the relation of teacher and student being at least 
partially reversed. For the institutions involved I should mention twenty years at 
the University of Illinois, twenty-five years at the Washburn Observatory of the 
University of Wisconsin, fifteen summers at the Mount Wilson Observatory,.and 
the past two years at the Lick Observatory of the University of California, under the 
auspices of the United States Office of Naval Research. 

Instruments.—The present is a fine opportunity to open the floodgates of 
reminiscence, but I shall go back only to the last time I spoke in this room, during 
the meeting of the International Astronomical Union in 1925. After that meeting 
had dispersed I had the pleasure of going to Oxford and calling on Lord Cherwell, 
then Professor Lindemann, and I discussed with him the troubles we were having 
with a Lindemann electrometer which we had obtained for experiment a short time 
before. After a few minutes’ conversation he told me that the drift or instability 
we had found must be caused by a defect in the instrument. I returned to 
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Wisconsin, obtained another electrometer from the Cambridge Scientific Instru- 
ment Company and from the first night that this second electrometer was used on 
the telescope the previous installation was discarded, and Lindemann electro- 
meters were used in all photoelectric work at Madison and Mount Wilson for the 
next six or eight years. It would have paid my institution to send me across the 
Atlantic simply to learn that we had a defective sample of a very useful instrument! 

This and other forms of electrometer were finally supplanted, in our case, by 
the thermionic amplifier which Whitford successfully applied to the measurement 
of the minute current produced by the light of a star falling on a photosensitive 
surface, down to 5 x 10-!8 amperes. The series of photometers developed during 
the past twenty years represent the struggle to obtain greater effective sensitivity, 
either by finding a better photocell to increase the light signal or by reducing the 
irregularities or ‘‘ noise’ in the electrical circuit. Actually, more progress has 
been made in reducing the noise than in increasing the signal ; the advances can be 
illustrated by specific cases as we proceed. 

The final items in the improvement of a single-stage amplifier include mounting 
both photocell and amplifying tube in an evacuated chamber to improve surface 
insulation and to eliminate the effects of stray ions caused by cosmic rays, and 
refrigerating the cell down to the temperature of dry ice, — 80 deg. C., to reduce 
the thermionic dark current. All this is standard procedure in a laboratory but is 
somewhat out of ordinary at the end of a telescope. Also, whereas formerly it 
was our practice to make visual readings of a spot of light on a galvanometer scale, 
for several years past most of the American workers in this field have been using a 
recording potentiometer, which lightens the work at the telescope, enables one 
observer to operate alone, and gives a permanent record ; but like the chronograph 
for time observations, the recording instrument simply piles up more work for the 
daytime. 

Colour Systems.—If a photographic or visual observer publishes a set of 
stellar magnitudes or the light-curve of a variable star, his work may be taken as 
it is, but if the same task is accomplished with a photocell, then it seems that a 
determination of the colour curve or effective wave-length of the receiver is expected. 
This requirement is not onerous because experience has shown that there is little 
value in measuring the brightness of a star in fewer than two colours, so that 
comparison with any other system like that of the International can be made. 
Likewise, light-curves of a variable star in two colours always give more informa- 
tionthanasinglecurve. After one has once tried colour photography for snapshots 
with a hand camera the old-fashioned black-and-whites have little appeal. So for 
ten years or more we have followed multicolour photometry of stars and nebulae, 
always working in at least two colours, but more often using combinations of a cell 
and filters which isolate six spectral regions ranging from 3530 A. tu 10,300 A. (x), 
With filters the limits of such regions are not sharp; there may even be some over- 
lapping of adjacent regions, but for faint diffuse objects like the nebulae the saving 
of light makes up for this drawback. To isolate the spectral regions of an 
extragalactic nebula with a slit spectrograph would reduce the intensities below 
practicable measurement with a photocell. 

Another shortcoming of the method of filters is that the effective wave-length 
of the whole combination: star, atmosphere, telescope, filter, and receiver, 
depends upon the temperature or colour of the source. As we go along the spectral 
sequence from hot to cool stars the effective wave-length for any filter increases 
progressively, but, fortunately, for any two filters this change may be a nearly 
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parallel shift, the base line for colour index remaining of about the same length. 
Any method which integrates the intensities over all wave-lengths of a given 
interval also takes in the absorption lines and bands, and this condition holds for 
any system of colour index. 

The foregoing remarks apply with even greater force to the International 
photographic magnitudes and colours. As Seares (2) has shown, the photovisual 
system includes a range of about 1000 A. from 4900 to 5900 A., but the photographic 
range is about twice as great, from 3000 to 5000A. Thus, while the computed 
effective wave-lengths of Pv for black-body sources vary from 5420 A. at 30,000 deg. 
to 5480 A. at 3000 deg., the corresponding change in Pg is from 4160. to 4560A., a 
change of 400 A. as against 60a. for Pv. Similar changes in effective wave-length 
are present in all our six-colour work, but because of the parallel displacements 
for different filters, unless absolute measurements are involved, the relative 
colours for ordinary giant and dwarf stars are as useful as the results in any two- 
colour system. ‘The reason that such crude methods are used with the photocell 
and filters, or with the photographic plate, is of course to reach faint objects. 
The spectrophotometry of the Greenwich gradients, where definite regions of 
continuous spectrum can be selected for each star, ends at about the 5th magnitude 
with a 36-inch reflector. To go ten magnitudes fainter with any telescope requires 
a change of procedure. 

7 Aquilae.—As a sample of recent work in six colours we give the light-curves 
of the well-known variable star, 7 Aquilae, determined with the 36-inch Crossley 
reflector of the Lick Observatory. While in the similar curves for 8 Cephei (3), 
obtained with the 60-inch at Mount Wilson, we were able to smooth out most of 
the suspected secondary fluctuations—actually this smoothing had already 
been done by Dr Smart (4) with a photoelectric photometer at Cambridge—in 
n Aquilae the secondary “hump” still remains, and there are evidences of minor 
waves on the main curve, but not strong enough to be convincing. As in 5 Cephei, 
these new curves show decreasing amplitudes with increasing wave-length; the 
range of variation is about four times as great at 3530A. as at 10,300A. The 
overall colour class of 7 Aquilae varies through the cycle from about cf7 to cg 4; 
in 5 Cephei the range is from cf 4 tocg2. As might have been expected, the same 
retardation of phase with increasing wave-length is found; the maxima and 
minima come about % period later in the infra-red than in the ultra-violet. 

As pointed out by others, this effect is readily explained in a pulsating star by 
the combination of varying diameter and surface brightness of the apparent disk. 
At maximum light, when the radial-velocity measures indicate that the star is 
expanding, the increasing size of the disk is more than counterbalanced by the 
decreasing intensity in the short wave-lengths, and the star begins to fade; but 
in the long wave-lengths the same increase of area outweighs the much smaller 
decrease of intensity and the star still brightens. The reverse conditions obtain 
at minimum. As far as I am aware, no satisfactory complete theory of 7 Aquilae, 
accounting for the hump in the light-curve, has been worked out. Since 5 Cephei 
and 7 Aquilae are prototypes of their respective subdivisions of classical Cepheids, 
it is no great generalization to predict that the same phase effect is probably present 
in most other variables of these types, and also in the £ Geminorum stars. In an 
incomplete study we find this same effect in the cluster-type variable, RR Lyrae. 

Algol.—Another star studied in six colours, this time an eclipsing system, 
is our old friend Algol, but because of other problems the observational data were 
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Fic. 1.—Light-curves of n Aquilae. 





not quite completed in the season 1949-1950. Suffice it to say that the ranges 
of variation at 3530 A. and 10,300. are roughly 1-28 and 0-86 mag., respectively, 
for primary minimum; 0-03 and 0-13 mag. for secondary minimum; and 0-05 
and 0-08 mag. between minima; with proportional values for the intermediate 
wave-lengths. The data at hand would probably be sufficient for a solution for 
all the elements of the system, including the spectral type and relative size of the 
companion—the brighter component has a B8 spectrum—but it seems inadvisable 
to undertake such calculations before the data are definite, especially at the critical 
times of beginning and ending of the eclipses. 

M8 and Trifid Nebulae.—The six-colour measures of stars can often be 
compared with conventional two-colour results for the same objects, and we 
revert to some previous studies of reddened B stars near the centre of the galaxy. 
In the grouping of heavily obscured, strongly reddened B stars within 2° or 3° 
of the galactic equator, there is a gap at / = 335°, b= —2° (Lund pole), where a 
few stars at distances greater than 400 parsecs are of normal or near normal 
colour, the only such gap in 40° of longitude (5). It cannot be a mere coincidence 
that just in this gap are found the two conspicuous objects, M8 and the Trifid 
nebula. More likely is it that we see these two nebulae simply because there is 
a large hole in the intervening absorbing clouds, and no doubt many other such 
nebulae are present along and near the galactic equator, illuminated by the high- 
temperature B stars but obscured by the interstellar dust. Both M8 and the 
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Trifid are involved with reddened stars and it is the foreground that seems to be 
relatively clear. 

Globular Clusters.—However, not far north and south of these nebulae are the 
globular clusters, NGC 6440 and 6553, separated by the 7° zone of no globular 
clusters and reddened almost as strongly as any objects inthe sky. In Table I are 
some previously unpublished six colours of these and other globular clusters 
taken at Mount Wilson in 1940 to 1945. ‘The headings of the first three columns 
are self-explanatory. The fourth column contains the total photographic 
magnitudes by Christie (6). The values of C,, our first colours with the short 
base line, are from a previous paper (7). The six colours from U to I, referred 
in magnitude to the mean of blue, green and red, are on our standard system for 
stars. In the faint reddened clusters the violet and ultra-violet intensities were 
too weak to be measured on the background of the sky. Probably most of these 
clusters have varying integrated colour from centre to outside. We have found 
as much as 0-17 mag. difference in the V-I scale, inner redder than the outer 
region, but a study of this difference will involve a check upon the possible 
instrumental effects. 
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Fic. 2.—Colours of normal and reddened globular clusters. 











The data of Table I are shown graphically in Fig.2. The normal clusters 
have the curves to be expected from mixtures of stars, brighter in the U and I. 
A straight, horizontal line represents the mean colour of stars of spectrum dG6. 
Values for C,, drawn properly to scale, show that even with such a short base 











*L.,S wSeryderp {yout-og YIM sainseoUw [[Y 


60.0+- gl.o-+- 
pousppey Lo.o+- gl.0+ 
¥0.0+ go.1-+ 


¥o0.0+ go.0-++- ¥o0.0— Z1.0— 
Z0.0-+ go.o-+ 10.0— So.0 — 
[euLION ’ Z0.0-+- ¥o0.0-+- zZ1.0— ZZ.0— 
10.0+ go.0-+- 00.0 Lo.o— 
Z0.0-+- $0.0 -+- Z0.0— Lo.o— 


SL.1 So.z LE.z fg.c=y/t 
te lahat: | j Do a A nN ye Jaissayl 


s4ajsn] avpngojs fo sanojoy 


] FTaV,], 


3 
= 
ia~) 
2 
8 
¢ 
s 
o 
> 
> 
5 
Ss 
= 
a 
a 
S 
2 
& 
) 
& 
© 
= 
K 


No, 5, 1950 








422 Joel Stebbins Vol. 110 


line a good indication of the amount of reddening was possible. The globular 
clusters are good qualitative indicators of space reddening, but because of the 
dispersion in their intrinsic colours they are not as good as B stars for standards. 
As has been pointed out by Becker (8), the integrated spectral types of non- 
reddened globular clusters, determined by Mayall (9) at the Lick Observatory, 
show a strong correlation with the photoelectric colours, the difference between 
equivalent types A5 and G5 being about 0-10 mag. on the C, scale. Moreover, 
there is small chance of estimating the spectral types of the strongly reddened 
clusters. The difference between average O and Bz2 stars is only about 0-02 mag. 
on the same scale; and the dispersion within any subclass of the O’s and early B’s 
is very small indeed. 

Andromeda Nebula.—We have published elsewhere (10) the six-colour 
results on a number of bright extragalactic nebulae, but there are some differential 
measures in the Andromeda nebula which I reported at the meeting of the I.A.U. 
in Copenhagen in 1946. The observations were taken with the 60-inch reflector 
at Mount Wilson. With the same photocell as for six colours, the special system 
C, has two light filters giving effective wave-lengths at about 4250 A., near photo- 
graphic, and at 7500A., a base line about double that of the International system. 
The infra-red was shortened to eliminate the strong emission at 10,400 A. in the 
night sky. The photograph in Fig. 3 (Plate 7) has been retouched to bring out the 
fainter regions. The colour index for each area is referred to the mean of five 
such areas along the major axis; as usual the + sign means redder than the mean. 
The greatest difference, between + 0-30 in the dark lane and — 0-10 in the opposite 
region, is 0-40 mag. in C, or 0-20 mag., International. 

If we assume that the light of the nebula is symmetrical in these inner regions, 
except for intranebular absorption, then with the blue filter giving photographic 
magnitudes we have at once the ratio of total to selective absorption. To get more 
light, this determination was made with diaphragms twice the diameter of those in 
the figure. The relation between the scales of C,; and C,,, was found from the 
North Polar Sequence. If AA,, is the change in total photographic absorption 
from one region to another, and AC, the corresponding change in colour, then we 
have : 


AA,,,/AC, = 1-82 + 0:06 (p.e.), 
AC,/ACin, =2°13 + 0°05, 
MA jg/ACint = 3°87 + 0°16. 


Another result = 4-17 + 0°53. 


Mean AA,,,/ACint | 


=4°0 + 0°2. 
or Agg/ Eine 

The second result is from an independent but poorer determination made 
previously with other filters giving two-thirds the scale of C;. This ratio of about 
40 is all right for M31 if we place the dark lane on the near side of the nebula and 
permit the assumption of symmetry. We may note also that Holmberg (11) has 
found A,,,/ Ej =4°3 from a similar study of NGC 5195, the companion nebula to 
Ms1. Either value is not too different from various results found for the interstellar 
absorption in our galaxy, though I prefer a figure nearer 5 as determined from B 


stars. 
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Colours in M31. 
Retouched photograph by Mayall and Chappell, Lick Observatory. 
Diameter of circles : 2'-04. 
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Standard Lamp.—To calibrate the six-colour results on an absolute scale, 
Dr Kron and I have recently undertaken some comparisons of the stars with a 
standard lamp of known temperature and spectral distribution. The procedure 
was to measure with the Crossley reflector an artificial star at different distances 
on the peaks of Mount Hamilton. The star was produced by a hole 5 inch(0-3 mm.) 
in front of the enclosed lamp with tungsten ribbon filament, operated by storage 
batteries at a colour temperature near 2500 deg. K. as calibrated by the U.S. 
Bureau of Standards. At a distance of 1100 feet the hole subtended an angle of 
o”-2, and visually the star appeared red and between zero and Ist magnitude; 
brighter than Antares but fainter than Arcturus. With a wire screen absorbing 
five magnitudes, the light was reduced for convenient comparison with real stars 
of the 4th or 5th magnitude. 

For determination and elimination of the horizontal atmospheric extinction, 
the lamp, equipment and batteries were carried by jeep to Kepler Peak, 3200 feet 
from the Crossley ; at the greater distance a 2-5-mag. screen again equalized the 
artificial and real stars. Contrary to expectation the elimination of the extinction 
gave little difficulty. Tie most we could find for the lamp was a difference of 
0-02 mag. for 4220. between the longer and shorter distances, or 0-01 mag. for 
1100 feet. ‘The extinction for longer wave-lengths was negligible ; because of the 
red colour of the lamp we did not attempt to measure the weak ultra-violet at 
3530 A. 

Similarly, the regular atmospheric extinction was determined on the same 
nights by measuring stars at high and low altitudes. With an air mass, sec 32, 
varying from 1-2 to 3 or 4, the reduction to no atmosphere for each colour could be 
made with little uncertainty. All told, the accidental errors of the comparison 
between lamp and stars on different nights were satisfactorily small, though such 
measures would be difficult if not impossible with a sky less favourable than that of 
Mount Hamilton. 

The difficulty of the comparison turns out to be the uncertainty of the effective 
wave-lengths of the combinations of the photocell and filters, which are still in 
process of determination in the laboratory. We have found differences in colour 
sensitivity from spot to spot over the surface of the photocell. As used on the 
stars a positive lens in the focal plane forms an image of the main mirror about 
I cm. in diameter on the cathode of the cell, and the problem in the laboratory is 
to form such an image and measure its energy with a thermocouple. Until this 
calibration is accomplished, the observations cannot be reduced to an absolute 
standard. 

Magnitude of Sirius.—We turn for the moment from colours to stellar 
magnitudes, where again the photocell has the great advantage of a linear relation 
between the response and the intensity of the incident light. Actually the 
brightest star we have observed is the Sun, which at the 1937 eclipse in Peru was 
used as comparison for the corona, and to determine the extinction at the low 
altitude of 5° at totality. An optical “sun reducer’’, essentially a compound 
microscope, gave an out-of-focus image on the photocell of the same order of 
intensity as the corona or the full Moon. 

Of the regular stars we begin with Sirius, the brightest of all. Naturally it 
does not require a 36-inch telescope to measure Sirius, but we were interested in 
the colours, and incidentally we determined the visual magnitude, which is not too 
easy toget withtheeye. The light of Sirius was reduced with a 5-mag. wire screen, 
and the visual magnitude was interpolated between our blue and green magnitudes. 
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In TablelII are the results from four polar stars and five southern stars in the 
vicinity of Sirius. The magnitudes in the second column are Harvard visual. 
The spectra are for information only; the photoelectric colours were used in the 
reductions. Each visual magnitude in the fourth column is based upon that of the 
corresponding star. On any night a number of later type G or K stars were used to 
get the proper interval for interpolating between our blue and green, just as polar 
stars are used for collimation or azimuth in time determinations. The average 
deviations of the photoelectric measures on different nights indicate that the 
discordances in the column for Sirius must be in the magnitudes of the comparison 
stars. I see no way out of the conclusion that the generally used visual magnitude 
of —1-58 is too bright, and that — 1-42 is nearer right when referred either to the 
stars near the pole or to the other white stars in the vicinity of Sirius. Incidentally 
the effective wave-length of the Harvard visual system is found to be 5380. for the 
stars in the table. 

Tas_e II 
Visual Magnitude of Sirius 
Sirius 
m 

4°44 1°32 
5°86 1°34 

‘28 1°44 
6°33 “1°47 

‘99 “36 

“39 “42 
5°19 *52 

“30 é “45 
5°93 é ~ 1°48 

Four polar stars "39 +006 
Five southern stars ~1°45 +0°05 


m 





All —1°42 +0°03 
Harvard —1'°58 va 


Faint Magnitudes.—From the bright end of the magnitude scale we now go to 
the faint limit of our observations. The advent of the RCA 1P21 photomultiplier 
has introduced a new epoch into photoelectric photometry of the stars. With only 
a few precautions such a multiplier on a telescope will give satisfactory results over 
a wide range of magnitudes. With the refinement of refrigeration of the multiplier, 
and the use of a recording instrument, the effective sensitivity, or ratio of signal to 
noise, may be something like 10 to 50 times better than what was previously 
attainable with an ordinary photocell. 

A record has been taken by Whitford and H. L. Johnson of a star of photo- 
graphic magnitude 19-3 with the 100-inch telescope in 1949. Througha diaphragm 
of diameter 11 seconds of arc the light of the sky background was more than ten 
times that of the star; the zero for the dark reading had to be depressed off the 
scale. Also, the noise or irregularities at the tops of the deflections are caused 
mostly by the sky light. A smaller diaphragm might be used, but under ordinary 
conditions of seeing and guiding we prefer to take no chances of losing some of the 
light. In the present case the star was kept centred as in photography, with an 
auxiliary eyepiece and guiding star. The measures required about 30 minutes 
plus the time for finding the field, setting on the star, selecting a good guiding star, 
etc. ; so the rate of observing for these faint objects is ratherslow. Since the average 
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absorption of the yellow and blue filters is about 0-9 mag., the intensities in those 
colours were slightly less than for a clear magnitude of 20-0, but both magnitude 
and colour must be measured to compare with any other system. 

Someone said (I believe it was Paul Merrill) that the best way to take an infra- 
red photograph is simply to wait. After another year or so the manufacturers will 
bring out a new red-sensitive plate which will reduce the exposures to half or less, 
and the work will be made that much easier. In 1907 with a selenium cell on a 
12-inch telescope I managed to get a barely detectable response from the first- 
magnitude star, Aldebaran. By waiting 40 years I have seen useful measures go to 
magnitude 19 ‘with the 100-inch, and mere detection, say, to magnitude 21. I can 
mention this fact with due modesty because my own direct contribution toward 
the improvement of some 20 magnitudes has been nil. ‘The Radio Corporation of 
America developed the photomultiplier, and Whitford and Johnson perfected, 
installed and operated the photometer. In the very first trials with selenium, F. C. 
Brown and I got nothing at all from Jupiter, so we turned to the Moon and measured 
it withoutatelescope. Ina way, the improvement has been perhaps 23 magnitudes, 
which is a measure more of the crudeness of the first attempts than of the perfection 
of the latest achievement. The use of the 100-inch telescope under the California 
sky, accounting for 4 or 5 magnitudes, has been an important if not a major factor 
in the advance. 

Selected Areas.—In a paper soon to be published (12) we are giving a series of 
photoelectric magnitudes, reduced to the International system, of stars in Selected 
Areas 57, 61 and 68, with a comparison with the North Polar Sequence. The 
agreement with the polar stars, taken with the 60-inch, is very good from magnitude 
7 to 15, but with the NPS as standard the discordances from the Mount Wilson 
Catalogue (13) in all three areas are surprisingly large both in accidental errors and 
in scale. The details are given in the paper, and we seem to have exhausted the 
possibilities of finding systematic errors in our work; the internal accidental 
errors of the photoelectric magnitudes are negligible in the comparison. Ina word, 
on the basis of our results the zero points of the catalogue magnitudes in the three 
areas agree with NPS somewhere between magnitudes 9 and 14, different for each 
area, but the catalogue magnitudes at 18-5 are about 0-5 mag. too bright. From 
these increasing discordances toward the faint end of the magnitude scale we believe 
that any further conclusions based upon magnitudes in these Selected Areas 
should be postponed until the disagreement is straightened out. 

Hubble’s (14) surveys of nebulae with limits of Pg magnitude 18-5 to 21-0, 
based in part upon Selected Areas, are mostly fainter than our limit with the 
photocell, and we can say nothing about the scale from 19 to 21; but from com- 
parisons with SA68, Baade (15) has given corrections to Hubble’s magnitudes 
of stars in the region of M31, ranging from + 0-35 at 18-3 to + 0-63 at20-0. The 
exact relations between our magnitudes and those of Baade in SA68 are not yet 
available. I have heard that some philosophers say that the only time you can be 
sure you are talking sense is when you are counting. In cosmology you can make 
sense only if you know exactly when and where you have stopped counting. 

Nebulae.—We come to the observations of faint nebulae. In 1947, in coopera- 
tion with Messrs. Hubble and Baade, Whitford and I started, in addition to the 
work on Selected Areas, the measurement of more than one hundred nebulae 
for magnitude and colour with the reflectors at Mount Wilson. Most of these 
objects of various types were on a routine programme, but we also concentrated 
on a few clusters of nebulae where the red shifts had been determined by Humason. 
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The clusters were in Virgo, Coma Berenices, Corona Borealis and Bootes, the 
last a single nebula with a red shift giving a velocity of 38,g00km./sec. The 
results were published two years ago (16). By limiting the discussion to the E 
types, the dispersions in the intrinsic colours were minimized, and we came to the 
unexpected result that these faint nebulae were reddened more than would be 
expected from their red shifts. Some of the conclusions were :— 


1. The energy-curves of the nebulae are those to be expected from mixtures 
of stars; they show excess in both the infra-red and ultra-violet over that of a 
black body or of asingle star. The ultra-violet excess must be taken into account 
in evaluating the effect of the red shift. 

2. The E nebulae show a reddening proportional to the red shift. This fact 
may be taken as corroborative evidence for a linear velocity-distance relation. 
About 60 per cent of the reddening is a colour excess not explained as a consequence 
of the red shift. The figures of the Bootes nebulae are : colour excess, 0-52 mag. ; 
amount attributable to red shift 0-22 mag., amount unaccounted for, 0-30 mag. 

3. If the excess reddening is a distance effect caused by selective absorption 
by internebular material similar to the interstellar material in our galaxy, then, 
with a ratio of four to one for total to selective absorption, the photo- 
graphic absorption for the Bootes nebula becomes A,,,= 4 x 0-30 = 1-2 mag. for an 
uncorrected distance of 70,000,000 parsecs, corrected to 40,000,000 parsecs 
because of absorption. Any such absorption would mean a virtual impenetra- 
bility of space within only a few times the distance to the Bootes cluster. Also, 
with the run of the nebular counts, which Hubble has found to indicate nearly 
uniform distribution in transparent space, this reddening would mean a strong 
outward increase in density of distribution. 

4. An alternative to the absorption effect is a time effect, namely, that the E 
nebulae were formerly redder than now, and what we see is their progressive 
change in colour in the interval of 2 x 108 years from the Bootes nebula inward. 
The effect could be produced if the E nebulae formerly contained large numbers of 
high-luminosity red giant stars which have faded since the light left the distant 
ones. The corrections to the photometric distances on the time-effect. hypothesis 
would be less than for the space-reddening hypothesis and of the opposite sign. 

The foregoing alternatives of absorption effect and time effect summarize the 
situation as Whitford and I saw it three years ago after one season of observations. 
We were aware that a further test would be to measure the colours of spirals in 
the clusters of Corona and Bootes, if they could be found. The difficulties are 
that because of greater dispersion in their intrinsic colours the spirals are poorer 
standards than the elliptical ones, and with less concentration of light toward 
their centres their red shifts are more difficult if not impossible to determine. 

In 1949 Whitford measured the colours of seven nebulae of types Sb and Sc 
in the vicinity of the Corona cluster. The results have been referred to privately 
as indicating that the Sb’s exhibit the same excess as the E’s, but that the Sc’s 
have colours agreeing with the red shift of 22,000km./sec. One of the compli- 
cations is whether some of the objects should be classified as Sb or Sbe. Also, 
since the red shifts have not been measured for these spirals, any inferences from 
their colours should wait until they are shown really to have the velocity of the 
Corona cluster. 

The case for the Bootes cluster is different. Within the past month, in 
1950 September, Whitford has sent me the latest measures by himself and 
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Johnson of three additional E nebulae in Bootes, taken with the 100-inch telescope. 
All of magnitude 19-0 or fainter, they form a little group not far from the previously 
measured No. 1 nebula ; and while their red shifts are lacking, there can be little 
doubt that they belong to the cluster. The preliminary data are in Table III. 
Tas_e III 
Magnitudes and Colours of Nebulae in Bootes 
Pg, Cy A.D. of C, Obs. Remarks 
18°23 +1°36 +0°07 3 
19°33 +1°39 +0°10 3 I poor 
19°17 + 1°31 +0°05 2 
18°97 +1°I7 “ I 
Mean C, “+131 +0°07 
Standard, E nebula +o'84 
AC, obs. +0°47 
AC, comp. +0°22 
Excess +0°25 


Good 


The unexplained excess is reduced to 0-25 mag. from the 0-30 mag. given by the 
single No. 1 nebula. The agreement of the first three nebulae is too good to be 
true; and a discordance like that of No: 4 is to be expected from the dispersion in 
colour of the E’s in Virgo and elsewhere. ‘The evidence for the reddening of the 
nebulae is much strengthened by having four objects instead of one at the farthest 
limit; and whether the excess is 0:25 or 0-30 mag. makes little difference; we are 
back to where we were in the discussion. 

Now that the observed reddening of the Bootes nebula has been confirmed, 
the only way to reduce the infra-red colour excess would seem to be to revise our 
calculated effect of the red shifts. The procedure, depending upon the energy 
curves of the Sun and the elliptical nebula, M32, can be found in our original 
paper. To the criticism that far-reaching conclusions should not be based upon 
observational results on a single cluster of nebulae, the answer is that up to recently 
only one other velocity near 40,000 km./sec. has been measured, that of a nebula 
in Ursa Major No. 2. ‘Ihe assumption that the linear relation between red shift 
and distance holds to magnitude 21, involves an extrapolation of 3 magnitudes 
from the last measured red shift at magnitude 18; that is, an extrapolation by a 
factor of 4 in the distance and by a factor of more than 60 in the volume of space 
considered. ‘These uncertainties, combined with the uncertainties in the zero 
point and scale of the magnitudes indicate that the whole picture of the observable 
region of space can be subject to serious modification by further work on faint 
nebulae. 

In his studies of the distribution of the nebulae, Hubble (17) considered both 
the possibilities of space absorption and of the change of luminosity of the nebulae 
with time, but he decided that the effect of the red shifts was more important in 
modifying the calculated departures from uniform distribution. Whether this 
view is to be maintained or not, the assumption that the nebulae radiate like 
black bodies at 6000 deg. must be changed, and, as mentioned, the zero point and 
scale of magnitudes down to the limit of discussion will have to be checked beyond 
question. Lest there be any misunderstanding about Hubble’s attitude on these 
matters, I repeat that it was on his initiative and with his cooperation that we 
undertook the photoelectric observations of the nebulae herein described. 

We should like to apply the photocell to something constructive in the problems 
of cosmology. One possibility could be the testing of faint field nebulae for cases 
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of large colour excess. No matter what the interpretation of the reddening may 
be, if the linear relation between colour index and red shifts holds, the search for 
non-cluster nebulae with large red shifts could presumably be made more rapidly 
with the photometer than with the spectrograph; say something like one hour 
per object as against several hours or even all night for the same limit, with the 
probability of a negative result most of the time. Here again, however, we can 
refer such problems to photography where hundreds of objects are recorded 
simultaneously and where there should be little difficulty in searching on yellow 
or red photographs down to photographic magnitude 21, because the reddened 
nebulae will be more than a magnitude brighter on the photovisual scale. If such 
a search for very red nebulae were made photographically, the suspicious objects 
could be checked with a photocell. 

Conclusion.—In closing this account of some of the problems we have attempted 
to attack with the photoelectric cell, I am somewhat appalled, even embarrassed, 
by the range of the problems, as well as by the range of the brightnesses of the 
objects we have observed. From the Sun to the Moon to Sirius to a star of 
magnitude 19, is about 45 magnitudes, or a change of intensity of ro! to 1; and 
throughout this whole interval the constant effort has been to get a large value for 
the ratio, signal to noise, also expressible as the ratio of data to errors; or, if you 
like, the ratioof what you want to what you do not want. During the improvement 
in the sensitivity of the apparatus by a factor of some hundreds of millions, the 
increase in relative precision has been something like a factor of only ten, from 
errors of a few hundredths to errors of a few thousandths of a magnitude; but the 
achievement of another decimal place is always welcome in any division of physical 
science. From my own considerable experience in the past, perhaps I should 
have ideas about some of the next steps in astronomical photometry; but I can- 
not predict what will happen in this field even in the near future. In any event 
it now seems truer than ever that in astronomy as elsewhere the applications of 
the photoelectric cell are limited only by the imagination of the experimenter. 


The recent and current investigations at the Lick Observatory, University of 
California, have been carried on as a project under the United States Office of 
Naval Research, and have also been supported by grants from the American 
Philosophical Society. 
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A PHOTOGRAPHIC SURVEY 
OF BRIGHT SOUTHERN PLANETARY NEBULAE 


David S. Evans and A. D. Thackeray 
(Received 1950 June 5) 


Summary 
Twenty-six objects believed to be planetary nebulae with diameters 
greater than 8 seconds of arc, south of declination — 40°, have been investigated. 
Descriptions, and photographs of those which definitely are planetary nebulae, 
are presented. Some comments on the classification of planetary nebulae are 


offered, and the proportions exhibiting various types of symmetrical form are 
tabulated. 





The present project originated in a suggestion by Professor J. H. Oort, made 
during a visit to the Radcliffe Observatory in 1948. It was felt that, although the 
proper motions of planetary nebulae would probably be very small, it was 
nevertheless desirable that plates which might provide first-epoch positions 
should be obtained as soon as possible. The project was later extended, and 
much longer exposure plates were taken, in order to provide, in addition, a 
descriptive catalogue of bright southern planetaries. A list of objects believed 
to be planetary nebulae, having diameters of more than 8 seconds of arc, south 
of declination — 40°, contained 26 objects. In the course of 1948-50, a total of 
123 plates of these objects was obtained. This included 16 plates of IC 4406 
(already described elsewhere by one of us, DSE)* and of the remainder, 89 
plates are by DSE and 18 by ADT. 

Table I lists the objects observed, with references to the literature, the 
numbers of the plates in the Radcliffe files which are recommended for possible 
proper motion studies, and the numbers of the accompanying illustrative 
photographs. 

At first sight it appears that there is a large literature on the planetary nebulae 
described here, but, in actual fact, a great many of the works cited are compilations 
from earlier publications. In the NGC (1), based to a considerable extent on the 
work of the Herschels, the classification of southern planetaries as such is by 
appearance, and rests on Sir John Herschel’s “Results of Astronomical 
Observations made in 1834-8 at the Cape of Good Hope” (London, 1847). 
Campbell and Moore’s work (2), done with the spectroscope, represents the most 
reliable source. A number of identifications by Mrs Fleming, quoted in the 
Index Catalogue, also have a spectroscopic basis. For the southern nebulae the 
source of the magnitudes and diameters given by Vorontsov-Velyaminov and 
Parenago (3) is the examination of the Franklin-Adams plates. The results 
obtained are repeated by Vorontsov-Velyaminov in his compilation (4). In this 
work it is stated that all the images were stellar, and that the appearance of the 
stellar images was used to assign integrated magnitudes. In addition the image 
diameters were used to assign diameters to the nebulae. This seems to be a 
double use of the same material, and does not seem very reliable. Indeed, 


* M.N., 110, 37, 1950. 
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Vorontsov-Velyaminov himself does not regard the results of any of these 
southern nebulae as meriting more than assignation to his class of third quality 
results. In addition all ‘stellar’? nebular images were arbitrarily assigned a 
diameter of ?2”. Most of the other references are to particular individual 
nebulae. It follows, therefore, that in fact the observational basis for previous 


TABLE I 


s eee Galactic PM Figure Re- 
. Designation Position (1950) pr Literature siees Me. tke 


ies s*49"4 - 244° —32° (7),(9),(xz) —... 1 PNP 
IC 2448 9° 061-6 —69 253° —15° (x), (2),(3), A648 2 P 
(4) A684 
NGC 2792 101-6 —42° 234° + 5° (xz), (2), (3), A649 P 
(4), (8), (14) A665 
NGC 2867 20M-qg —s58° 246° — (x), (2), (3), A685 P 
(4), (48) 
NGC 3132 10!04™-9 , —40 240° + (x), (2), (3), A252 P 
(4), (12), (13) A253 
NGC 3195 10" 10™- -80° 37’ 264° — (1) A265 
NGC 3211 10% 16™-2 —62° 26’ — 5° (x), (2), (3), Azo1 
(4) 
rr 26m-2) — f 458° 4 (2), (4) A297 
A298 
NGC 3918 11" 47-8 ’ 263° + 5° (x), (2),(3), A704 
(4) 
NGC 5189 1329-9 —6s5 : - (x), (2), (4), 
(12) 
NGC 5307 13% 47™- 5 > (x), (2), (3), 
(4) 


IC 4406 14" 19™- 43°55 
va 15) 30m. 59’ 7 (3), (4), (5) 


15h 47m. , (6) 


16) ;oM-5  — 50° (4), (5) 


161 73M-5 5 (4), (5) 


NGC 6153 165 28m. > + (1), (2), (3), 
(4), (10), (11) 
IC 4642 172 97M. , 302° (x), (2), (3), 
(4), (6) 
NGC 6326 175 16™-8 - ; 306° - (x), (2), (3), 
(4) 
ca 182 92.8 50 310 > (6) 
NGC 6630 18227M-7 —63 299° - (rz), (6) 
IC 4723 182 31M-; —63 299 (1), (6) 
NGC 6707 18% 51M-3 53° 310 (x), (17), (19) 
NGC 6708 Nf 6707 by 6’ 
NGC 6902 20%21™-2 —43° 325° — (x), (15), (16), 
(17) 
eee 22h 35™-5 —57° 29’ 298° —53° (3), (4), (5) 
NGC 7408 22% 52™-7 —63° 58’ 289° —49° (x), (3), (4), 
(5) 
In the Remarks column: P denotes an accepted planetary. 
NP denotes “‘ not a planetary ”’. 
PNP denotes “ probably not a planetary ”’. 
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discussions of southern planetary nebulae has been very slight, and is by no means 
as extensive as might be inferred from our reference list. 

We now proceed to a short discussion of each object in turn :— 

No. I: 

Innes (9), studying three plates taken with the 10-inch Franklin-Adams cameras, 
gives the three descriptions :— 

“Large but faint spiral with stars: a small cluster forms the nucleus.” 

“ Fine circular spiral with stars.” 

“Very fine circular spiral with stellar nucleus.” 

In 1927 the same author (11) says: ‘‘ The stars in this spiral nebula are shown 
as a letter S but the nebulosity is virtually invisible’. Jenka Mohr (7) discusses 
plates taken with the 24-inch Bruce, the 10-inch Metcalf and the 3-inch Ross 
Fecker instruments at Bloemfontein. The diameter is given as 420 seconds of arc, 
and it is stated that the nucleus “‘ appears to be a central star of magnitude 12-6, 
but it is possible that this star is merely superposed”’. It is pointed out that 
the nebula is almost certainly superposed on the Magellanic Cloud, since, at the 
distance of the cloud, it would have the very large diameter of 50 parsecs. Estimates 
based on statistical. discussions of planetaries, assuming the selected star to be the 
central star, give distances of 1200 parsecs and 2000 parsecs. 

Our plates show very faint nebulosity (invisible to the eye in the telescope) 
with a diameter of 420 seconds of arc. The nebulosity is much brighter on the 
preceding side. There is a wispy structure with some suggestion of a spiral or 
helical form. The nebula does not seem to be a typical planetary, and it is quite 
impossible to designate any particular star as the nucleus. Our photographs 
show the central area to be occupied by two distinct and compact groups of stars, 
some of which appear to be directly mingled with nebulous wisps. The galactic 
latitude is greater than that of any accepted planetary in our list, and in view of 
the great variety of nebulosities with complex forms appearing throughout the 
Magellanic Cloud, the explanation in terms of a galactic planetary superposed 
on the Cloud is not clearly demanded. Evidence in favour of its belonging to 
the Cloud is strengthened by the presence of the double cluster in the central 
region, which is unlikely to be galactic. If it belongs to the Cloud it must be 
rejected as a true planetary, and the problem arises of how such regularity can 
be exhibited by nebulosity 50 parsecs in diameter. Purely on appearance the 
nebula might be thought to resemble the final stage of a supernova outburst. 
The photograph, Plate 8, Fig. 1, has been intensified to bring up the faint 
nebulosity. 

No. 2: 1C 2448. 


This planetary (1) was found on the basis of its spectrum by Mrs Fleming, who 


, 


designated it as ‘‘ Planetary, stellar’. Campbell and Moore (2) give a diameter 
of 8”, magnitude 9, spectrum N,, N,, H8 (10:4:1), velocity = —24 km./sec. 
Vorontsov-Velyaminov (3), (4) gives diameter 8”, integrated magnitude 11-5. 
Our plates, Fig. 2, show a bright central star, with a small oval patch of luminosity 
placed slightly eccentrically ; dimensions 9” x 7”. 
No. 3: NGC 2792. 

Dreyer (1) says: ‘Remarkable planetary, pretty bright, equivalent to 9th 
magnitude star, very small, round, among stars”. Campbell and Moore (2) 
give a diameter of 10”, 10th magnitude, spectrum N,, Ne, velocity = +14. 


31 
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Innes (8) says: ‘‘ An 11th magnitude planetary about 20” in diameter. Is Np a 
pair of 10-5 mag. stars. No stars within 3’”. Gregory (14) has: “ Planetary, 
pretty bright, small, very little extended in P.A. 165°+. Appears as a nearly 
uniform disc without a central star”. Vorontsov-Velyaminov (3), (4) assigns 
diameter 10", magnitude 13-5. We find a typical, almost symmetrical, ring 
nebula, with central star. The ring is slightly brighter to the south. Diameter 
13”. (Plate 8, Fig. 3.) 

No. 4: NGC 2867. 


Dreyer (1) says: “‘ Very remarkable planetary, equals 5th mag. star, very small, 
round. 15th mag. star at P.A. 59°, distant 13””’. Delisle Stewart (18) gives 
a diameter of 18”. Campbell and Moore (2) have: “‘ Diameter 8”, round, edges 
sharply defined. Centre appears fainter. May be very small ring nebula. 
Mag. 8-5. Velocity=+18 km./sec.”. Vorontsov-Velyaminov (3), (4) has 
diameter 8”, magnitude 9-7. The object is a bean-shaped ring nebula, the ring 
being broken up into spots of brightness, brightest in the Sp sector. The 
central star is faint. About 13”x11". (Plate 8, Fig. 4.) 
No. 5: NGC 3132. 


Dreyer (1) says: ‘“‘ Very remarkable planetary; very bright, very large; little 
extended; star of gth magnitude in the middle, 48 diameter”. Campbell and 
Moore say: “‘ Ring nebula with 8-8 magnitude star in centre. Diameter of ring 
about 30”. Spectrum: N,, N,, HB, N,(10:4:1:1); V=—8km./sec.”. There 
is a drawing by Wilson. Vorontsov-Velyaminov (3), (4) gives diameter 30”, 
integrated magnitude 8-2, central star, magnitude 10-6. Knox Shaw (13) has: 
“Very bright, pretty small, little elongated in P.A. 150° + : annular ring, fainter 
on preceding side. B star not quite central and thus probably unconnected”’. 
A note says that the spectrum of the bright central star is continuous. Shapley 
and Paraskevopoulos (12) say: ‘“‘A series of photographs of varying exposures 
would be necessary to bring out the intricate details . . . . It could well be named 
the ‘‘8-burst”’ planetary from the number of distinct arcs on the boundary of 
the main disc or shell. The class A star, HD87892, mag. 9:5 is centrally 
superposed”. A photograph is given. 

The object has a complex structure, including a bean-shaped ring, very bright, 
with fainter, helical, extensions. Dimensions about 84” x53”. The object 
is illustrated in Plate 8, Fig. 5. Photographic manipulation has been employed 
to bring up the centre on the print, but the structure has not been falsified. There 
is a nucleus or foreground star on the ring in the Sp sector. 





LEGENDS TO PLATE 8. 
Fic. 1. 
Object: 5h 43™-4, —67° 53’ 
x 2:7. 
Kodak 103-a O, 30 mins. 


Fic. 2. Fic. 3. Fic. 4. 
IC 2448, X7-0. NGC 2792, X7°0. NGC 2867, 7-0. 
Ilford Astronomical Ilford Astronomical Ilford Astronomical 
Zenith, go secs. Zenith, 5 mins. Zenith, 60 secs. 
Fic. 5. Fic. 6. Fic. 7. 
NGC 3132, x5°8. NGC 3195,  6°3. NGC 3211, x 6-3. 
Kodak 103-a O, 16 mins. Kodak 103-a O, 15 mins. Ilford Astr 7 
Zenith, 10 mins. 
Note.—The magnifications are relative to the original plate scale of 22”:5/mm. 
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Fic, 6, Fic. .7 


David S. Evans and A. D, Thackeray, A Photographic Survey of Bright Southern Planetary Nebulae. 





Montuty Notices or R.A.S. Vor. 110, Prate 9. 





Fic. 8, ' Fic. 9. 
Object: 11h 26™-2, —52° 39’ NGC 3918, x5. 

<9°0. Ilford Astronomical 
Kodak 103-a O, 30 mins. Zenith, 20 secs. 


Fic. 10. Fic. 11. 
NGC 5189, 7:0. NGC 5307, 7:0. 
Kodak 103-a O, 30 mins. Ilford Astronomical 
Zenith, 5 mins. 
David S. Evans and A. D. Thackeray, A Photographic Survey of Bright Southern Planetary Nebulae. 
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Fic. 12. Fic. 13. 
Object: 15% 47™-4, —51° 21’ Object: 15% 30™-2, —58° 50’ 
“8-8. x 6°3. 
Kodak 103-a E, 20 mins. Kodak 103-a E, 20 mins. 





Fic. 14. Fic. 14 5. 
Object: 16% 10M-5, —54° 50’ Object: 164 1om-5, —54° 50’ 
x 10. X 23. 
Kodak 103-a O, 40 mins. Kodak 103-a O, 40 mins. 


David S. Evans and A. D. Thackeray, A Photographic Survey of Bright Southern Planetary Nebulae. 
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Fic. 21. 


David S. Evans and A. D. Thackeray, A Photographic Survey of Bright Southern Planetary Nebulae. 
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No. 6: NGC 3195. 
Dreyer (1) says: “‘ Remarkable planetary, pretty bright, small, little extended; 
13° diameter; three small stars near’. We find (Plate 8, Fig. 6) an oval nebula 
with central star. There are two sectors of great brightness on the oval to east 
and west, and faint extensions in the form of a ring to north and south. Dimen- 
sions about 44” x 33”. 
No. 7: NGC 3211. 


Dreyer (1) says: “‘ Planetary, equals roth magnitude star, round, among 150 stars’’. 
Campbell and Moore (2) give: “ Planetary, round, diameter 8”, roth magnitude. 
Spectrum N,, N,, HB, Ng (10:4:1:2). V=—16 km./sec.”. Vorontsov- 
Velyaminov (3), (4) has diameter 8”, integrated magnitude 11-8. The object 
(Plate 8, Fig. 7) is a bright disk 14” in diameter, with a slight bump on the 
northern edge. The central star is barely visible. 

No. 8: 
Campbell and Moore (2) give the description: ‘‘ Faint patch of nebulosity slightly 
elongated: trace of nucleus: diameter about 15”, but limits are quite indefinite. 
Spectrum; N,; V=+2-8 km./sec. All the plates are much underexposed’’. 
Vorontsov-Velyaminov (4) merely repeats these details. 

Positive prints show an image precisely like the face of a circus clown, but 
a negative is more revealing. ‘The central star is surrounded by an oval extended 
roughly east and west, and having dimensions about 45” x 30”. There is a faint 
dark lane concentric with the main oval, so that in fact we have an oval outer ring. 
On the following side is a dark indentation (the clown’s mouth), while the “eyes” 
are dark areas extending almost north and south from the central star. The ring 
is brighter on the north preceding side, while on the preceding end is a bright 
nucleus (? star). There is a very faint equatorial hoop extending north and 
south (the clown’s “ears’’), having a total extent of about 60”. The contrast 
of the negative print (Plate 9, Fig. 8) has been increased. 


No. 9: NGC 3918. 


Dreyer (1) says: ‘‘ Planetary, remarkable, small, round, blue, equals 7th magnitude 
star. Diameter 18-5’. Campbell and Moore (2) say: “‘ Very bright planetary, 
magnitude 8-5, diameter 10”, spectrum :—N,, N,, HB, N3, 44363, Hy, Hb, 43768; 
V =—16-4 km./sec. Intensities, N,, N,, HB, Nz (10:3:1:1)”. Vorontsov- 
Velyaminov (3), (4) gives a diameter of 10” and integrated magnitude 8-4. We 
find a simple disk with very slight equatorial structure. No central star is visible. 





LEGENDS TO PLATE 11. 
Fic. 15. Fic. 16. Fic. 17. 
Object : 16h 13™-+5, —51° 52’ NGC 6153, X9°0. IC 4642, x65. 
x 6°5. Ilford Astronomical Ilford Astronomical 
Kodak O-E, 20 mins. Zenith, 10 mins. Zenith, 10 mins. 
Fic. 18. Fic. 19. Fic. 20. 
NGC 6326, x 6°6. NGC 6630, x 6°6. IC 4723, < 6°6. 
Kodak O-O, 3 mins. Kodak 103-a E, 15 mins. Kodak 103~a E, 15 mins. 


Fic. 21. 
Types of symmetrical planetary nebulae. 
a b c d 
e f g h 
j R l m 
Note.—The magnifications are relative to the original plate scale of 22”-5/mm. 


37° 
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Efforts by one of us (ADT) to reveal it, using green filters and F-type plates to 
cut out the nebulium lines, have met with no success. Diameter 13”. (Plate 9, 
Fig. 9.) 
No. 10: NGC 5189. 

The NGC says: “Remarkable, blue, pretty large: considerably extended: 
brighter towards the middle: curved axis: 4 stars involved”. There is a 
drawing by Sir John Herschel (loc. cit.). Campbell and Moore (2) say: ‘‘ Nebulous 
field about 1’-3inlength. Composed of three nebulae. NGC 51809, the following 
one of the three is about 20” in diameter, and slightly fainter than the spindle shaped 
nebula in the centre. Spectrum, N,, N, (10:3). V=—6 km./sec.”. There 
is a drawing by Wilson. It is evident from this that Campbell and Moore regard 
NGC 5189 as being one of a group of nebulae, whereas in fact the whole mass is 
a single complex nebula which is not a typical planetary. Vorontsov-Velyaminov 
(4) simply quotes Campbell and Moore (2). Shapley and Paraskevopoulos (12) 
give a drawing and say: ‘‘ NGC 5189, near the galactic equator is a gaseous nebula 
of such remarkable knotted structure that it is here best represented by a drawing. 
No bright star is involved in the nebulosity, but 6’ to the south is the 7th magnitude 
star HD117694, Bo. In the HD catalogue there is no star with a spectrum 
earlier than B5 within a degree of the nebulosity. Overall dimensions 3'-0 x 2’-0”’. 
Our photograph, Plate 9, Fig. 10, shows no central star; the object is clearly 
not an ordinary planetary. Dimension, east-west 185”; north-south 130”. The 
structure is helical or S-shaped*. 


No. 11: NGC 5307. 
Dreyer (1) says: “‘ Planetary, or very faint extended small double nebula’’. 


Campbell and Moore (2) say: “‘ Planetary, elongated in north and south direction, 
8” x10". Magnituder1. Spectrum: N,, N,, HB (10:3: ?). V=+42km./sec.”’. 
Vorontsov-Velyaminov (3), (4) gives the same dimensions, but integrated 
magnitude 12-1. We find (Plate 9, Fig. 11) a small gaseous parallelogram 
elongated roughly north and south with a central star. Dimensions about 
15" x10”. 


No. 12: IC 4406. 


This nebula has been discussed at length by one of us (DSE).t 
No. 13: 

Found by Menzel (§) ona Bruce plate taken at Arequipa and classified as planetary 
on the basis of its appearance. Menzel states that it is bi-nuclear and gives an 
angular diameter of 24”. Vorontsov-Velyaminov (3), (4) gives a diameter of 24”, 
integrated magnitude 12-5. Our photograph, Plate 10, Fig. 13, shows a very 
complex structure. A bright ring 34” in diameter and of oval shape, shows two 
stars within it. The appearance is as if this bright ring formed the base of a 
cylinder with generators extending roughly northward, with the exception of one 
line, which also extends slightly south. The nebula terminates in the north with 
a second, fainter and slightly smaller ring, seeming to form the top of the cylinder. 
This ring has three bright spots on it, and also has a bright nucleus ( ? faint star). 
Total length about 55”. 

* NGC 5189 is referred to by Radlova, Katz and Dokychaeva, Russian Astronomical Fournal, 
26, 160, 1949 (in Russian). The data are however merely derived from Franklin-Adams and Wolf 
Palisa charts. 

t M.N., 110, 37, 195°. 
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No. 14: 


Shapley (6), who found this fine ring nebula on plates of Milky Way field 233, 
gives a magnitude of 3-16 and a diameter of 72”. He says: “‘ The width of the 
ring, which is circular, and nearly uniform all the way round, is 15 seconds of arc”’. 
Plate 10, Fig. 12, in which the contrast of the faint ring has been increased, shows 
these data to be correct. The ring is brighter on the preceding side. 
No. 15: 

Found by Menzel (5) and classified on appearance. Diameter 30”, quoted by 
Vorontsov-Velyaminov (4). Menzel says that it appears to be annular. The 
object consists of a ring of diameter 28”, breached to north and south, so that 
in effect, it is composed of two comma-shaped masses of gas. From the tails of the 
commas extend curved wisps of faintly illuminated gas with small radii of 
curvature. The northern one has a bright knot on it, almost in line with the 
breach in the main ring. Total extent, north to south, about 45”. Plate ro, 
Fig. 14a is anormal print: the contrast of Plate 10, Fig. 146 has been increased, 
the better to show the fainter parts of the nebula. 


No. 16: 


Found by Menzel (§), who gives a diameter of 12” and says that it appears to be 
annular. Vorontsov-Velyaminov (4) quotes this diameter. This peculiar 
object consists of a central star with luminous masses to north and south (Plate 11, 
Fig. 15). That to the north is almost globular and appears slightly to overlap the 
central star. ‘The mass to the south is concave towards the star, and shows some 
trace of an absorption shell. Extending southward from this second mass is a 
faintly luminous pointed “beard”. Dimensions: north-south, 35”; east-west 
about 16”. 

No. 17: NGC6153. 

Dreyer (1) describes this planetary, discovered by Copeland, as “‘ planetary, 
stellar’’. Campbell and Moore (2) observed the N, line, and assigned V = + 39. 
Van den Bos (11), observing with the 26-inch refractor of the Union Observatory, 
described it as a planetary, with major axis 28” in P.A. 7°-187°, minor axis 22”. 
Finsen (10) describes it as a ring nebula. Vorontsov-Velyaminov (3), (4) 
assigns his conventional diameter ?2”, and magnitude I1°5. 

We find (Plate 11, Fig. 16) an object with a central star surrounded by acircular 
area of luminosity of diameter 20”. On this is superposed a curious arched 
structure resembling a geometrical diagram. It is a line of gas in the form of 
part of an elongated ellipse terminated by an asymmetrical chord. The curved 
part, which extends beyond the circular area, is faint. The sharper (preceding 
side) join of the chord and the ellipse is slightly extended in the form of a faint 
point projecting a little beyond the circular patch of luminosity. Dimensions: 
east-west, 21”; north-south, 28”. 

No. 18: IC 4642. 

Dreyer (1) takes this from the list of planetaries discovered by Mrs Fleming 
from their spectra. It is “planetary, stellar”. Campbell and Moore (2) say: 
‘“‘ This nebula is remarkable for the strength of the N, radiations. ‘The exposures 
are too short. V = +44 km./sec.”’. 

Vorontsov-Velyaminov (3), (4) assigns the conventional diameter ?2”, and 
integrated magnitude 12-4. Shapley (6) gives magnitude 14-4 and diameter 30”. 
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We find (Plate 11, Fig .17) an oval ring nebula with central star. The principal 
axis is in the direction Nf-Sp. The illumination of the ring is not uniform, 
and dark areas form a cross on the major and minor axes of the oval. Dimensions 
18” x15”. 

No. 19: NGC 6326. 


Dreyer (1) says: “‘ Very, very remarkable: pretty bright: very small, round”. 
Herschel (loc. cit.) has a drawing. Campbell and Moore found the lines N,, N,, 
HB (10:3:1+) (HB very very faint), and give V= +11 km./sec. Vorontsov- 
Velyaminov (3), (4) assigns a magnitude 12-2 but no diameter. The object 
(Plate 11, Fig. 18) consists of a roughly rectangular area of luminosity about 
12” x 10” with a central star. This forms the main “ring’’ of the nebula and 
in the south following corner it has a bright spot. There is a smaller, fainter, 
rectangular extension from the southern half of the preceding side. Dimen- 
sions about 16” x 11”. 


No. 20: 
Shapley (6) reports a nebula, magnitude 11-9, diameter 36” at position (1900) 
17" 58™ 538, —50°37’:1. We find nothing at the corresponding 1950 position. 
No. 21: NGC 6630. 


Dreyer (1) says: “Planetary, faint, small, round, gradually brighter in the 
middle”. Shapley (6) gives magnitude 15-3, diameter 30”. The classification 
is by appearance. We find an object with a quadrilateral of four stars or nuclei 
surrounded by a nebulous envelope. The object is difficult, and it seems open 
to doubt whether the object is correctly classified as a planetary nebula. 
Dimensions about 19” x15". (Plate 11, Fig. 19.) 

No. 22: 1C 4723. 

The object enters the Index Catalogue (1) as one classified as a planetary on 
appearance by Dudley Stewart. His description is: ‘very faint, very small, 
round”. Shapley (6) gives magnitude 15-4, diameter 27”. 

We find (Plate 11, Fig. 20) a curious round object of diameter about 19” with 
three distinct nuclei running along a north-south line across it. The designation 
as a planetary again seems doubtful. 

No. 23: NGC6707, 6708. 

For the former, Dreyer (1) has: “ Faint, small, very little elongated, gradually 
brighter in the middle”. For the second: “ Pretty faint, small, round, gradually 
pretty much brighter in the middle, last of group”’. Reynolds (17) has NGC 6707 
as “Faint spiral, 3'x1’'”. The classification as planetaries is by Frost (19): 
“very faint planetary elongated about 2’”, and “ Round planetary, magnitude 
14”. The objects prove to be definitely extragalactic. 

No. 24: NGC 6902. 

Dreyer (1)says: “‘ Faint, considerably small, round, brighter in middle”. Shapley 
and Ames (16) make it extragalactic: they give no description. Reynolds (17) 
also calls it a spiral, without description. The classification as a planetary 
nebula is by Madwar (15), who says: “ Faint, round, 30” diameter; planetary 
with a bright, almost stellar, nucleus”’. 

Short exposures do show a structure closely resembling a planetary nebula, 
but longer ones reveal spiral arms. These are very faint, but the object is 
extragalactic, presumably a late-type spiral. 
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No. 25: 


Discovered by Menzel (5), who gives a diameter of 30". Vorontsov-Velyaminov 
(3), (4) makes it a planetary fainter than the 15th magnitude with diameter 30”. 
Our plates show a faint ring structure broken up into bright spots surrounding 
a blurred nucleus. The field is full of extragalactic nebulae, and we regard the 
classification as a planetary as incorrect. 


No. 26: NGC 7408. 


The NGC says: “ Pretty bright, pretty small, round, very gradually a little brighter 
in middle’. Vorontsov-Velyaminov gives a diameter of 72” and a magnitude 
of 12-8. Menzel (5) assigns it to the planetary nebulae on the basis of appearance 
and assigns a diameter of 72”. The object is a barred spiral with a line of 
condensations or giant stars crossing it. 





It is not the purpose of the present paper to offer a theoretical analysis of the 
forms of these planetaries, but we are undoubtedly in a somewhat favoured 
position in having the original material to hand, and there would seem to be no 
harm in mentioning some of the theoretical points—avowedly speculative—which 
we have discussed between ourselves. It is possible that these may suggest 
sound lines of investigation to theoretical workers. 

A first problem is the question of the stability of ring structures: the classical 
form of the planetary nebula is a star surrounded by a uniform shell such as 
No. 14. Is this form stable, or if not, will the ring later be breached at roughly 
opposite points as in No. 15? 

A second problem is provided by the frequency of axially symmetrical forms 
such as Nos. 6, 8, 16 and 18. This type of form would seem to link up with 
the fairly common helical form such as Nos. 5 and 13, where the impression is, 
for all the world, as if matter had been squirted out from two opposite poles of a 
rotating star and left behind as it went along. This crude picture is palpably 
ridiculous: the ejected matter must have the momentum of the parent star. 
Can it then be that this type of form results from the interaction of two bodies 
or of the star with an interstellar medium? No. 13 is designated in the literature 
by ‘‘? two nuclei’’, while No. 5 has a second star or nucleus superposed on the 
nebula or involved in it. We can say nothing more here on this point. 

Finally, the form of No. 15 is especially intriguing. One is almost compelled 
to see the structure as one in which matter is being driven away through the two 
breaches by radiation pressure. Then when it gets clear of the inner ring it 
becomes less luminous, since it is now partly screened by the central ring. 
Possibly the bright knot to the north represents a region crossing back across the 
full stream of stellar radiation escaping through the northern breach. But why, 
we may ask ourselves, does this wispy material show so marked a curvature? 
The track is impossible for matter moving under a purely gravitational force, 
and radiation pressure would reduce the curvature, not increase it. Is it merely 
fanciful to see a resemblance to solar prominences in these arched structures? 
Qualitatively one might see the central star, plus the inner ring, as a region from 
which energetic electrons might readily escape to form a space charge around the 
structure, while the wisps of gas themselves must be regions in which, outside 
the main radiation blast of the central star, recombination is the dominant process : 
that is to say, the wisp must be rich in ions. Have we here a situation somewhat 
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similar to that proposed by one of us (DSE)* in connection with prominences, 
where, it was suggested, ions move in electric or magnetic fields? If this were so 
it might account for the opposite senses of the curvatures of the northern and 
southern wisps. 

Now that we have good photographs of planetary nebulae available all round 
the Milky Way, it may be useful to offer some comments on the subject of the 
classification of planetary nebulae. Curtist proposed a classification on the 
basis of form, but, as has been pointed out by Stoy {, this classification is not very 
satisfactory. It does not satisfy the essential conditions that the system must 
be so devised that each planetary goes into one and only one class, and that 
independent workers, classifying the same material, must arrive at substantially 
the same results. A classification on the basis of the observed spectrum is also 
possible, and has been proposed by Wright.§ However, it seems desirable that 
both types of classification should be attempted, since the relations between the 
two may throw light eventually on the origin and evolution of planetary nebulae. 

One of the most striking things about planetary nebulae is that a high 
proportion of them possess symmetry of various types. Any system which has 
symmetry properties immediately becomes much more amenable to the 
mathematics of theoretical discussion, and it therefore seems desirable in order 
to assist theoretical development to concentrate on the symmetrical forms which 
occur. ‘This may be only a temporary expedient, but, as the result of trial and 
error, the authors have found that the use of symmetry properties as a basis for 
classification is one (and indeed the only one) which enables a substantial measure 
of agreement to be secured between independent workers attempting to classify 
the same material. 

In Plate 11, Fig. 21 we illustrate some of the symmetrical forms which are of 
common occurrence. Plate 11, Fig. 21 a and 5b illustrate what we call the 
spherically symmetrical forms, the simple disk and the simple ring respectively. 
Fig. 21 c, d, e, f, g, h and j are forms having biaxial symmetry about two axes at 
right angles. The first two we call the simple disk and the doubly indented disk. 
The next two are the simple ring and the doubly indented ring. The last three 
represent three examples of what we call complex biaxial forms. Plate 11, Fig. 21k 
is the helical form. Finally / and m represent centrally symmetric forms in which 
the structure is unchanged by reflection of each point of itself in the central star. 
We have proceeded on the system that no nebula which does not strictly fall into 
one of these classes should be classified as regular. The results of this classifica- 
tion for the 81 nebulae in our list and in Curtis’s list having dimensions over 
8 seconds of arc are as follows :— 


Spherical symmetry: I0 per cent, 
Biaxial symmetry : 45 per cent, 
Helical : 5 per cent, 
Central symmetry : I0 per cent, 
Irregular : 30 per cent. 


That is to say, 70 per cent of planetary nebulae exhibit a high degree of symmetry. 
Further, we find that classification on this basis is unambiguous. If desired, it 


* M.N., 106, 300, 1946. See also Eddington, Internal Constitution of the Stars, Cambridge, 
1930, p. 389, footnote. 

+ H. D. Curtis, Lick Obs. Publ., 13, 57, 1918. 

} R. H. Stoy, The Observatory, 56, 269, 1933. 

§ W. U1. Wright, Lick Obs. Publ., 13, 193, 1918. 
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is possible to split the first two classes into the sub-classes already enumerated, 
e.g. simple disk, simple ring and so forth, but then the classification by 
independent workers becomes more ambiguous. 

One final point of interest is provided by the question of projection. The 
maximum axis ratio exhibited by symmetrical nebulae is about 2:1. If this is 
the true form of planetaries of this type, i.e. if they are prolate spheroids with 
axes in the ratios 2:1:1, what proportion will appear circular as the result of 
projection? It can readily be shown that the proportion of nebulae with observed 
axis ratios in the range R to R+ dR is 


RdR 
V{3(4—R*)}" 
where R runs from ito 2. The integral of this, giving the proportion with axis 
ratios between I and I1+4, is 
I — 4/(1 —2x/3 — x?/3). 

This gives 3-5 per cent of cases with axes differing by less than 10 per cent and 
7:6 per cent of cases with axes differing by less than 20 per cent. That is to say 
that the observed percentage of spherically symmetric cases is much too high 
to be accounted for by projection alone. 

We hope that the material and comments presented here will be found to be 
in a form particularly adapted to the needs of theoretical workers. 


Radcliffe Observatory, 
Pretoria : 
1950 May 11. 
References 


(1) J. L. Dreyer, NGC and IC, Memoirs R.A.S., 49, 1888; 51, 1895; 59, 1908. 
(2) W. W. Campbell and J. H. Moore, “‘ Velocities of Bright Line Nebulae”’, Publ. Lick Obs., 
13, 77, 1918. 
(3) B. Vorontsov-Velyaminov and P. P. Parenago, Russian Astronomical F., 8, 206, 1931. 
(4) B. Vorontsov-Velyaminov, ‘‘General Catalogue of Planetary Nebulae”, Astronomical 
Journal of the Soviet Union, 11, 40, 1934. 
(5) D. H. Menzel, ‘‘ Five New Planetary Nebulae”, Harvard Bull., No. 777, 1922. 
(6) H. Shapley, ‘‘ Five Planetary Nebulae and a Globular Cluster’’, Harvard Bull., No. 902, 
26, 1936. 
(7) Jenka Mohr, ‘fA Planetary Nebula superposed on the Large Magellanic Cloud”, 
Harvard Bull., No. 907, 13, 1938. 
(8) R. Innes, Union Observatory Circular, No. 41, 345, 1918. 
(9) R. Innes, Union Observatory Circular, No. 61, 233, 1924. 
(10) W. S. Finsen, Union Observatory Circular, No. 69, 374, 1926. 
(1x) R. Innes, Union Observatory Circular, No. 73, 419, 1927; W. van den Bos, ibid., p. 424. 
(12) H. Shapley and J. S. Paraskevopoulos, ‘“‘ Photographs of Thirty Southern Nebulae and 
Clusters’’, Harvard Reprint, No. 184, 1940; and Proc. Nat. Acad. Sci., 26, 31, 1940. 
(13) H. Knox Shaw, ‘Observations of Nebulae, 1912-14”, Helwan Obs. Bull., No. 15, 129, 
1915. 
(14) C. C. L. Gregory, ‘Fourth List of Nebulae Photographed with Reynolds Reflector”’, 
Helwan Obs. Bull., No. 22, 219, 1921. 
(15) M. R. Madwar, ‘‘ Sixth List of Nebulae Photographed with the Reynolds Reflector”’, 
Helwan Obs. Bull., No. 38, 1935. : 
(16) H. Shapley and A. Ames, “‘ Survey of External Galaxies Brighter than 13th Magnitude”’, 
Harvard Annals, 88, No. 2, 1932. 
(17) J. H. Reynolds, “Spiral Nebulae in the Zone —40° to —go° (from the Franklin- 
Adams plates)”, M.N., 81, 598, 1921. 
(18) Delisle Stewart, Harvard Annals, 60, 174, 1908. 
(19) R. H. Frost, Harvard Annals, 60, 190, 1908. 











ON THE LUMINOSITY-VELOCITY RELATION OF CEPHEIDS 
Svein Rosseland 
(Received 1950 April 17) 


Summary 


The purpose of the paper is to clarify the ideas with respect to the phase 
relationship between luminosity and velocity of pulsating stars, which was 
recently discussed by Milne. Using a generalization of Milne’s procedure and 
known results of non-adiabatic pulsation theory, it is shown that for small 
deviations from adiabatic equilibrium and a constant ratio of specific heats, the 
luminosity is closely in phase with the negative displacement, irrespective of 
the stellar model employed. 





1. Introduction.—Since the inception of the pulsation theory of Cepheids 
various authors have attempted to explain quantitatively the observed approximate 
coincidence of maximum luminosity with maximum velocity of approach, which 
contrasts sharply with the common-sense expectation of maximum luminosity 
occurring at the time of highest compression. It was obvious that a small lag 
of the kind might be expected because of dissipation of wave energy by viscosity, 
conduction and similar causes. A Jarge lag like the observed amount of about 
70° was unexpected and inexplicable without postulating special conditions not 
envisaged in a standard star model, for example. 

In several papers Eddington (1) showed that the hydrogen ionization zone 
provides physical conditions of the right character to produce large phase 
differences, though it seems to be stretching the matter too far to say that this 
idea has been proved (or disproved) by quantitative means. One might also 
think of other means of transforming wave energy into thermal energy like 
prominence activity, turbulence in a wide sense, shock waves, etc. None of 
these possibilities has been investigated in a quantitative manner as yet. 

In a recent paper Milne (2) has reconsidered the phase-lag problem from a 
new angle. Stripped down to essentials Milne’s contribution to the problem is 
to show that, under certain assumptions with regard to the manner in which 
non-adiabatic conditions influence the oscillation frequency and the displacement, 
it is indeed possible to have the luminosity variation in phase with the surface 
velocity, a quarter wave later than the negative displacement. 

The mathematical technique used by Milne affords an interesting avenue of 
approach to the phase-lag problem. In the present paper the method has been 
worked out in greater generality, and it is shown that the physical assumptions 
underlying Milne’s results are not such as we usually associate with real stars. 
It is shown in a completely general way that the common-sense view of the matter 
is correct—that is: for small deviations from adiabatic conditions and with a 
constant ratio of specific heats, the luminosity will be very nearly in phase with 
the negative displacement. 

2. The luminosity.—The luminosity of the star L is given as the decrease in 
total energy E per unit time: 
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Let K, —Q, H and E, denote kinetic energy, potential gravitational energy, 
heat energy and nuclear energy of the star. Then 
L=-K+Q-H-E,. 
Using the equation of the virial 
}=2K-Q243(y—)H, 


where J is the star’s moment of inertia with respect to the centre and y the ratio 
of specific heats, the heat content H may be eliminated, giving 
I 


L= -B,+ 7 —5l(3y—4)2-(3y—5)K + 41], (x) 


3. Loss of nuclear energy.—Let « denote the loss of nuclear energy per unit 


time and mass. Then 


n 


where M is the total mass of the star and dm a mass element. Owing to the 
pulsations « will separate into a steady part, «9 say, and a variable part de. 
Assuming ¢ to depend on temperature T and density p by the power law 


«xp"T”, 


de pdp vdT 
— SS —- + “=> » 
€o p T 


this gives to the first order 


where dp and 87 are the first-order Lagrangian variations of p and T. In the 
far interior of the star conditions are very nearly adiabatic, and we may write to 
a sufficient approximation for the present purpose 


8T dp 
F one? 
which gives for 5e 
2 €o9 
Be = [n+ fy —1)] “S*. 


But 5p/p depends on the displacement ¢ by the equation of continuit 
P P y q y 


dp/p = —div&. 
Hence follows 
de = —[p+ ly—1)] diveg 


and finally for the rate of loss of nuclear energy 


‘ M M 
-E,= { codm— | [w+ Ay—1)] divGegdm 
0 0 


— E,, =Ly— [p+ y—1)] div Lo, 


M 
0 
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is the average luminosity of the star, and where a bar denotes an average taken 
over the stellar core where « differs sensibly from zero. 

This result shows that the variation in luminosity due to variations in the 
output of nuclear energy is essentially in phase with the displacement in the 
stellar core. Under nearly adiabatic conditions the displacement preserves a 
nearly constant phase throughout the star, making the oscillation very nearly a 
standing wave. In that case this term will be nearly in phase with the negative 
surface displacement and a quarter wave out of phase with the surface velocity. 
The only way of reversing this result is to introduce non-adiabatic processes in 
the part of the star outside the core and of such a nature as to convert the 
standing wave into a running one, as was indicated in the introduction. 

4. The dynamic terms.—Consider next the remaining terms in the right-hand 
side of (1). Introducing the proper expressions for K, Q and J we find for these 
terms 

: I -M m ; 2 
ea eS | : | @r-4) = — 437-5)? +3 ae) | dm. 


Here m denotes the mass inside a sphere of radius r concentric with the star, 
while G is the constant of gravitation. The radius r is here the distance of an 
individual mass element from the centre in the Lagrangian sense. Hence m 
and dm are constants during the motion. 


Let a denote the average value of r for an individual mass element and write 


r=ae". (3) 


Performing the differentiations involved, equation (2) assumes the form 


-M 
BL = = | "| (ay—4) Sat i+ Gy—DiP + y+ 2H | Adm.) 


Restricting the consideration to harmonic oscillations of infinitely small 
amplitude and of frequency o/27, it follows from (3) that » becomes equal to the 
relative displacement ¢/r, and (4) reduces to 


eo Gm 
Sis nce ery ee 
3(y— I) | ™ | (sy 4) r3 oc ] C m 


Under adiabatic conditions the right-hand side integral of (5) must, of course, 
vanish. ‘That this is so is easy to show directly from the adiabatic wave equation. 
Conversely, a non-vanishing 5L is conditioned by additional terms in o* and 
which are the consequences of non-adiabatic conditions. 

5. Non-adiabatic perturbations.—The question now arises as to how frequency 
and displacement are affected by non-adiabatic conditions. This problem has 
been discussed by various authors and the results are summarized in the book 
Pulsation Theory of Variable Stars, where, incidentally, the full formal solution of 
the problem is discussed (Chap. V, Section1). It was there shown that ordinary 
non-adiabatic perturbations do not change the oscillation frequency to the first 
order, but may produce a damping factor, and will always produce additional 
terms in the displacement. These terms produce a variable phase of the 
displacement, the amplitude remaining unaltered. The damping factor may be 
disregarded here, as it is believed to be a characteristic of Cepheid pulsations that 
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they are very nearly undamped. That leaves the phase as the essential perturbation 
element. Writing the displacement in the form 
{=£ cos (ot — 6) 


we are entitled to give both € and o their adiabatic values, while @ is the 
perturbation element, which may be assumed to vanish in the adiabatic limit. 
Since 6 is a small quantity we may write to the first order 


{= — o€ sin (ot — 6) = — of sinot + c£€0 cos ot = — of sinot + off. 


Introducing this expression for £ into (5) it follows that the integral over the 
term —o€sinot will vanish, since it corresponds to adiabatic conditions, while 
the integral over the second term gives 


o rae NR. 
sated on [ aya) SF ot | roca. (6) 
This gives a luminosity variation in phase with the displacement and not with the 
velocity. This result is seen to be independent of the particular stellar model 
under consideration. 

The course to follow in the further investigation of the phase-lag problem 
would now seem to be clear. It consists in a continued study of such conditions 
as will change the stationary adiabatic oscillation into a running wave. Perhaps 
a little more attention might be given to the possibility that the processes 
responsible for the energy generation in Cepheids may be of the delayed-action 
type. In that case all three elements of the wave—frequency, amplitude and 
phase, will be affected by the perturbation. However, it would seem more 
rational first to investigate quantitatively obvious sources of dissipation effects 
like the hydrogen ionization zone, shock waves in the atmospheric fringe, or 
similar phenomena. 


Blindern, Oslo : 
1950 April 11. 
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Summary 


Russell’s relationship between the planets and the white dwarfs is 
explained on the basis of the atomic theory of solids. A single mass-radius 
diagram can be constructed which refers to both the planets and the white 
dwarfs. The radius is an increasing function of the mass in the region of 
planetary masses, and it is a decreasing function of the mass for stellar 
masses. The radius is a maximum at some mass intermediate between 
planetary and stellar masses. The maximum radius depends on the 
chemical composition, and it is 85,000 km. for hydrogen. 





1. Introduction.—At the International Astronomical Union Meeting which 
was held in Paris in 1935, H. N. Russell (1) pointed out that the planets and the 
white dwarfs form a single family, assuming that they are chemically homogeneous 
and that the chemical composition is the same for all. This family has only one 
free variable, which may conveniently be chosen as the mass. All the other 
bulk properties, such as the radius and the mean density, are then determined 
by the mass. The mean density of a planet is low because its mass is small. 
The mean density of a white dwarf is very high because its mass is large. It 
should be possible to construct a single mass-radius diagram which refers to 
both the planets and the white dwarfs. 

Russell’s reasons for grouping the planets and the white dwarfs together in 
this way are both simple and convincing. He points out that both the planets 
and the white dwarfs are cold in the sense that the density at any internal point 
is determined by the pressure at that point. In other words, the influence of 
temperature is so small that it can be neglected to a good approximation. Thus, 
in the accepted theory of the white dwarfs it is assumed that the electrons 
constitute a degenerate Fermi gas at absolute zero temperature. If the density p 
at absolute zero temperature were known at every pressure p, then the bulk 
properties of the planets and the white dwarfs could be deduced on the assumption 
of hydrostatic equilibrium. The equation to be solved is 


dp/|dr = —pg, (I) 
where r is the distance from the centre, and where g is the local acceleration due 
to gravity. ‘The admissible solutions of this equation form a one-parameter 
family. A solution can be specified unambiguously by the central pressure, or, 
as is generally more convenient, by the total mass.* Thus the mass determines 
the radius and the mean density uniquely. Russell has gone one step further. 
Among the planets the radius increases with the mass. But the opposite trend 
prevails among the white dwarfs; the greater the mass of a white dwarf, the 
smaller is its radius. Russell therefore asserts that, somewhere between planetary 
and stellar masses, there must be a maximum radius for acold body. Observations 


* Under certain exceptional circumstances the mass does not determine the solution uniquely (2). 
But this is not important in the present connection. 
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on the planets and the white dwarfs are consistent with the view that this 
maximum radius is not much larger than the radius of Jupiter. 

This ingenious argument of Russell is irrefutable. It establishes beyond doubt 
a simple connection between the planets and the white dwarfs. But it does not 
go far enough. The relationship between the planets and the white dwarfs 
must be understood on the basis of atomic theory. On the one hand, the planets 
are solid (or liquid). The properties of a solid can be reduced to those of its 
constituent atoms or molecules. Solids have many individual characteristics 
because atoms and molecules differ widely. For this reason some solids are 
metallic, and others are non-meiallic. Also, solids can differ greatly in density ; 
for example, at low pressures gold is more than two hundred times as dense as 
solid hydrogen. Again, some solids are much more compressible than others ; 
at low pressures the compressibility of solid hydrogen is about a thousand times 
as great as that of gold. On the other hand, there is no such variety in the interior 
of a white dwarf star. This is because atoms do not exist at these pressures. The 
electrons form a degenerate Fermi gas, and an individual electron is not attached 
to any particular nucleus. All solids are metallic at these pressures. Also, 
the density at a given pressure is insensitive to the chemical composition; it 
depends only on the ratio of atomic weight to atomic number, which is about 
two for most elements. Again, the compressibility is a function of the pressure 
only, and it does not depend on the chemical composition at all. The properties 
of materials at the pressures in the interiors of the white dwarfs are comparatively 
well understood. The properties of solids at laboratory pressures are also, in 
principle, understcod. But little attention has been paid to the intermediate 
pressures, such as occur in the interiors of the planets. It is this range of 
intermediate pressures which is of greatest interest in the present paper. It will 
be shown that an understanding of the behaviour of materials at these pressures 
automatically explains Russell’s relationship between the planets and the white 
dwarfs. 

2. Transitions to metallic phases.—If the pressure on any non-metal is 
continually increased, a critical pressure will eventually be reached at which there 
is a drastic change of phase.and the material becomes metallic. ‘Transitions to 
metallic phases have been discussed in considerable detail in other papers (3, 4). 
The critical pressure for a phase transition of this type is usually of the order of 
magnitude of a million atmospheres, but in some cases it is much less. Both 
metallic and non-metallic allotropes of a few elements exist under laboratory 
conditions, notably tin, arsenic and phosphorus. In addition, detailed 
calculations have been made on metallic hydrogen (5, 6). Hydrogen would be 
a metal of the alkali type at pressures above 0-8 million atmospheres. The 
transition to the metallic phase is distinguished from other polymorphic transitions 
by the large jump in density. It is not uncommon for the density to double at a 
transition to a metallic phase. The density of solid hydrogen increases at the 
critical pressure from 0-35 to 0-77g./cm.*. The density of metallic arsenic 
is 5-1 g./cm.*, as compared with 2:0g./cm.* for yellow arsenic. In ordinary 
polymorphic transitions the density changes by only one or two per cent. In 
these transitions the atoms or molecules are merely rearranged into a different 
crystal pattern. But in a transition to a metallic phase the atoms or molecules 
themselves are partially destroyed, and electrons are set free to wander through 
the crystal. The author (3) has shown that the jump in density at the boundary 











446 W. H. Ramsey Vol. 110 


of the Earth’s core is due to the material becoming metallic. The pressure at 
the boundary of the core is 1-4 million atmospheres, and the density increases 
by more than 50 per cent. 

Consider a solid composed of ions with closed shells (for example, an ionic 
crystal) or of saturated molecules (for example, solid hydrogen). The electrons 
in such a solid cannot conduct electricity; a net current, even in the presence of 
an electric field, would violate the Pauli exclusion principle.* At high pressures 
these solids are difficult to compress because of the strong repulsions between the 
closed shells or between the saturated molecules. The solid could collapse into a 
dense metallic state if these units were broken, but to break thern would require a 
large energy, often in the neighbourhood of ten electron volts per molecule. 
Nevertheless, if the pressure is increased, a critical pressure P will be reached at 
which the work done by the hydrostatic forces during the collapse to the metallic 
phase supplies just the energy required. This is expressed by the usual condition 
for the equilibrium of two phases at absolute zero temperature: 

U,—U,=P(V,—V,), (2) 
where Uand Vare the internal energy per mole and the volume per mole respectively, 
and where the subscripts 1 and 2 refer to the initial and final phases respectively. 





log p 


/ 
/ Fermi Gas 











log p 


Fic. 1.—A schematic pressure-density relationship. 


The metallic phase is stable at pressures higher than the critical pressure P, but not 
at pressures indefinitely higher. In general, the metallic phase also contains ions - 
with closed shells. The transition to the metallic phase released the electrons in 
the outermost shell, but the inner shells remain intact. ‘The inner shells will be 
broken at still higher pressures by phase transitions of the type described above. 
Thus, if the pressure on any material is increased indefinitely, it will undergo a 
series of drastic phase transitions. Each phase transition will increase the number 
of “ free ’’ electrons by breaking the saturated bonds of molecules or the closed 
shells of ions. This series of transitions will continue until every atom and mole- 
cule has been completely broken down. The electrons will then form a degenerate 
Fermi gas, which is imperfect only because of the electrostatic interactions. 

A schematic pressure-density relationship is shown in Fig. 1. The density 
increases at each phase transition, and also in each phase. It has been shown in 


* A survey of the atomic theory of solids has been given by Seitz (7). 
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another paper (4) that, for any particular phase, the bulk modulus k is approximately 
a linear * function of the pressure p: 


k=ky +ap. (3) 

The bulk modulus is the reciprocal of the compressibility. The constant kp is 

usually in the neighbourhood of a million atmospheres, and it is of little importance 

at the very high pressures which are of interest in this paper; it will therefore be 

neglected. On integration, equation (3) gives the following pressure-density 
relationship : 

log p=a~! log p + constant. (4) 


Strictly speaking, the factor « is not constant but decreases slowly with pressure. 
But the total decrease is usually small in the range of pressure in which any one 
phase is stable. The segments referring to each phase in Fig. 1 are therefore 
approximately straight lines. ‘The constant « is closely related to the size of the 
ions. Itis large for large ions, and it is small for smallions. For solids composed 
of only very heavy atoms, such as caesium iodide, « has the value 5. The smallest 
possible value of « is 3, and this is attained only in a perfect Fermi gas. Every 
material at high pressures has a value of « which lies between § and 5. This 
correlation between « and the size of the ions means that materials with large atomic 
weight are less compressible than those with small atomic weight. Also, this 
difference in compressibility persists until all the atoms have been completely 
broken down. These conclusions are the result of a detailed investigation (4), 
but the necessity for them can be easily demonstrated. At low pressures solids 
with large atomic weight are up to 200 times as dense as those with small atomic 
weight. But at the pressures in the interiors of the white dwarf stars the density 
is almost independent of the atomic weight. The difference in density can be 
wiped out only if the denser materials are less compressible. The difference in 
compressibility must be maintained to very high pressures indeed. At four 
million atmospheres, which is the pressure at the centre of the Earth, gold is still 
more than thirty times as dense as metallic hydrogen (cf. Table I). Also, gold 
will undergo phase transitions at higher pressures, but hydrogen will not. The 
compressibility is independent of the atomic weight only at pressures more than 
a million times as great as that at the centre of the Earth. This result is important, 
as Bullen (8) has recently suggested that the compressibility is insensitive to the 
chemical composition at the pressures prevailing in the Earth. Bullen advanced 
this hypothesis because of the similarity of the compressibility-pressure relation- 
ships in the Earth’s core and mantle, which were generally believed to be chemically 
distinct. It has been shown in another paper (4) that this similarity is due to the 
fact that the materials of the core and mantle are chemically identical. The general 
ideas behind Bullen’s compressibility-pressure hypothesis can be reformulated, 
and the hypothesis may become a very useful guide in interpreting seismic data. 
A transition to a metallic phase automatically decreases the size of the ions, as it 
releases electrons by breaking up closed shells. For this reason the constant « 
always decreases at a transition to a metallic phase, but the decrease is usually 
small (4). According to equation (4), this means that the slope in Fig. 1 always 
increases slightly at a phase transition. ‘The material is still not a perfect Fermi 


* Deviations from the linear relation (3) occur at pressures below about 0°5 Rp. But this is of no 
importance in the present paper. 
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gas when the last closed shell has been broken. ‘The electrons interact electro- 
statically with the nuclei and with one another. Asa result, the density is always 
greater than that of a Fermi gas at the same pressure, and the compressibility is 
always less. ‘The magnitudes of all the electrostatic interactions are inversely 
proportional to the lattice spacing. But the Fermi energy, which arises from the 
Pauli exclusion principle, is inversely proportional to the square of the lattice 
spacing. Consequently, as the pressure increases, the Fermi energy becomes 
more and more preponderant. But, at high pressures, the Fermi energy is just 
the energy of a perfect Fermi gas. Hence the pressure-density relationship 
asymptotically approaches that of a Fermi gas, as shownin Fig.1. The equation 
of state of a Fermi gas at absolute zero temperature is * 


p=KZ*V-5, (5) 


where V is the molar volume and where Z is the number of electrons per molecule. 
The constant K is independent of the chemical composition, and its value is 


K=10-04 x 10" dynes cm.°, (6) 


If the pressure p is expressed in atmospheres and the volumeV incm.°, the numerical 
value of K is 9-9 x 10°. The mass of one mole is A grams where A is the molecular 
weight. The pressure-density relationship is therefore 
p=K(Z/A)P*p%. (7) 

The only influence of the chemical composition on this equation is through the 
ratio(A/Z). This ratio is about two for most materials, but it is unity for hydrogen. 
For an element, Z is the atomic number and A is the atomic weight.t 

Hydrogen becomes metallic at a pressure of 0-8 million atmospheres, and above 
this pressure atoms and molecules do not exist. But very much higher pressures 
are required to break down the atoms of heavy elements completely. An order of 
magnitude estimate of the pressure required can be made as follows. Consider 
the last phase transition in an element of atomic number Z. The critical pressure P 
for the transition satisfies equation (2). Ions will not exist at pressures above this 
critical pressure P. As the density always increases by a large percentage at a 
transition to a metallic phase, the change in volume (V,-V) will be comparable 
with the final volume V,. Equation (2) may therefore be simplified to 


U,— U,=PV;3. (8) 


Apart from the electrostatic interactions, the final phase is a Fermi gas. The 
equation of state (5) for a Fermi gas may therefore be used to express the volume V, 
in terms of the critical pressure P. Equation (8) then gives the critical pressure P 
in terms of the energy difference (U,— Uj): 


U,— U,=K**ZP%, (9) 


The last phase transition corresponds to the breaking of the K shell, which is 
the most tightly bound shell in the atom. The energy required to break this 
shell by excitation is approximately Z* Rydberg units{(z0). The energy 


* Chandrasekhar (9) has given a detailed account of the Fermi gas approximation and its 
application to the white dwarfs. 

+ At the highest pressures in the white dwarfs equations (5) and (7) should be replaced by the 
corresponding relativistic expressions ; cf. Chandrasekhar (9). 

t One Rydberg unit is the energy necessary to ionize a hydrogen atom; it is about 
13°5 electron volts. 
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difference (U, — U,) is therefore about Z? Rydberg units per atom. One Rydberg 
unit per atom corresponds to 1-3 x 10! ergs per gram-atom; hence 


U, — U, =1-3Z? x 10" ergs. (10) 
Equations (6), (9) and (10) lead to the following expression for the critical 


pressure P: 
P=Z®"2 x 10 dynes/cm.?. (11) 


One atmosphere is approximately 10° dynes/cm.?. The atomic number Z is 
nearly a hundred for the heaviest elements. A pressure of 1012 atmospheres is 
therefore required to break down the heaviest atoms completely. Such pressures 
do not exist in the planets, but pressures even thousands of times greater occur 
in the interiors of the white dwarf stars. The pressure at the centre of Jupiter, 
which is the highest in the planets, is 3 x 107 atmospheres. This pressure could 
break down only the lightest atoms completely, possibly only hydrogen and 
helium. 
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Fic. 2.—A schematic mass-radius diagram. 


When the density is known as a function of the pressure, the bulk properties 
of the planets and the white dwarfs can be computed on the basis of the hydrostatic 
equation (1). The resultant mass-radius diagram is shown schematically in 
Fig. 2, and it has all the characteristics predicted by Russell (x). The internal 
pressures are small when the mass is small, and so atoms exist as units. For small 
masses, such as planetary masses, the radius increases with the mass. On the 
other hand, the internal pressures are enormous for masses of stellar magnitude, 
and atoms exist only in a thin surface layer. The equation of state (7) for a Fermi 
gas may be used for these masses, except in the surface layer in which atoms 
exist. The thickness of this surface layer decreases rapidly with the mass, and 
so the asymptotic behaviour of the mass-radius relationship is given by the usual 
theory of the white dwarfs. The mass-radius relationship corresponding to the 
Fermi gas approximation* is indicated by the broken curve in Fig. 2. As Russell 
pointed out, there is a maximum radius for a cold body. The value of the 
maximum radius depends on the chemical composition. 


* Cf. S, Chandrasekhar (9). ‘ 
22 
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This general treatment will now be illustrated by precise calculations on 
hydrogen. These calculations are more refined than is necessary for the present 
purposes, but they will be of value in a future paper on the internal structures and 
constitutions of the major planets. It turns out that Jupiter and Saturn are 
predominantly hydrogen ; the proportion of hydrogen in these planets is about 
the same as that in the Sun’s atmosphere. 

3. Metallic hydrogen.—From the point of view of atomic structure, metallic 
hydrogen is the simplest possible solid. It is simpler than the other alkali metals 
in that the ion core is a proton, which rigorously has a Coulomb field. 
Consequently many of the properties of metallic hydrogen can be predicted 
accurately by not too laborious calculations in wave mechanics. Calculations 
on metallic hydrogen were performed first by Wigner and Huntington (5), and 
more recently by Kronig, de Boer and Korringa(6). The numerical predictions 
of the two sets of calculations are not essentially different. The calculations of 
Wigner and Huntington are very accurate in the region of low pressures. But the 
author has found the approximations developed by Kronig, de Boer and Korringa 
more suited to the high pressures prevailing in the interiors of the major planets. 
The present calculations do not contain any new mathematical methods or 
approximations. The calculations therefore need not be described in detail. 
A numerical pressure-density relationship is derived by the methods developed 
by the previous workers (5, 6). 

The calculations on metallic hydrogen will now be briefly outlined in order 
to clarify the connection between metallic hydrogen and a Fermi gas. To facilitate 
a comparison, the notations of the original papers (5,6) will be used. The 
energy per atom of metallic hydrogen will be written as E Rydberg units, where E 
is dimensionless. The volume per atom will be written in the form 


v=(47/3)r;a), (12) 
where 4a, is the radius of the first Bohr orbit of the hydrogen atom, and where 7, is 
a dimensionless parameter which is proportional to the lattice spacing. The 
energy per atom of a metal is usually written as the sum of the three terms (7) : 


E=E,+ F+W. (13) 


The first term E, is the lowest energy which an electron can have in the periodic 
field of the lattice; it corresponds to an electron moving through the lattice 
with zero velocity. The second term F is the Fermi energy, which corresponds 
to the kinetic energy which the electrons must possess by virtue of the Pauli 
exclusion principle. Finally, W is an electrostatic term which arises partly from 
the exchange of electrons and partly from the correlation of their positions. The 
Fermi energy F is positive, and it prevents the lattice from collapsing. The last 
term W is also positive, but it is at least an order of magnitude smaller than the 
Fermi energy. The first term E, is negative, and it provides the binding of the 
lattice. Approximate expressions for Ey, F and W are* 


Ey= —3/Te (14) 
F =2-21/r, (15) 
W =(0-284/r,) —0°58/(r,+ 5-1). (16) 





* Equations (14), (15) and (16) above should be compared with equations (13), (29) and (37) 
of reference (6). 
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These expressions are accurate only when r, is small, that is, when the pressure is 
high. Anumber of corrections must be applied at low pressures ; the corrections are 
given in convenient form in the original papers on metallic hydrogen (5,6). The 
importance of the corrections decreases rapidly with pressure. They are already 
small at a million atmospheres, and they are of little importance above ten million 
atmospheres. The calculations are therefore most reliable at high pressures, 
The region of low pressures is of little interest, as metallic hydrogen is unstable. 
It will be noted that, at high pressures, the electrostatic energies E, and W are 
proportional to r,-, but the Fermi energy F is proportional to r,-*. The Fermi 
energy therefore becomes more and more preponderant as the pressure is increased, 
But the Fermi energy (15) is identical with the energy of a perfect Fermi gas. 
Hence, as the pressure is increased indefinitely, the energy of metallic hydrogen 
tends asymptotically to that of a perfect Fermi gas. In the limit of very high 
pressures, the equation of state of metallic hydrogen is therefore identical with that 
of a Fermi gas. 

The best available calculations on solid molecular hydrogen are those of Kronig, 
de Boer and Korringa (6). De Boer (1x1) has shown that the interaction energy 
between two hydrogen molecules can be written as follows: 


$(x) =(A/x"*) —(B/x*), (17) 
where x is the distance between the two molecules. If the energy ¢(x) is expressed 


in Rydberg units and the distance x in units of the Bohr radius ap, the numerical 
values of the constants A and B are: 

A=7:76x10°8, B=27-0. (18) 
The internal energy of the solid is obtained by summing the interaction energy (17) 
over the different pairs of molecules in the crystal. The present calculations 
differ from those of Kronig, de Boer and Korringa only by a refinement ; the zero- 
point energy of the molecules has been taken into account by an approximate 
method due to Debye, cf. reference (7). The zero-point energy exerts a strong 
influence at low pressures. If it is neglected, the interaction (17) leads to a density 
of 0-144 g./cm.® at zero pressure, which is more than 60 per cent greater than the 
empirical density of 0-089 g./cm.’. The correction to the density is much less at 
higher pressures, and over most of the pressure range which is of interest the 
correction is about 7 per cent. The zero-point energy of the protons in metallic 
hydrogen has also been taken into account. 

The critical pressure at which molecular and metallic hydrogen can co-exist in 
equilibrium is given by equation (2). According to the present calculations the 
critical pressure* is 0-8 x 1012 dynes/cm.?. At the critical pressure the densities of 
the molecular and metallic phases are 0-35 and 0-77 g./cm.® respectively. These 
results do not differ appreciably from those of Kronig, de Boer and Korringa.t 
These workers estimated the critical pressure to be 0-7 x 10! dynes/cm.?, and 
at the critical pressure the densities of the two phases were found to be 0-4 and 
o-8 g./cm.*. The results of the present calculations are given in Table I, where 
the density is shown as a function of the pressure. Below the critical pressure 
the densities refer to the molecular phase, and at higher pressures to the 
metallic phase. As regards the reliability of these data, the uncertainties in the 


* One atmosphere is 1-01 X 108 dynes/cm.?. 


+ The small differences between the present results and those of Kronig et al. arise, of course, 
from the inclusion of zero-point energies. 
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densities of the metallic phase are too small to be of any consequence. The 
densities of the molecular phase are less reliable, and an error of five per cent would 
not be surprising. The critical pressure is the most uncertain quantity in Table I, 
and it may be in error by ten, or even twenty, per cent. 


Tas.e I 
Pressure-Density Relationship for Hydrogen 
Molecular Phase Metallic Phase 
Pressure Density Pressure Density 
10!* dynes/cm.* — g./cm.? 101? dynes/cm.2_—_g./cm.? 
° 0089 o'8 o'77 
orl 0°225 1° 0°81 
o'2 0°261 2°0 0-98 
0°3 0°284 3°0 1°12 
0-4 0°302 4°0 1°25 
o's 0°316 50 1°36 
06 0°329 100 1°84 
o'7 0°340 20°0 2°53 
08 0°*350 30°0 3°07 

4. Hydrogen planets.—The bulk properties of planets composed of hydrogen 
have been computed on the basis of the pressure-density relationship given in 
Table I. It has been assumed that the interiors of the planets are in hydrostatic 
equilibrium, so that equation (1) is applicable. ‘This equation has been solved 
for a number of arbitrarily chosen central pressures, using step-by-step numerical 
integration. The solution of equation (1) determines uniquely the mass, radius and 
mean density of the planet. A selection of the results of the calculations is given 
in Table II. The unit of mass in this table is the mass M, of the Earth. The 
central pressure in a small planet is too low to support a core of metallic hydrogen. 
A core is possible only if the mass of the planet exceeds 88 My. A planet whose 
mass lies between 88 M, and 95 M, can exist either with or withoutacore. This is 
a consequence of the fact that the density of hydrogen jumps by more than 50 per 
cent at the critical pressure ; a full explanation has been given in another paper (2). 
A planet always has a core if its mass exceeds 95 Mg, and the fraction of the total 
mass which is contained in the core increases very rapidly with the mass of the 
planet. Fora planet of mass 97 My, the core contains 33 per cent of the total mass ; 
but for a planet of mass 1800 Mg, the core contains over 99 per cent of the mass. 
The radius of a planet of mass 1800 M, is 83,000 km., but the thickness of the 
mantle of molecular hydrogen is less than 3000km. The properties of molecular 
hydrogen have therefore little influence on the mean density of a planet of this mass, 
and even less influence for more massive planets. The radius of the planet is a 
maximum when the mass is about 1000 My. The maximum radius is 85,000 km., 
and this is the largest possible radius for a cold body composed of hydrogen. As 
Russell (1) suspected, the maximum radius is not much larger than the radius of 
Jupiter, which is69,900km. For masses exceeding 1000 Mg, the radius decreases 
with the mass and the mass-radius relationship asymptotically approaches that 
given by the usual theory of the white dwarfs. 

5. Kothari’s theory of the planets.—As long ago as 1936, Kothari (12) pointed 
out that atoms cannot be treated as structureless particles at the pressures prevailing 
in the planets. An atom consists of a nucleus surrounded by a complex system of 
electrons, and this internal structure must be taken into account. Kothari 
developed a simple theory which can be regarded as a second approximation to the 
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Fermi gas treatment. Allowance is made for the tendency of the electrons to 
cluster around the nuclei. This clustering is of little importance at the pressures in 
the interiors of the white dwarfs, but it is more pronounced at lower pressures. 
Kothari’s treatment leads to the formation of atoms when the lattice spacing is 
large; the Fermi gas treatment, of course, does not. Kothari has shown that his 
theory leads to a mass-radius diagram which is similar to that illustrated in Fig. 2. 
He therefore claims that his theory explains Russell’s relationship between the 
planets and the white dwarfs. In Kothari’s theory the maximum radius for a cold 
body depends on the chemical composition, and it is about 125,000km, for 
hydrogen. 
TasLe II 
Hydrogen Planets 
Mass Radius Mean Density 
(Mp=1) (km.) (g./cm.*) 
10°7 46,000 0°16 
60 73,000 0°22 
97 75,600 0°32 
208 77,100 0°65 
317 795400 9°99 
623 82,600 1°57 
1000 84,500 2°36 
1400 83,700 3°42 
1800 83,000 4°53 


Kothari’s theory has not been generally accepted as its predictions are 
incompatible with the observed properties of materials, not only at the pressures 
attained in the laboratory but also at the pressures in the interior of the Earth. 
The principal weakness in the theory is that the interactions between neighbouring 
atoms are ignored completely. For definiteness, let us take hydrogen as an example. 
The internal energy of metallic hydrogen is given by equation (13). The 
electrostatic energies E, and W are ignored in the Fermi gas approximation. 
Kothari’s theory is an improvement on the Fermi gas treatment in that an attempt 
is made to include the energy Ey, which corresponds to the interactions between 
the electrons and the nuclei. No attempt is made to include the term W, which 
arises from the mutual interactions of the electrons and is the metallic equivalent 
of interatomic forces. If the lattice spacing of metallic hydrogen is increased 
indefinitely, the result will be an assembly of unexcited hydrogen atoms. Kothari’s 
theory also leads to hydrogen atoms when the lattice spacing is large. But Kothari 
does not find solid molecular hydrogen, which is the stable modification at pressures 
below 0-8 million atmospheres. The failure of the theory is due to the fact that 
the interactions between atoms are ignored. If Kothari had taken the interactions 
between the atoms into account, he would have found the transitions to the metallic 
phases which are described in this paper. A transition toa metallic phase, like every 
other phase transition, is inherently a property of groups of atoms, and it cannot 
be understood if the interactions between the atoms are ignored. 

In his calculations Kothari paid special attention to iron, as it is one of the most 
abundant elements in the Earth. His theory predicts that the density of iron will 
be 70 g./cm.’ when there is one “free” electron per atom. This is nine times 
the actual density of iron. A density in the neighbourhood of 70 g./cm.* is not 
attained even at the pressures in the Earth’s core, which Kothari assumed to be 
iron. Kothari therefore replaced the density of 70 g./cm.® which appears in his 
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theory by the empirical density of iron. He then compared the “ theoretical ”’ 
mass-radius relationship with the data for the terrestrial planets.* Such a 
comparison is valueless when the theory was in error by an order of magnitude at 
planetary pressures. Kothari’s theoretical densities refer to the metallic phase in 
which the atoms are completely broken down. According to equation (11), this 
phase is stable only at pressures more than ten thousand times as great as the 
pressure at the centre of the Earth. Kothari’s calculations have therefore no 
bearing on the constitutions of the terrestrial planets. The density of iron at low 
pressures is comparatively small because of the strong repulsions between the 
closed shells of the ions. ‘The compressibility of the material of the Earth, which 
can be derived from seismic data, shows that the repulsions between the ions 
predominate at all terrestrial pressures (4). But these repulsions are ignored 
completely in Kothari’s theory, and so it is inapplicable. Kothari’s work has been 
followed by a number of investigations (13, 14, 1§) along the same lines, but, 
so far as the planets are concerned, no essential improvement has been effected. 

In general, Kothari’s theory is inapplicable at planetary pressures. But his 
numerical calculations on hydrogen have a bearing on the structures of the major 
planets. This is because a comparatively small pressure breaks down hydrogen 
molecules completely. It has been shown that the properties of molecular 
hydrogen are of secondary importance for planets whose mass is a few thousand 
times the mass ofthe Earth. Also, Kothari’s treatment gives the density of metallic 
hydrogen correctly to about 25 per cent at planetary pressures. Kothari’s estimate 
of the maximum radius for a cold body composed of hydrogen is 125,000 km., and 
this is in error by only 50 per cent. 

Part of this investigation was carried out at the Institute of Theoretical Physics 
in Copenhagen, and the author would like to thank Professor N. Bohr for his 
hospitality. He would also like to express his appreciation to Professor L. 


Rosenfeld and to Professor Harold Jeffreys for the interest which they have shown 
in this work. 
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DOUBLE STAR MEASURES—3rd SERIES 


W. P. Hirst 
(Received 1950 June 5) 


The following measures were made, unless noted otherwise, with the 7-inch 
refractor at the Royal Observatory at the Cape by kind permission of H.M. 
Astronomer, Dr J. Jackson. A few were made with the 26$-inch refractor at 
the Union Observatory, Johannesburg, by courtesy of the Union Astronomer, 
Dr W. H. van den Bos. These have been noted accordingly. 

The columns from left to right give: Name and ADS number (if any); 
right ascension and southern declination (equinox 1900.0); date (omitting 
century); position angle; distance; magnitude difference (estimated); number 
of nights; remarks. 





Date 
1900 


Star Position AM Nights Remarks 





45°616 
48°134 cay p (est) 0”°5 
45°964 
45°964 . . , p (est) 1°-2 
45°950 
I 27 45°930 
h 3423 45°913 
h 3494 45°964 . “ae p (est) 1°-0 
Jc 8 AB 45°045 3 
46-088 p (est) 0°*5 


h 3556 AC 46°088 


A 16 45°036 
f Erid 46-088 





B 1004 45°044 


Gale 1 45°048 





Vol. 110 





Position 


4AM Nights Remarks 





B 744 
ADS 3159 


Riimker 4 


B 184 


h 3683 


h 3752 


I 276 


4 23 
Sirius 
h 3997 
B 101 
Hu. 115 
= 1216 0°45 
I 489 0°45 est 
4 205 


0°55 est 


B 208 1°97 


1°95 


% Arg 


I 292 


113 


p unreliable 
p (est) o”"9 


U.O. 263 inch 


Rough measure 


U.O. 26} inch 


U.O. 264 inch 


U.O. 263 inch 


0"4 


One measure 
U.O. 264 inch 
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Double star measures—3rd series 





Star Position 


Date 
1900 


AM 





I 294 


I 78 


h 4478 


a Crux 


I 298 


B 343 


Howe 28 


R 227 


B 1197 


Slr 19 


h 4687 


a Cent 


I 236 


B 239 


46°386 


47°435 


47°435 


46°386 


46°447 
47°435 
49°541 
46-386 
47°435 
49°538 
46°548 
46548 
47°435 
49°537 


46°548 
46°546 
46°534 
49°537 
46°548 
46°548 
46°548 
46°548 
46°434 
49°533 


46°548 


46°546 
49°543 


° 
75°6 


p (est) ”-o 





W. P. Hirst 





Position 0 p AM Nights Remarks 





h m ° ” 
h 4728 14 58°3 75°6 1°68 
46° 40° 74°6 1°76 


h4753 AB 15 11°6 142°9 
47° 30° 142°8 


h 4753 AC 129°2 
h 4753 BC » 129°8 


h 4757 15 15°4 48°1 
58° 58’ 47°3 


A258 AB 15 55°4 284°4 
57° 30° 277°0 
268-2 


h 4825 AC 245°3 


A 268 176°5 


Cor 197 147°2 


HII 19 _ 
349°9 


Sh 243 170°0 
167°9 


Brs 13 222°7 


Melb 4 AB 286°3 
283°8 
280°1 
278°2 
274°3 
256°6 
257°6 
257°3 
112°5 
III‘3 

4°2 
4°9 
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Double star measures—3rd series 





Star Position 


Date 


AM 
1900 


Nights Remarks 





h 
Riimker 26 20 


62° 


h 5246 21 


54° 


h 5258 21 


53° 


B 766 21 


41° 


= 2909 22 


°° 


Jc 20 23 


44° 


Cape Town: 
1950 April 4. 


m 
43°3 
48’ 
03°! 
59° 


12°7 
52° 


18-0 
26’ 
23°7 
31’ 
O1'2 
04’ 


45°868 


45°868 


45°835 


43°931 


45°818 


43°931 


p unreliable 


p doubtful 


p doubtful 





' 
' 
i 
eH 
i 
' 
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DYNAMIC EFFECTS OF A LIQUID CORE (II) 


Harold Feffreys 
(Received 1950 May 20) 


Summary 


An approximate method is applied to the Eulerian and forced nutations to 
take account of possible effects of differences of density within the Earth’s 
core. It is found that all such effects are adequately taken into account by 
replacing the actual core by a homogeneous one with the same moments of 
inertia. 





1. Heterogeneity of the core.—It was shown in an earlier paper* that the 
inertia of a liquid core in the Earth, with a rigid shell, would lead to too short 
a period for the Eulerian nutation and to a reduction of the theoretical amplitude 
of the 19-yearly nutation to substantially less than the observed value. Both 
effects were afterwards shown to be mitigated by allowance for a finite rigidity 
of the shell.t In both papers the core and shell were taken as homogeneous. 
The variation of mechanical properties within the shell could be taken into account 
by a numerical method, without, I think, introducing any new questions of 
principle. Variation of density in the core might introduce complications of a 
new type. The displacement of the core boundary is a rotation. Within the 
core there is an additional motion involving distortion. This gives rise to no 
gravitational energy when the density is taken as uniform, since it does not alter 
the density at any point. But such a motion, with the density increasing toward 
the centre, would be a distortion of a fluid mass from a stable steady state, and 
would be opposed by gravity. Hence it appears that with a heterogeneous core 
gravity might reduce the motion of the core relative to the shell, and that the 
results might again approach those for a completely rigid Earth. 

As the time of travel of a longitudinal wave through the core is very short 
compared with a day, it is safe to treat the core as incompressible. At present 
we return to a rigid shell. The motion in the core will consist of three parts: 
the rotation of the shell; a linear distortion producing no normal motion at the 
boundary; and a superposed distortion varying with distance from the centre 
in such a way that the normal displacement due to it also vanishes at the boundary. 
This allows the displacement to contain cubic terms in the coordinates. With 
certain obvious changes from the notation of Paper II, the most general cubic 
forms satisfying these conditions are, to the first order, 


a a 2 
f=x+ (14 2h)2+22(48-1 =): 

a a 2 
H=y+ (m +$m)2+$2(mS—m, 5), 

c c c x? + y? 


c c c x? + y? 
—y{ m+ £m, +m, © S — ms(x - a@ )}. 








* H. Jeffreys, M.N., 108, 206-209, 1948, referred to as Paper I. 
t H. Jeffreys, M.N., 109, 670-687, 1949, referred to as Paper II. 
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24 v2 
where S == ~ + a (2) 
These do not include displacements about the axis of z because they are of the 
second order and would lengthen the writing greatly. The /,m displacement 
represents the general rotation of the shell, /,, m, the linear distortion. /,, mg give 
displacements at all points along similar ellipsoids. To represent the kinematic 
possibility that the motion might be along level surfaces we need an additional 
displacement capable of giving motion with a component normal to similar 
ellipsoids, and this is provided by /,,m,. ‘Then the system (1) is sufficiently 
general to represent the main features of the actual motion. It does not satisfy 
the equations of motion exactly, and a full solution would be very difficult. But 
it is probably near enough. For small variations from the true motion, 
Hamilton’s integral taken over a period varies by a quantity of the second order 
in the variations; and Rayleigh’s method of deriving the principle of stationary 
periods from Hamilton’s principle applies equally well to gyroscopic systems. 
Hence the error in the computed free periods will be very small. For a forced 
motion some of the displacements found will have errors of the first order, but 
as the quantities observed are only indirectly related to those that are in error 
it may be again expected that the form (1) is accurate enough. But then we 
can use Lagrange’s equations as before. 

We can indeed go a little further. We are interested in the effects of variation 
of density in the core and of the differences of the ellipticities of surfaces of equal 
density. The former are of the order of 20 per cent on the model that we shall 
use, the latter of the order of 8 per cent. Without these variations the cubic 
terms would not occur. It seems likely, then, that the omission of the terms in 
1,, mg, /,,m3 would itself be a first-order error in Hamilton’s sense, and that any 
error due to their omission would be of the order of a few per cent. If a large 
correction is found it may be worth while to make a closer approximation. 

In Paper II the displacements were evaluated to the second order, on account 
of the presence of terms w&y(— €y), 4w*X ("+ y”) in the kinetic energy. If 
we restrict the adopted form of the displacements to the parts containing 
l,m, 1,,m,, the second-order terms will be the same as before. For the other 
parts, if it should be necessary to retain them, we should have to add second-order 
parts to &, n, { so as to make 0(£, n, £)/0(x, y, 2) unity to the second order and leave 
a boundary particle on the boundary, also to the second order. This, however, 
is unnecessary. In Xyéy the extra term is Dyx*(0"y/dxdt), which is a derivative 
with regard to the time and does not affect the equations of motion. Zé? will 
contain terms in /?, m?, which will be the same as for a rigid body and can be written 
down at once; but also terms linear in /, m and in /,,/,..., and terms in squares 
and products of /,,/,.... Since these terms contain no derivatives with regard 
to the time, it is indifferent whether we place them in the kinetic energy or the work 
function. Now a constant angular displacement in space of the whole body, 
originally in steady rotation, leaves it in steady rotation. Hence, when the 
equations of motion are formed for /,,/,..., the coefficient of / must vanish for 
disturbances of speed —w. This effectively determines the coefficients of ll, and 
analogous terms. Again, the steady rotation is a stable state and the surfaces 
of equal density are the surfaces of constant geopotential U + 4w*(x?+y*). The 
part of T that depends only on the displacements can be added to the work function; 
but if this is done the only contribution to the work function from 4, 4, ...comes 
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from displacements normal to the level surfaces, and can be found to the second 
order as for departures from a spherical gravitating mass, since the ellipticities 
are small. 

If » is an element of mass, we introduce the following constants: 


Ay, =Z px? =4C); A33 =X ps" =}(24,—C)), 
Ag, = A,;(I — 2¢€,), c=a(I-e), C,-—A,=«,C\, e—e€,=€', (3) 
Fy =25 Ay, =A,(I-«’), Fy=2 © Ay =Ay(1 +e’). 


For a homogeneous core F,=F;; this is no longer true. We neglect e? and 
similar quantities. ‘Then, if we retain only /, m,1,,m,, 


' ; «;*¥° Os « 
27 ge = S34 (i + ; i) + (i + - i \ 
S x2 € ;\? , € . \? 
+= px {(i+ £4) + (i+ San) | 
, a me Fe : 8 
+2wEpet{(1+ 1) (i ae it) = (i+ = i,)(m ca - m,)} 
—w(C, — A,)(2 + m*) + 20? (s ps? ; —X px? <)t +mm,) 


+term in /?+ mi. (4) 
Adding the kinetic energy of the shell we have 
2T = A(I2 + m?) — w(C — A)(2 + m?) + w(2A — C)(Im— Im) 
+ 2A, (Il, + mm,) + wF (Lm + ln, — lm, — 1m) 
+Cya(1 +2¢’)(1,n, —1,m,) + A, (2 +m?) 
+ 2w*e’ A,(Il, + mm.) + term in [? + m?. (5) 
2. Forms of «' and of the normal displacement in the core.—We take 1, for the 


mean radius of a level surface in the core, and its ellipticity is temporarily denoted 
by e. Take 


9% ‘ s 
=pPy— Pi e=@,+e, = 
P=Po- P1 a2’ ot &1 Fe» 


so that e of (3) is eg+e,. The mean moment of inertia is 


87 a 
1 == (b0—$nr)a’ 


ra r2\ d r2 
Ji (rom os) ae{r (coer is) 


Bie ied 5 7 5 
= — (oP + €1Po— F€0P1 — §1P1)@”s 


T5 
C,-A, _ 4 Pr (9) 
I 63 Po— 71 
To the first order in the ellipticity, if a’, c’, e’ are the semi-axes and ellipticity of a 
level surface in the core, the normal displacement is 
2(€o + ey m e’) 
a’ 


ce’ =@)+e, — 


a(Lx+myy). (10) 


a’ {= aa eee PB mz — Fae (x+my)} = 


ac a’*¢ 
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Since xz/r? and yz/r® are surface harmonics of degree 2, this can be written 
r 
2(@é9+e,—e’ )rY,=2e, (x- “) rY. (11) 


3. Gravitational work function of the core.—In a body slightly deformed from 
a sphere, so that the surface of uniform density and mean radius 7, is 


r=1r(1+¢,Y9), (12) 
where ¢, may depend on 7,, the gravitational potential is 


@ /1a®_ 1a’5¢’ gt, Sete eg 
U= af | "aa (a5 eid Ys) da +4nf | Pag (ta*+ ar%6 Y,)da 


= — 4nf¥o ["p {= - = ~ S73da’ d (a'54' )} a’ + $nfY. 2) "nth da 


+ function of 1. (13) 


Here we take 


#r:)=2¢, (1- 3) (14) 


in the core, ¢=0 in the shell, so that the upper limit in the second integral is 
the mean radius of the core. Then the part depending on Y, is 


2 rt oe ry 
8ufe, Ya — 4pori (: i 3) + 76P171 + 35P1 7a — Bp, 4}. (15) 


Associate with Y, a number « between o and 1, and consider a redistribution of 
mass such that « increases by 5a. The matter at a place is replaced by matter 
from a place deeper by 7,4, Y.5«, and its density is accordingly increased by 


dp T\ 3. 
—nbi Ys 7 - 7 8 = s(x 3) 3 Y,5a. (16) 


The part of U depending on r, yer does not affect the work function, since the 
total mass is unaltered. ‘Then the increase in the gravitational work function is 


‘aa te | r? a ax ré 
32nfada. ep, | | } Fd | {- Sport (: *¢ 1) + oP171 + 35P1 7a - Bea 
2 
1 
2 


r 
x (x- 1) dr, = — 1 nfo Saetp, (7 +m). ana?” (4P0— $p1) (17) 
and the work function for « =1 is 
. 512 
W,= 
1 3.15-693 
since mean gravity at the surface of the core is 


: 16 A 
a fe? — 35,)(I2+ m2) = — — 1 ge? 1, (2 +m?), (18 
‘feip,a"(po =P1)( 1 1) 693 a, 81 pa 1 1), (13) 
1 = $nfa(py— 3p). 

We can write (18) as W, = — A,w*e" (17 + mi), (20) 
and e” is another small number. . 

4. Work function due to the external field.—We take the external potential to 
be (c,x+c¢,y)z. The corresponding work function in the displaced position is 
found to be 


a c 
(Ag3 — Ay;)(,/ + cgm) + (: As33 — - An) (cl, + Cgm) (21) 
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and adding the contribution for rotation of the shell we have 

Wy = —(C— A)(eql + cgm) + €’C,(c,l, + comp). (22) 
The Lagrangian function is now the sum of T from (5) (the terms in /? + m? being 
omitted), W, and W,. 

5. Values of e,, <’, «”.—As we are attempting only a rough treatment it seemed 
that Bullen’s original solution* would be adequate. But on inspection of this 
solution a slight numerical inconsistency was found. The solution for p is not 
far from the form adopted here, and as, for uniform composition, p,/p) depends 
mostly on the velocity of longitudinal waves in the core, it can be regarded as well 
determined. All Bullen’s trial values make it about 0-190. This makes J/Ma? 
very near 0-390, whereas he gets 0-395. _ Using 0-390 and interpolating between 
his solutions for the shell I get the density at depth 400 km. = 4-091 g./cm.', 


M(core) = 189-7 x 10% g,, 
I(core) =8-91 x 10® g./cm.?. 

n for the core boundary is now 0-076, as against Bullen’s 0-04, and 

€,/€9 = 0-038. (24) 

The adjustment of course does not affect his determination of the ellipticity 


correction in seismology to the accuracy needed, since an accuracy of 1 part in 10 
for this is quite adequate. With e+e, =0-00260 this gives 


(23) 


e’ =0-0000014. (25) 
Also e” =2°8 x 10-8. (26) 
The smallness of the numbers ¢’, «” does not prove that their effects are negligible, 
because small divisions occur in the solution. 

These estimates ignore the inner core. This has little effect on J, but may 
increase M by something like 3 per cent. This will multiply «’ by about 2 and 
e” by about 4. 

6. Free periods.—We put 

l+im=(, lL, +im,=¢,, C,+1¢,=k, (27) 
with a time factor e. Then the equations of motion are equivalent to 
(y+ w){Ay —w(C— A)}0 + (y+o)(Ary +€'A,w)o, =(C— A)k, (28) 
(y+w)(Ayy +e’ Aw) + {Ayy* + wyC\(I + 2€’) — 2€” Aw}, = —e’ CR. 


There are two free periods with y/w small and two with it near —1. For the first 
pair, treat (C—A)/C, (C,—A,)/C, and y/w as small quantities of the first order 
and neglect their squares and their products with «’ ande”. Then to this order 
the solutions are 
Y _ a’, C,(C — A)+2AA,e”" —2€'’A? (29) 
w AC, - A? 
For the homogeneous core there is a zero solution, corresponding to a free steady 
displacement of the core with respect to the shell. We now find that gravity 
converts this into a free oscillation with a period of about 5 x 10° years. The 
second solution is the Eulerian nutation. It differs from its value for the 
homogeneous core by about I part in 10‘; the change is therefore negligible. 





* K. E. Bullen, M.N., Geophys. Suppl., 3, 395-401, 1936. 
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As expected, one period makes y/w exactly —1. The other (corresponding 
to a resonance noted in Paper II) makes 
A ei dale Bae C(C,— Ay) +2(e' +") AT (30) 
w CA, — Aj 
It is a retrograde motion in about 350 days, and the corrections again amount to 
about I part in 10, 

7. Forced motions.—Except for periods very near the free periods, since the 
changes in the latter are so small, we can drop «’ and e” on the left sides of the 
equations of motion. The new feature is that the external field has a direct 
effect on the motion of the core relative to the shell. ‘There are two interesting 
cases; if y= —w-+n, m/w is small in both. In the first, the 19-yearly nutation, 
it is about — gy and therefore small compared with (C—A)/C and (C,—A,)/C,. 
In the second, the fortnightly and semi-annual nutations, it is about +34 or + 4, 
and therefore larger than these two numbers. The correcting factor for the 
amplitude of { is, roughly, 





n A,C,e’ (31) 
w (C—A\(C,— A, + Ayn|w) 3 


As there are two small divisors it looks as if this correction might be appreciable ; 
but for the 19-yearly period it is about — 3 x 10-6 and for the fortnightly one about 
+4x10-5. The effects of the changes of the coefficients on the left are even 
smaller. 

The conclusion is that all effects of heterogeneity in the core are taken into 
account with sufficient accuracy if the core is taken as homogeneous and given 
its correct moments of inertia: Even if, as seemed possible in Section 5, the 
factors «’ and e” have been somewhat underestimated, this remains true. This 
is to be regarded simply as an analytical device, for a homogeneous core with the 
correct moment of inertia would have less than the correct mass. 

8. Outstanding discrepancy in the 19-yearly period.—This analysis justifies the 
use in Paper I of the correct moments of inertia of the core. The effect on the 
amplitude of the nutation, with a rigid shell, was there estimated as —6 parts in 
1000. In Paper II the data used made the moments of inertia of the core too 
large; but the important point now is that allowance for elasticity of the shell 
made a reduction from 9g parts to 6-6 parts in 1000. It seems probable, therefore, 
that with the correct treatment of the core and the same elasticity of the shell 
the effect would be about —4 parts in 1000. Again, in the solution for the 
variation of latitude I found ¢=0-515y, ¢ representing the elastic displacement 
at the surface and y that at the core boundary. For the nutation the ratio is 
0-725. This indicates that distortion near the base of the shell is relatively more 
important in the variation of latitude than in the nutation, and therefore that it 
would be appropriate to use a lower rigidity in the latter. It appears possible 
that with a full treatment of the actual shell the effect might be as small numerically 
as —3 parts in 1000, making the theoretical amplitude about 9”-200 against the 
rigid body value of 9”:227.  G.M. Clemence* has pointed out that in Spencer 
Jones’s observational determination of 9”:213 (and in Jackson’s previous one) 
the correction found was applied with the wrong sign, and that the data actually 
lead to 9"-2066. If these considerations are correct the observed value may agree 





: Re 


*G. M. Clemence, Astr. Fourn., 53, 179, 1948. 
33* 
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better with the calculated value for a fluid core and elastic shell than with that for 
a rigid Earth; at any rate they are near enough for a full treatment to be worth 
while. 

The fortnightly nutation has recently been determined from observations by 
H. R. Morgan. The amplitude for a rigid Earth would be 0-088. Morgan * 
gets 0”-098, revising a preliminary determination that gave 0”-078. In this range 
a fluid core with a rigid shell would produce an increase of about 3, so that 
Morgan’s revised value is in good agreement. 

Morgan has also made an analysis for the 19-yearly period, introducing a 
possible phase shift as an unknown ; it is small and doubtfully significant. I no 
longer think that imperfections of elasticity can produce a noticeable phase shift. 
Since the whole effect of the departures from perfect rigidity is expressed by a 
change of the coefficient by about I part in 400, and at the worst dissipation could 
only shift the phase of this part by 90°, the maximum phase shift for the nutation 
as a whole would be about 0°-2._ The relevant tidal component has speed —w+n 
relative to the Earth, and is therefore the lunar K,, a diurnal, not along-period, tide. 
For the semi-diurnal tide we can get the order of magnitude of the phase shift in the 
tide itself by attributing the whole of the tidal friction to the bodily tides, and it 
would be well under1°. For the diurnal tide the shift will be of the same order of 
magnitude, and it is therefore unlikely that the resulting phase shift in the nutation 
will be over 0°*002. 

160 Huntingdon Road, 
Cambridge: 
1950 May 28. 


* Private communication. 
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Summary 


A method is described for the routine calculation of the aberration day 
numbers, C and D, from the components of the Earth’s velocity referred 
to the centre of gravity of the solar system. The presence of lunar terms in 
these day numbers gives rise to difficulties in checking the apparent places 
of stars at ten-day intervals, and two alternative methods are suggested for 
dealing with this problem. 





1. Introduction.—At the meeting of the International Astronomical Union 
in Ziirich in 1948, a number of questions concerning the volume Apparent Places 
of Fundamental Stars was discussed by Commission 4 (Ephemerides). The 
increasing need for greater accuracy and for consistency in the calculation of 
these apparent places led to the adoption of a recommendation :— 


«|. , that the Day Numbers C and D should be computed directly 
from the motion of the Earth referred to a suitable fixed frame of reference 
(instead of from multiples of the sine and cosine of the Sun’s longitude 
as at present);.... It is further recommended that the terms in the 
aberration depending on the eccentricity of the Earth’s orbit should 
continue to be omitted.” 


The correction for aberration, as computed at present in terms of the Sun’s 
longitude, is not the true value; apart from other approximations it does not 
include the constant terms due to the eccentricity of the Earth’s orbit. The true 
value may be computed from the ratios x'/V, y'/V, 2'/V, where x’, y’, 2’ are the 
components of the Earth’s velocity referred to the centre of gravity of the solar 
system, and V is the velocity of light; it includes not only the effect of lunar and 
planetary perturbations, but also the constant terms already referred to. These 
terms, which are constant for any given star, and which have always been omitted 
from the aberration corrections, may be considered as incorporated in the star’s 
mean place. The removal of these constant terms, as recommended by the 
I.A.U., is clearly necessary, since their inclusion would necessitate a revision 
of the mean places of all stars. 

These terms assume considerable importance in the case of photographic 
observations of planets and comets. As Cowell pointed out (1) as long ago as 
1900, the position of the planet is antedated by the amount of its light-time and 
this affords a complete correction for aberration, including the terms depending 
on the eccentricity. The mean places of the stars to which the planet is referred 
already incorporate these terms, and an error may thus be produced which can 
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amount to xe=0"-342. This appreciable source of error is removed in recent 
ephemerides, such as that of Pluto in the Nautical Almanac and American 
Ephemeris, which give astrometric places, i. e. mean places to which the constant 
terms have been applied. This type of ephemeris, which is immediately 
comparable with photographic observations, will also be given in future in these 
publications for the four minor planets Ceres, Pallas, Juno and Vesta. 

2. Formulae for day numbers C and D.—The method of computing aberration 
corrections from the components of the Earth’s velocity was first mentioned by 
Turner (2) in 1909, and the problem was discussed in greater detail by 
Plummer (3), who later (4) gave the first table of accurate values of C and D for 
the year 1910. In this early work the velocities were referred to the Sun as centre, 
while the effect of the planets Jupiter and Saturn was introduced by simple 
approximations. In routine calculations it is better to use the published values 
of the planetary coordinates, and the day numbers may be computed from the 
expressions 


C+constant terrm=+y'/V, D+ constant term =—x'/V. (1) 


The value of V used in these equations must be derived from the accepted value 
of the aberration constant, 20"-47; introducing the expressions for the constant 
terms 


C= +1189"-80 (y’ + Ke cos 7 cos €), = —1189"-80 (x’ + Ke sin 7), (2) 


in which = is the longitude of perihelion, e the eccentricity of the Earth’s orbit, 
¢ the obliquity of the ecliptic, x’ and y’ are measured in astronomical units per day, 
and K is the aberration constant divided by the factor 1189-80. It is to be noted 
that this factor is not consistent with modern values for the velocity of light, which 
would require a figure larger by about two units. Any future revision of the 
fundamental constant « will necessitate a revision of the factor in (2). The 
constant terms may be evaluated (for equinox 1950-0) from the numerical 
quantities 


e =0:016730, m =282°-081, € =23°-446, 
leading to 
Ke cos 7 cos «= + 553.107’, Ke sin 7 = — 2815. 107’. (3) 


These terms are not, in fact, constant, being subject to small changes in e and e, 
but the variations have no significant effect over long periods of time. 
Equations (2) may now be written 


C= +1189"-80 (y’ +553). 107’, D = —1189"-80 (x’—2815).10-7, (4) 


in which x’ and y’ are to be referred to the centre of gravity of the solar system 
and are in units of the seventh decimal of an astronomical unit per day. Now if 


X, Y, Z are the solar coordinates referred to the Earth, 
Xn» Yn» Z, are the planetary coordinates referred to the Sun, 


then the coordinates of the Earth, referred to the centre of gravity of the whole 
system, will be 


X,= -X- 2 mx,/(I+ Xm), (5) 





No. 5, 1950 Stellar aberration 469 


where m represents planetary masses in terms of the Sun’s mass, and similar 
expressions hold for Y, and Z,. 

The values to be adopted for the masses of the planets are not critical, and 
in the present paper the values given in the volumes of Planetary Coordinates 
are used without change. The maximum effect of a planet on x’ or y’ due to the 
use of (5) is, in the usual notation, kman=km/,/a, or, in units of the seventh 
decimal, 

Mercury 0-04 Jupiter 72:0 
Venus 0°50 Saturn 15°9 
Earth 0°52 Uranus 1-7 
Mars 0:04 Neptune 1-6 


Since one unit in these values will only affect the fourth decimal of C or D, 
it is clear that Jupiter, Saturn, Uranus and Neptune alone need be considered. 
The denominator of (5), however, must contain the masses of all the planets, so 
that this equation may be written in the form 


X,=—-X-—=m'x,, (6) 


where m’ is now used to represent the planetary masses divided by the total 
mass I:001345. The “ reduced’ masses to be used in this equation are thus 


Jupiter 0-000953 Uranus 000044 
Saturn 0-000285 Neptune 0-000052 


The velocity components of the Earth, including all the planetary and lunar 
perturbations, may thus be obtained from the differences of X,, Y, and Z,, and 
the aberration corrections are computed then from the usual equations: 


‘ 


2"; 
cos §6 Ax = — =sin a+ ~, cos a, 


V V 
; (7) 


x’. f e 2’ 
Aé = — ~ sin 8 cosa—~2-sin8 sin a+ —cos 8, 


V V V 
Practical methods for computing the corrections in this form have been 
discussed by Zagar (5), whose tables for 1940 and 1941 (6 and 7) give three day 


numbers for use with the second of equations (7). The conversion of this 
equation into the usual form 


Ab = Ce’ + Dd’ 


may, however, be accomplished without any appreciable loss of accuracy. 
There are three sources of error in replacing 2’ in (7) by y’ tane: firstly, the 
lunar and planetary perturbations in the Sun’s latitude, giving rise to a maximum 
error of about 10 in the seventh decimal, corresponding to 0”-0012 in the 
declination ; secondly, the secular rotation of the ecliptic, giving rise to an error 
of about 4 (number of years from 1950), or 0”-oor in the declination in about 
20 years; finally, the inclinations of the planetary orbits whose contribution, as 
shown below, is negligible. With the usual notation, for a particular planet, 


P A : : k , 
y= Three Aen 2,= 7g (Becos E— A, sin £), (8) 
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from which it follows that z/,—y/, tane cannot exceed 


ky/a ree 
—~— sec € sin ?, (9) 
where 7 is the inclination to the ecliptic. The error is thus less than two units of the 
seventh decimal for Jupiter and even smaller for the other planets. 
The desired standard of accuracy may therefore be attained with the usual 
formulae 


Aa=Cce+Dd,  A8=Cc'+Dd’. (x0) 


3. Calculation of day numbers C and D.—The computation is divided 
naturally into two steps as follows: 

(a) The daily differences of the Sun’s X and Y lead to velocity components 
X’ and Y’; taking third differences into account, and using central-difference 
notation, 


X’ =p(3—483)X,  Y’=p(3—489)Y. (11) 


The second term in these expressions cafinot exceed 10 units of the seventh 
decimal, and may be considered as determinable from the Sun’s ephemeris for 
any year. ‘These small terms may therefore be thrown, together with the 
constants of equations (3) and (4), into the second step, involving the reduction 
to the barycentre. 

(b) The second step, like the first, involves mechanical differentiation. The 
quantities Xm'x,, and Um'y,, which form the reduction to the barycentre in (6), 
are obtained from the heliocentric coordinates x, and y, of the planets (taken, 


for example, from Planetary Coordinates) at an interval of 100 days. These are 
multiplied by their respective “‘ reduced” masses, the sums of these quantities 
are tabulated for each date, and from the differences the daily velocities are 
found at 50-day intervals. Denoting these values by (mx)’ and (my)’, and adding 
the small term in (11) and the constant terms in (4), the quantities 


AY’ =(my)'— §u3°¥—553, AX’ =(mx)' — 4n8°X — 2815 (12) 


may be formed at intervals of 50 days. The day numbers C and D can then be 
formed from the expressions 


C = —1189":80(u8Y+AY’), D= + 1189"-80 (pdX + AX’), (13) 


which are now in a form suitable for machine working. In the case of D, for 
example, daily values of 2(u5X + AX’) are formed in one operation on a National 
machine by integration from second differences of X, with simultaneous 
interpolation of AX’ to fiftieths. Multiplication by 594-90 then gives the required 
values of D, three decimals being retained in the final values of both C and D. 

4. Lunar terms in aberration corrections.—Reference has already been made 
to the fact that C and D computed in this way contain small terms due to lunar 
perturbations. ‘The presence of these terms causes no difficulty in a daily 
ephemeris, but the calculation of the apparent places of stars at every tenth transit 
(such as in Apparent Places of Fundamental Stars) will be difficult to check, 
since the higher differences fail to run smoothly. Occasional fourth differences 
exceeding 100 in units of the third decimal will be found in the ten-day values of 
C and D. However, it is possible to estimate their value at any given date with 
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some accuracy; the effect of lunar perturbations on the Earth’s velocity may 
be shown to be approximately 


AX’ = + mvp sinG AY’ = — mvqcos(, (14) 


where m is the mass of the Moon (Earth=1), ( is the Moon’s mean longitude 
and wv its mean velocity. The factors p and q are functions of the inclination and 
longitude of the node of the Moon’s orbit, and have a range of values from 


0-88 to 1-00. These equations may be written in the form of variations in 
C and D: 


AC = —0”-0086 qcos (, AD = —0"-0086 psin (. (15) 


They may be expressed, with sufficient accuracy, in units of the third decimal 
of a second of arc in the form 


AC =-—9cos(, AD= —gsin(. (16) 
Thus the variations in C and D are actually of the type 


Asin(B+ 4), 


where @ is the mean motion of the Moon in 10 sidereal days, namely 131°-4. Now 
the fourth differences of this expression are 


16A sin‘ $6 sin (B + 8), 


being 16sin‘@/2, or approximately 11, times as great as the original function: 
their amplitude will thus be about 100. The variations in the apparent place 
caused by the lunar terms are thus 


Aa =cAC + dAD, AS =c’'AC + d’AD. (17) 
We may combine the first of these equations with (16) to give 
Aa = — 08-0006 sec cos ((—«). (18) 


The maximum value of this expression is 0*-o012 for a star of declination 60° 
and 08-0035 for a declination of 80°. For the declination, the second equation 
of (17) may be put in a similar form by the transformation 


c’=rcoss, d'=rsins, (19) 
so that 


As = —0”-009r cos ((—s), (20) 


in which r cannot exceed 1-09, so that Ad has a maximum value of 0”-o10. 

The lunar terms, although small, thus affect the last figure retained in an 
ephemeris place and for the most accurate work must clearly be included. 
Interpolation of a ten-day ephemeris by ordinary methods, using second 
differences, provides adequate accuracy, in spite of the fact that the phase of the 
lunar terms changes by 131°-4 in the ten-day interval. For many purposes the 
inclusion of these small terms may be considered unnecessary, and it is then 
clearly simpler to remove the lunar terms given in (16) from C and D before 
computing the apparent place. The ‘“‘ smoothed’’ values of C and D so obtained 
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will give star places which can be checked in the usual way by differencing, while 
if the highest accuracy is required, the lunar terms can be added to the computed 
place by means of the corrections (18) and (20). 

Alternatively, the computer may prefer to compute apparent places directly 
by means of the correct values of C and D, and use multiples of (18) and (20) to 
check the fourth differences within the limits normally imposed by rounding-off 


errors. A suggested method of tabulating the lunar terms and of employing these 
two methods is given below. 


TABLe I 
Day numbers C and D, 1956, for every tenth sidereal day 





True values ** Smoothed ”’ values 





54C ; d!C 





o7 | — 2°984 -¢ +20°197 — 2°990 + 20°204 
6°213 19°318 | 6214 - 19°309 
9°249 17°808 9°241 6 17°813 
30°6 11°966 +74 15°757 - 11°975 26 15°759 
9°6 14°330 —S5I 13°229 + 14°326 7 13°221 


196 | —16°233 +10°277 —16'230 —19 +10°285 
29°6 17°618 +63 77047) + 17627. —33 7°044 
10°5 18°485 —89 3609 4 18477 + 1 3°605 
20°5 18771 — 2 + 0068 - 18-773 —26 + 0077 
30°5 18502 +36 — 3°429 ‘ 18-507 —20 — 3436 —I15 


951-1770 -97 ~— GOe7 —17697 — 1 — 6826 — 6 
19°4 16°375 +48 10:006 — 16°381 —22 10000 — 6 
29°4 14598 — 3 12°862 + 149598 — 7 12871 —23 
9°4 12.416 —77 15°364 12°409 - 2 15°358 — 7 
19°3 9°873 +81 17°405 9882 —22 17°403 —12 


293 | — 77090 —43 —18:948 + — pols +236 ~—20055 —87 
8-3 4107 —34 19°984 4104 — 9 19975 —14 
18-3 | — 1°007 +79 20°437 + 33] — rors —10 20°441 — 9 
28-2 + 2°093 —70 —20°341 + 22| + 2-101 +20 —20°345 —2!I 











5. Checking the lunar terms.—The first of these methods will necessarily 
involve the publication of ‘ smoothed” values of C and D which, by way of 
illustration, are given in Table I for part of the year 1956. The data are given 
at o sidereal time for every tenth sidereal day, and the “ smoothed ”’ values 
have been obtained by the application of (16) to the true values. The Moon’s 
mean longitude has been taken as 136°-6 on the first date, and increments of 
131°-4 have been added for each subsequent entry. It would probably be simpler 
in practice to use a series of terms such as those given in Tables II and III, since 
it is not necessary to remove completely the oscillations in the fourth differences. 

The star positions may then be computed in the normal manner, using the 
“‘ smoothed ” values of C and D; the aberration corrections computed in this 
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way are shown in the first columns of Table II for a typical star (FK3 No. 1397), 
one extra figure being retained in order to reduce the rounding-off errors. The 
lunar terms to be added to the smoothed values to give the final values are 
tabulated in the next columns of Table II and are taken from Tables III and IV. 
Table III gives the values of (15), (18) or (20) for four different amplitudes and 
for increments of 132° in the angle. The four amplitudes chosen (0-9, 3, 6 and 9) 
are sufficient to cover all possible cases, and the use of 132° to represent the 
motion in ten days is sufficiently close to the truth to give ample accuracy, while 
having the advantage of giving a series recurring after 30 terms; the unit is the 
third decimal in all cases. Table IV gives the series and number of the term in 
Table III to be used to represent the correction on the first date of the year’s 
ephemeris. Subsequent corrections at intervals of ten sidereal days are then 
given by succeeding terms of the same series, taken in the order printed, and 
repeating from the first term when the thirtieth term is reached. Table IV must 
not be interpolated; it is entered with the nearest hour of R.A. and the nearest 
tenth degree of declination, the intervals having been chosen so that there is very 
little divergence between neighbouring terms. It will be noticed that series P 
is to be used for the R.A. corrections for stars with declinations less than 60°; 
for stars nearer to the poles it is customary to drop the third decimal in the R.A. 
and series Q is to be used. The declination corrections have a zero value at 
65, +23° or 18", —23° and have a maximum value of 9°8, but the four ranges 
of amplitude given have been selected to cover these values with a minimum of 
error. Since the Moon’s longitude changes by 132°-7 in one Julian year, these 
tables may be used for future years bv increasing the numbers of Table IV by 
unity for each year after 1956. This process may continue for some fifteen years 
without the necessity of a revision. 

For star No. 1397 (R.A. 15" 05™; Dec. +54° 45’) it is seen that the R.A. 
corrections have P13 for their initial term, while those in declination commence 
with S15. The application of these corrections, followed by rounding off, gives 
the final values in Table II]. These may be compared with the true values in 
the following columns, the latter having been computed directly from the true 
values of C and D in Table I. 

The second method is also illustrated in Table II, where the computed 
aberration corrections are differenced. The oscillations in the fourth difference 
are compared directly with the terms of Table III, commencing as before with 
P13 and S15, but noting that the terms now represent units of the second decimal, 
since they refer to the fourth differences, approximately 10 times the function 
values. The example illustrates the order of agreement which is obtained, and 
which usually comes well within the accepted limits of six units in the fourth 
difference ; the signs of the differences are indicated with some certainty. 

The error involved in the use of these tables is very small. The discrepancy 
between the true motion of 131°-4 in ten days, and the assumed value of 132°, 
may give rise to a divergence from the correct value of the longitude at the end 
of the year of about 23°. By commencing the series with a value of 127°, instead 
of the true value of 136°-6, the errors are distributed over the year and 
do not exceed 0*-0003 in R.A. (at 60° declination) or 0”-004 in declination; 
they can be ignored. The errors due to the use of an approximate amplitude are 
also negligible. 
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Taste III 
Lunar terms in aberration corrections 





P Q R S Term P Q R 





16 —o's5 —4 
17 —o'2 —I 
18 -o8 +5 
19 —o'9 —6 
20 toq = 4 +3 


o’s 
Loz 
—o8 
+o'9 
—o-4 
—o'4 
+or9 
—o8 
o*2 
+o'5 


| 
| 


| 


21 +O'4 + +2 
22 —o'9 —6 
23 +o'8 +5 
24 —o'2 —I —1 
25 —o'5 —2 -3 
26 +o"9 +3 +6 
+O"7 27 ry OE 2 Hae 
oo ° 28 roe) ° ° 
OF... 8. ee SOO top... +0.. 8 
teg +3 6. +6): ei SES 


oe ON AM AW DN 
+1 ++14 +1 +4 
OM DW HU DN |B AU HS 


| 


—o'9 


i ee sie. ae Se at eB, oe de 
ONO NN WO FO ONWN 











The entries in the above table are in units of the third decimal. 


TABLE IV 


First terms in right ascension and declination for 1956 





Right Declination 
ascen- 
sion ‘ ’ 30 40 50 60 
PorQ 
12 Qi2 Q28 Rsy Ree RS S 6 
I Qi2 Q28 R17 R6_ R25 S14 S14 
9 Qi2 Q28 Q6 Ri4 Rri4q n.3 By 
28 Qi2 Q17 Qz5 Q3 R22 R22 S11 
17 Qi2 Q17 Gs Ges  Krs R30 =©6@R3o 
6 Qi2 P28 Q30 Qig Qri9 RS RS S88 
14 Qi2 Pi2 ray. Qazy Gaz R27 R2z S827 
3 Qi2 P26 Q13 Q24 Q5 as 2S RS 
22 Qi2 Q7 Q10 Q2 Ri R24 R24 S24 
II Qi2 Q26 Q2z9 Qari Rar Re Ba tes 
19 Qi2 Q26 Q18 R2zg Rio Ra: Sa: Sar 
8 Qi2 Q26 R7 R18 Rag R29 S10 Sro 
27 Qi2 Ris R26 R7_ R18 $18 S29 S29 
16 Qi2 R15 R15 R26 S7 9 BF Be 
24 Qi2 R 4 Ris S15 S26 $26 S26 S26 
13 Qi2 R 4 R4 S44 Srs5 Szs Sxs Sis 
2 Qi2 R23 Raz. Sas Sa S34. 84.83 
21 Qi2 Ri2 Ra3 S23. S823 S23 S23 S23 
29 Qi2 Ri2 Ria S:12 S12 S:a Ssa S12 
18 Qi2 R 1 Riz 3 Str 5 ea ee 


| Qi2 R 1 R20 S20 S20 Sao Sao So 
26 Qi2 R20 R20 Sg S9Q S9 So9 S828 
4 Q12 R20 R2o Rg S28 S28 S28 Sr7z7 8217 
23 Qi2 R2o Rg R28 S28 S17 Sxy 86 36 23 
For corrections in right ascension, use series P for all stars having declinations less than 
60°, and series Q for stars with declinations from 60° to 80°. 
For negative declinations, use R.A.-+ 12" as argument for the corrections in declination. 
For other years, increase the numbers in the above table by unity for each year after 1956. 
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The second of the two methods described is more accurate and more direct, 
and should be used for systematic computation of apparent places of stars unless 
the method of calculation used relies on the smoothness of the differences. In this 
case, the first method offers a simple alternative. 


H.M. Nautical Almanac Office, 
The Royal Greenwich Observatory, 
Herstmonceux Castle, 
Sussex : 
1950 October 4. 
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ON THE EXPULSION OF CORPUSCULAR STREAMS BY SOLAR 
FLARES 


F. D. Kahn 
(Received 1950 July 8) 


Summary 


A value is found for the minimum momentum density in a corpuscular 
stream which has been emitted from a flare and which causes a terrestrial 
magnetic storm. It is shown that radiation pressure alone will always be 
insufficient to expel such a stream from the Sun. 





1. Introduction.—Some twenty years ago Chapman and Ferraro (1) developed 
the corpuscular theory of magnetic storms, which is now generally accepted. They 
showed that terrestrial magnetic disturbances should take place whenever a cloud of 
particles emitted by the Sun travels past the Earth. In the case of the irregularly 
occurring magnetic storms, which are strongly associated with solar flares (cf. 
Newton (2)), it must be supposed that the particles are emitted by the direct 
action of the flares themselves. From the observed average interval between the 
two phenomena it is commonly deduced that the particles have speeds of between 
1000 and 1600 km./sec. 

Kiepenheuer (3) has proposed that the radiation pressure due to the flare may 
cause the expulsion of the cloud. This suggestion has also been subsequently 
made by Hoyle (4) and by the author (§). However, it will be shown in this paper 
that the momentum which can be communicated by radiation alone to the nearby 
particles will always be insufficient to expel a stream dense enough for a storm of 
observed intensity. 

2. The required momentum density of the stream.—To obtain a workable model of 
the expulsion it was assumed in a previous paper (5) that the particles are emitted 
radially from the surface of the flare, here taken to be the cap of a sphere of radius b. 
The flare is situated in the solar atmosphere and the emission continues from time 
t=o untilt=7,. Since storms always last much longer than the associated flares 
it was supposed that the velocity of emission is highest at the peak of the flare and 
then falls off gradually, thus giving a longitudinal expansion to the stream. With 
these assumptions the density near the front of a stream, t seconds after the flare 
peak, is given by 
b?N(0)7, 


"= oy— 0) 


whent>7,, 
b =hypothetical radius of flare (2 x 10°cm.), 
N(0) = rate of emission per unit area of flare, 
7, = duration of flare (10* sec.), 
U) = velocity of emission at flare peak (1-6 x 108 cm./sec.), 
v, =velocity of emission at end of flare (10° cm./sec.). 


(This formula was obtained in (5), Section 3.32. The values given in brackets 
in the list of symbols are typical, and will be adopted in the subsequent working.) 
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If mis the average mass of the particles, the momentum density near the front 
of the stream is given by 
b?N(0)mr, 


nae mg. = 
F=n(t)mvp = (vp —2,)¢3 ’ 


when t>7}. 


The momentum density at the emitting area is Fy = N(0)mv,, hence 


, 


_ Py  4x10F, 
 Up(Vp—2,)e2 _ ans 





When the storm reaches the Earth, t~10° sec., and 
F=4x10-! F,. 


Chapman and Bartels (6) give an estimate of the minimum value of F, sufficient 
for a magnetic storm. ‘They find that the least possible momentum density is that 
of a stream containing 10 H-atoms per c.c. and travelling at 1-6 x 10°cm./sec. 
Thus, if m, is the mass of a hydrogen atom, 


F > 10myv; = 4 x 10-7 dynes/sq. cm. 
and F, > 10° dynes/sq. cm. 


The rate at which momentum is given to the stream by the flare must then 
exceed 107 dynes/sq. cm. 

This is not a very large value, but nevertheless it will be shown to exceed the 
radiation pressure which the flare can exert. A much smaller value is required by 
Hoyle, who assumes that particles are emitted radially from the Sun rather than 
radially from the flare, and who neglects the longitudinal expansion of the stream. 
On this model the density on emission can be much lower, but, since a flare 
subtends only a small solid angle at the centre of the Sun, the resulting emission is 
confined within a very narrow cone. Now even a flare well removed from the 
centre of the Sun’s disk can cause a magnetic storm, so that we need to have cones of 
emission wider than is supposed on Hoyle’s model. 

3. The mechanism of radiation pressure.—T here are two ways in which radiation 
pressure can act on the particles : 

1. The atoms may absorb resonance radiation coming from the flare, and 
re-emit at random. 

2. The atoms may be ionized by radiation coming from the flare. 

In either case the absorption of a photon of wave-length A gives an impulse 
h/X away from the source. 

3.1. The resonance process.—Consider an atom with a resonance line at wave- 
length Aycm. Let the atom be exposed to a field of radiation contained within a 
solid angle Q at P and le: it have a velocity v away from the field. Light arriving at 
P within a small cone of solid angle dQ around the line RP (7 RPO =8@) can be 
absorbed by the particle if its wave-lengui is 


A=Ajy (x- = cos #) 
c 


or if AA=Ay~ cos 8, 


and the particle will be raised to an excited state. Each photon coming along RP 
has momentum //A, where A is Planck’s constant, and the component of this along 
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OP will be A cos 6/A. If the emitting region is symmetrical about OP, the side- 
ways components will cancel out. 

The atoms are accelerated to a final speed v, and the width of the spectral 
band in which radiation coming from the direction RP can be absorbed is given by 


The increase in outward momentum which is given by each absorption is 
approximately hcos@/Ay <h/Ay, and in order that a group of atoms may gain 
momentum p, the total number of absorptions and re-emissions must exceed 
Poro/A. Now the rate of gain of momentum of the stream near the flare is at 
least 10° dynes/sq.cm.; if this is entirely due to radiation pressure the rate of 
absorptions and re-emissions above a unit area of flare has to be greater than 
10°A,/h per sq. cm. per sec. 





y 


Fic. 1.—The radiation field at P. 


0 
Q 


Let J, photons/cm./sq.cm./sec. be the intensity of the radiation emitted by 
the flare near A, cm. wave-length. The maximum number of absorptions over 
a unit area of the flare within a band of width AA, is then 

T,AAy =I, ms photons/sq. cm./sec. 
Hence, if radiation pressure is the sole cause of the expulsion, 


Ao 


2 10°F 


, 


1, et 
Cc 


c 
or I,> 10° 7 = 3 x 10%! photons/cm./sq. cm./sec., 
Uo ; 


if we adopt the values 

h=6-6 x 10-*’ erg.sec., c=3x10!cm./sec., Up =1-6 x 10®cm./sec, 
Where the absorbing atoms have a resonance doublet such as the H and K lines 
of Call the minimum intensity need have only half this value. 

It will be shown that J, is much greater than the intensity of radiation to be 
expected from the flare, whichever group of atoms we choose to consider. Thus 
Partridge (7) has found that the intensity near the H and K lines is of the same 

34 
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order as black-body radiation at 6000 deg. K.—it may be up to five times as intense. 
In full radiation the intensity 3 is given by 


H(v) dv = ll exp(—hv/kT) dv, 


or B(A) dA = ced exp(—/Ac/ART) da, 


or 3(A) =7-6 x 10 photons/cm./sq. cm./sec. 
This is only 5 x 10-5 of the required least value of the intensity, and, at best, 
H and K radiation pressure fails by a factor of 2-5 x 10~*. In contrast, this 
pressure is strong enough on Hoyle’s model. 
Although there are, of course, no experimental data about radiation at 

A =1215A. it may be inferred that the Lz line is also insufficiently strong for the 
expulsion. To be effective the rate of absorption of L« photons in a band of 
width dA above the flare must exceed 3 x 10%! dA photons/sq. cm./sec. and so the 
rate of absorption by H atoms with velocities between v and v+dv exceeds 
3 x 108A, dv/c, where 

a = ad Ayg=I215A 

} sis c ’ 0 an ° 
Each absorption raises an H atom to the 2p state whose natural lifetime is 
1'6x10-* sec. It follows that the number of H(2p) atoms over a unit area of 
flare with velocities in the range (v, v + dv) will be 

3 x 10%! x 1-6 x 10° *Ay dv/c ~ 5 x 107*Ay dv/e, 
unless there is some process which de-populates the excited state at a rate 
comparable with or faster than that of spontaneous transitions to the ground 
state. Each of the 2p atoms can absorb in Ha light (A, =65634.) where the 
absorption cross-section is, after Unsold (8), 
_ TeX eo. ~il om 8 
ig as f=2:1 x 10-7! cm.?. 
Here e=electronic charge, m=electronic mass, c=velocity of light, 


f =oscillator strength of the transition, A, = wave-length of Ha. 


The atoms with velocities between v and v+ dv absorb over a band of width 

dX =, dv/c and the optical depth of the cloud for Ha light becomes, at least, 
a, x 5 x 1072 x =19, (Ao/Ay = 3). 
1 

The cloud of corpuscles is likely to have a great optical depth for La light 
and an La photon may be absorbed and re-emitted many times within the stream, 
although only the first absorption contributes to the momentum. The density 
of H(2p) atoms should thus be larger, and the estimate for the Ha optical depth 
should be too low. 

With such a great optical depth the cloud should cause a far greater asymmetry 
in the Ha emission of the flare than is commonly found. A rough comparison 
of the red and blue wings of the Ha line, as observed by Ellison (9), shows that the 
excess absorption on the blue wing is about 15 per cent. This is due to a cloud, 
with an optical depth of the order 0-15, receding from the flare. If the radiation 
pressure theory is to be maintained we must find a process which reduces the 
estimated density of H(2p) atoms by a factor of 10°. 

The cloud is too tenuous for collisional de-activation to be important, but it 
might be thought that ionization of H(2p) atoms by radiation of the Balmer 
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continuum is fast enough. According to Unsdld (xo) the total oscillator strength 
of the transition 2p—> continuum is f=0-19, and the total absorption coefficient 
for light of wave-length A < 3647 a. becomes «, =(7e*/mc)f (it is more convenient 
here to use the absorption coefficient per unit frequency interval) or 


a, =5°0 x 107%, 
In order that the removal of 2p atoms may be sufficiently fast we require 


6 x 10!° ionizations per H(2p) atom per second; if J, is the maximum intensity 
of the flare beyond the Balmer limit «,J, >6 x 10”, 


or J, >1-2 x 1018 photons per unit frequency interval/sq. cm./sec. 
’ p p q y q 


This is about 4 x 105 times the intensity of a black body at 6000 deg. K. at the 
Balmer limit, and is by many orders higher than the emission of a flare in this 
part of the spectrum. The low observed density of H atoms in the 2p state shows, 
therefore, that Z« radiation pressure cannot be the cause of the expulsion. 

3.2. The ionizing process.—Pressure is exerted by ionizing radiation in much 
the same way as by resonance radiation, but an essential difference between the 
two processes is that the ions here formed have to recombine with electrons 
before they can again be acted on by photons. The rate of recombination 
enables us to fix an upper limit to the effectiveness of radiation pressure. Here 
it is best to investigate the acceleration of the cloud by the radiation, and the 
previous model of the emission must be abandoned. 

To take the best possible case, consider a cloud of pure hydrogen. Let a 
typical unit volume of it contain n(I—x) H atoms, nxH* ions and nx electrons. 
If A is the total recombination coefficient of H+ ions to H atoms, the rate of 
recombination becomes An?x?<An?. For equilibrium An*x? is also the maximum 
rate of ionization of H atoms, and since an ionization by light of wave-length A 
transmits, on the average, an impulse 4/A, the maximum rate of communication 
of momentum to the group becomes An*h/A. In a time dt an upper limit to the 
increase dv in the velocity of the group is given by 

m,n dv < = dt 


A 


Anh 
or dv < —~ dt =3°3 x 107m dt, 
~ MyA 33 
where A=900A.=9 x1I0-* cm. (atypical value), 
A=7 x 10-% cm.3/sec. (cf. (§), Section 4.32, Table II). 

The acceleration produces a longitudinal expansion so that, while the cloud 
is near the flare, the particle density will decrease faster than by an inverse square 
law, or n< #2b?/r’. 

Here #2=density at the emitting area, )=radius of flare, r=distance from 
flare centre. 

The front of the cloud is accelerated to a final velocity vp in time 7,, while the 
acceleration decreases with increasing distance from the flare; if v is the velocity 
t sec. after the emission, 


if t <7}. 
At a time ¢ the group will be r—6 cm. from the surface of the flare, and so 


: j 
r—b= | vdt>,| 
0 J 


t¢ 


— t dt =4uot*/7,. 
071 


34* 
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dv< i dt 
MyA 


<3°3 x 10-1 $b? /r? dt 
< 33% 107° $25? 
~ (b+ 40et?/7,)? 
With the substitution tan 6 =t,/(v9/2br,) we find 


Ve tan= V (v97;/2d) 
dv <3°3 x10 © R i (=) [ cos? 6 dé, 
J0 0 


~0 
<3°3 x10" RN J (=) > 
0 


since tan™! 4/(v97,/2b)< 2/2. 
Insertion of numerical values gives 
Uy <4 x 108 RN. 
Now we require that vy =1-6 x 108 cm./sec., and this implies that 
N>R>4x 10%, 

Thus a cloud can be expelled by radiation of the Lyman continuum only if 
its initial density exceeds 4x10! particles per c.c. However intense the 
radiation may be, expulsion cannot take place in more tenuous clouds, since the 
rate of recombination could not keep pace with the required rate of photo- 
ionization. But the density of the solar atmosphere near the flare is much less 
than 4x10" particles per c.c. and this process also fails by some orders of 
magnitude. 

It seems hard to think of any other way in which radiation pressure could 
possibly be effective, and some entirely different mechanism of expulsion must 


be found. This should not prove too difficult, for the required excess pressure 
near the flare is only of the order of 10° dynes per sq. cm. 
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AN INVESTIGATION INTO THE POSSIBILITY OF OBSERVING 
STREAMS OF CORPUSCLES EMITTED BY SOLAR FLARES, II 


F. D. Kahn 
(Received 1950 July 8) 


Summary 
The ionization in a stream of solar corpuscles is re-calculated by means of 
two equations derived in a previous paper (M.N., 109, 324, 1949). The new 
values differ somewhat from those given there; in particular the Ca II content 
is found to be rather lower. This would make the spectroscopic detection of 
the stream even harder. 





1. Introduction.—The physical state of a cloud of corpuscles emitted by a 
solar flare has been discussed in a previous paper (x).* A simple model of the 
emission from a flare was used to calculate the variation of particle density in a 
typical stream. It was found that 


k+5 
n(t) = i (I, Section 3.32 (9)), 


where n(t)=density near the front of the stream at time ¢, 
t=lapse of time from instant of emission, 
1o* =density near the emitting area. 


‘This formula holds approximately when t>7,, the duration of the flare. 


The equations of ionization by solar radiation after the end of the flare were 
shown to be 


dx ‘I. Secti 
< = Bya(t—x)- Anz? yl, Section 4.33 (5)) (1) 


for hydrogen, and 


ES oe oe ; 
at = By\(I—x')—A’nx (I, Section 4.33 (6)) (2) 


for any other elements present. Here 


8 =dilution factor of solar radiation, 

x = proportional ionization of the element, 
Yi =1onization coefficient from ground state, 
A =total recombination coefficient, 
n,=electron density. 


The unprimed symbols x, y, and A refer to hydrogen, the primed symbols 
to any other element. 

Owing to the great abundance of hydrogen it is permissible to assume that 
the combined density of H atoms and H* ions is approximately equal to n, the 
particle density, and that n,~nx. When equation (1) has been solved to find n, 
the solution may be used in (2) for the state of some other, less abundant, 

* This paper will be referred to as I in the text. 
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elements. The ionization of sodium and of calcium was found to be of particular 
interest in Paper I, for NaI, Cal and Call have good resonance lines which may 
make the cloud detectable in absorption. 

Since some rather crude approximations were used in Paper I for the solution 
of equations (1) and (2) we shall give a better treatment of them here. 

2. The ionization of the cloud during the flare.—Before going on to solve the 
equations we must know the ionization of various elements in the stream when the 
flare ends. Observational evidence will help to find this. During the flare the 
ionization of any part of the stream will not tend to change. Near the emitting 
area the particle density in the stream will diminish as the inverse square of the 
distance from the flare (I, Section 3.31) in the same way as the intensity of the 
radiation received. Thus ionization and recombination will remain balanced 
at this time. 

Solar flares are known to emit the He resonance line (A = 584.) and are thought 
to be rich in radiation of the Lyman continuum, and thus the hydrogen will be 
highly ionized early on. In the later working the typical values assumed will 
range from 80 to 99 per cent ionization of H. 

Both NaI and CaI have low ionization potentials (y= 5-12 e.V. and 6-09 e.V., 
respectively) and may be expected to be almost entirely absent. CalII is harder 
to ionize; here observational evidence must be used to fix a lower limit to the 
ratio [CalII]:[CalIlI]. If there were much CalII present we should expect the 
H and K emission lines of the flare to show an asymmetry similar to that observed 
for the Ha line by Ellison (2). Partridge (3) has failed to find such an asymmetry 
for the H and K lines, and so there can be little absorption over the flare. 

Suppose that 10* is the density of the stream on emission, then the number 
of corpuscles in a unit column along the line of sight from the Earth to the flare 
will be MW, say, where 


a) 


b? : 
v~| 10* — dr = 10*b =2 x 10**9. 
b Tr 


Here b=radius of model flare =2 x 10° cm. 
On Strémgren’s composition (quoted in I, Section 4.1, Table 1) this contains 
a total number of Caions (Call or CaIII) equal to 


2 x 10*+9 x 1-6 x 10-§ = 3:2 x 10*+8, 


Since the ratio [Ca III]: [CaII]=x':1-—x’, the number of Call ions per unit 
column becomes 3:2 x 10*+3(1—x’). These ions may have a velocity spread 
of up to 1600 km./sec. owing to the longitudinal expansion of the cloud, and their 
absorption band will extend over a width of 224., or less. The number of Call 
ions per cm. spectral width in the column is then at least 1°5 x 10*+%(1—x’). 
Figures given by Unsdld (4) show that 10° Call ions per cm. width in a unit 
column cause about Io per cent absorption, and in the case of the H and K lines 
the depth must be of this order, or less, thus 


5 x 1o¥+10(7 — x’) < 10, 
(1—x') <7 x 10**, 


We find the following minimum values for the ionization of CaII while the 
flare is in progress. 
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TABLE I 
Minimum ionization of Ca II during a flare in some typical clouds 
x’ 
k=11 0°93 
k=12 0°993 
k=13 ©9993 

Now the anomalous absorption by Call ions should be strongest directly 
over the emitting area. During the flare no effect is visible here, by Partridge’s 
result, and it is useless to search for absorption on other parts of the Sun’s disk. 
A stream in which there is little recombination later will remain invisible even 
when the flare has ended. 

3. The solution of the equations for ionization after the flare. (a) The solution 
for hydrogen.—The flare ends when t =7, and at later times the front of the cloud 
will be several solar radii removed from the Sun. We may then take B = a?/4v2t?, 
where 

a=radius of Sun, vy,=velocity of the front of the cloud. 
Inserting this and the expression for n(t) into equation (I) gives 
dx  a*y,I-—x) 4x 10*t5Ax? 
dt 4u2t? “e ; 
We put 7'=t/7, and find 
dx  a*y,(I—x) 4x 10*t5Ax? 
dT = 4v3r, T? 7?T? 
A(i—x) Bx? 
= Aa ~ "s- (3) 


The initial condition is that x=x, when t=7,, or T=1. Adoption of the 
numerical values 








a=7x10"cm., vy=1-6x 10%cm./sec., 7, = 10% sec., 


¥1 =1-2 x 1073 sec! 


cf. I, Section 4.32, Table II. 


A=7 x 107% ¢m.3/sec. 
gives 
A=5:7x10-?, B=2-8x10"%, 

The cloud reaches the Earth when t~10°sec., or T~10?. 

There is a simple method of approximating to the value of x which allows 
some insight into the physics of the problem. With the initial condition that 
x=x, when T=1, the solution of (3) for any T>1 lies between the solutions 
Y, y of (4) and (5), where 

adY_ A(1-Y) 
tae 
oo. 
ee 

In other words x will be overestimated or underestimated according as 
photo-ionization or recombination is neglected. 


Y =x, when T=1, (4) 


y =x, when T=1. (5) 
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To prove that Y>«x we write (3) and (4) in the form 


2 
-A/T) e74AT Bx* |-air 
T? , 


aT © hn 


Sn(Ye~4t) = & ben e~AIT > Fn(xe-4. 


Thus Ye~4/7 > xe—4/7 + constant, for T>1. When T=1, Y=x=4,, so that, for 
T>1, Y>x. 

Similarly x>y. 

The equations (4) and (5) integrate to give 


Y=1—(1-—x,)e~40-¥), (6) 


r{Eo8(-#F 0 


The hydrogen will never be completely ionized or completely recombined, 
for Y steadily increases and tends to Y,=1—(1~—x,)e~4, while y steadily 
-3 

decreases and tends to y= (= + 3) . 
1 
If A is small it can be shown that (x—y) will remain small and at any given 
time x, the ionization, will lie close to y, the lower limit. To prove this we subtract 
equation (5) from equation (3) to find 


d , _ A(1—x) " 
aT (x —¥) es T2 £3: T3 (x J 
With x>y>y,, it follows that 


d A(I-Yea B 
a Are) — leet) 


(x-—y)< 


A(I—y., 2By., 
< Ape) ae (#— ~y). 





This may be re-written as 


- {(x¢ — y)e~BualT™) < A(I— yao) e~ Buell 


which gives on integration (with x—y=o when T'=1), 


- 1 A(I—-Yo) V7 > Py 
_— B ol T* — / _— 
(x—y)e~™ —. {erf V By. —erf Tf: 


But in the interval 0 <u <u, the maximum value of 
f(u) =e"'(erf uy—erfu) 
occurs foru=o. Applying this result to (8), with u* = By,/T?, 


x-y< <= AUN ext V By ‘ By a. 


V Byo 
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In the most unfavourable case to be considered B=2-8, x,=0-80 and 
Ya =0°37; this will give the largest value of (x—y). Here By,,.~1, and 


x-y<A(I — 0-37) erft 


x—y<0-027. 


This is a small margin of error, considering the nature of the problem. The 
table below gives the values of the ionization of hydrogen to be expected at two 
typical times, in clouds of different densities. 


TasLe II 
The ionization of hydrogen at two typical times after the end of the flare. 
t~10° sec. corresponds to the time when the stream reaches the Earth 
Distance 
from Sun k=10 k=1!1 
{cm.) 


Initial Time 
ionization (sec ) 


2x 10° 3°2 x 108 080 o°79 
080 
1°5 x 101% 0'80 0"79 


3°2 X 104 0°90 o'89 0°83 
1°5 x 10}% 0°90 : 0:80 
3°2 x 10" 0°99 . 090 0°50 


0°99 
10° 1*5 x 10}% 0°99 0°99 0°87 o*41 


With the formula n(t) = 4 x 10*+5/#3 for the density near the front of the stream 
we find the following values for ,, the electron density. 


Tas_e III 

The electron density at two typical times after the end of the flare for streams of various densities’ 

Distance 

from Sun 
(cm.) 


Initial Time 


Electron density (cm.~*) 
ionization (sec.) 


k=11 k=12 k=13 
2x 10% 3°2 xX 104 4:0 X 10° 3°7 X 107 2°2 X 108 
0:80 

10° 1*5 x 101% ° 3°2 X10 2°9 X 10? I*5 x 10° 
2x 10° 3°2 x 10! 4°5 x 10° 4°2 X 107 2°4X 10° | 
10° 1°5 x 1018 . 3°6X 10 3°2 X 10° 1°6X 108 | 
2x 108 3°2 X 1012 4°9 x 10° 4°5 X 107 2°5 x 10° i 
10° 1*5 x r0}* : 3°9 X10 3°5 X 10? 1‘6 xX 10° 
(b) The ionization of sodium and calcitum.—Some estimates for the ionization, 
and recombination coefficients of Nal and Call were given in Paper I 
(Section 4.32, Table II). Since then Ditchburn and Jutsum (§) have measured 


the continuous absorption coefficient for sodium. ‘Their experimental results 
show it to be only about 70 per cent of Rudkjobing’s calculated value (6), which 
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was assumed in Paper I. Further, rocket spectra taken at White Sands (7) 
have established that the intensity of solar radiation near the head of the sodium 
continuum (A~2400 A.) is only about one-tenth that of a black body at 6000 deg. K. 
In this way the estimate of y; for NaI is reduced to 2-6sec.-!. Even so this is a 
high value compared with that for hydrogen and, as before, we conclude that the 
sodium will remain highly ionized at all times. 
From data given by Bates and Massey (8) we find that for Cal 
A’ =2 x 107 cm.3/sec. and y; = 1-5 sec.~}, if the solar radiation is only Io per cent 
of full black-body intensity near the Cal continuum (A~2020A.) Once again, 
neutral calcium will be very rare, and it remains only to evaluate the relative 
densities of Call and Ca III. 
According to (2) the equation of ionization of calcium becomes 
dx’ _ ay\(1—*’) 
pre a A Tf ey. 2 
dt 407? si 
Now in Paper I we saw that y}=7-4 x 10~*sec.-1, and A’ =2 x 10° cm.3/sec. ; 
writing T'=t/r, gives 
=a Bey 4. =a a 2 x 107)37,n,x’. 
The first term on the right-hand side is small and so 
dx’ 
dT 
In the same way the equation of ionization of hydrogen may be written 
dx 
dT 
Division of (I) by (11) gives 





~—2x 107!37,n,x'. (10) 


~—7 x 107!37,n,x. (11) 


dx 2x’ 
cee AR pr 
Gz 2692 


, 2/7 
~—ahaek (12) 
x; xy 


Now x’ and x,’ denote the ionization of the Call present at time ¢ and at 
time 7,, respectively, and the ratio x’: x,’ is a measure of the importance of the 
recombination of CalIII to Call. In clouds for which k<11, x:x, remains 
approximately unity (see Table II), hence also x’ : x,’~1, and there will be little 
change in the ionization of calcium. There is sensible recombination in denser 
clouds, and some typical values of x’ are given in Table IV. Here we know 
that x,'~1, by the results in Table I. 

The abundance of CaII ions in the cloud may now be found. As in Paper I, 
Section 4.4, two different compositions are assumed for the stream. 

(c) Observation of the anomalous absorption.—There is effectively no 
recombination when k < 11, and such thin streams will remain invisible in calcium 
light both during and after the flare. It may be possible to observe denser streams 
of corpuscles and, using the same working as in Paper I, Section 4.4, we find 
the following values for the probable depth of the absorption band at two 
representative times. 
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TABLE IV 
Ionization of Ca II at two typical times after the end of the flare 
Time nig lonization of Ca II 
ionization of from Sun Pa hewe 
hydrogen (sec.) (cm.) ” 3 
2X 108 3°2 x 10! 0°98 0°84 
0°80 


10° 1°5 X 1018 0°97 0°80 


3°2 X 10! 
1°5 X 101% 
3°2 X 10! 


1°5 X 1018 


TABLE V 
The expected abundance of Ca II ions at two typical times after the end of the flare for clouds 
of two different densities 
Initial Time Distance Density of Ca II ions in the cloud 
ionization from Sun Composition A Composition B 
of hydrogen (ese.) (cm.) k=12 k=13 k=12 k=13 


2X 108 3°2x10" 46 3°5 x 107 1'2 xX 108 16 1*3 X 107 
0°80 


msxit0% 18x107-* 1°2 19X10 1°3x 10% 
3°aXx10"" 3°gx< 10° 1°3 x 10° 1°6 
rsxX10% 18x107% 1°53 1°99 X 1075 

2x 10° 3°2X 10% §«2°3 x 108 1°4 x 108 24 


108 asxz0"% 6 a%4xX10°° 1% 2°5 x 107% 
Aon Unséld’s composition, . eee 
B on Strémgren’s soma a i quoted in I, Section 4.1, Table I. 
, 


TasBLe VI 


Percentage depth of absorption near the H and K lines of calcium after a flare, 
for dense clouds 


Initial ; Distance Percentage depth of absorption 
ionization from Sun Composition A Composition B 
of hydrogen ‘ (cm.) k=12 k=13 k=12 k=13 


3°2 X 10 89 99+ — 20 
0°80 
1°5 X 1078 I 
3°2 x 104 99+ 
1°5 X 103% I 
2 X 108 3°2 x 10! 96 99+ 


10° 1°5 X 101% — o 


Values of less than 1 per cent are denoted by a dash. 
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The new values for the depth of absorption are rather lower than those found 
in the earlier paper (Section 4.4, Table VI). Even dense clouds for which k = 12 
can remain permanently invisible and so the chance of observing the absorption 
becomes very small. 
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AN ATTEMPT TO EXPLAIN THE POLARIZATION IN Ha AND D, 
FOR PROMINENCES 


H. Zanstra 
(Received 1950 July 29) 


Summary 


Using Heisenberg’s theory as developed by Ohman for problems of this 
kind, it is shown that, for the line Ha, the hydrogen atom acts as a space 
oscillator for 35 per cent of the incident energy, and as an isotropic scatterer 
for the remaining 65 per cent. For the D, line of neutral helium the figures 
are 36 and 64 per cent. 

For the photospheric radiation an approximate value for the coefficient 
of darkening §, 0-67 is derived for the centre of Ha from Abbot’s observations 
and Evans’ residual intensities. For the D, line a value 6,=f=1-36 is 
adopted. The theory, given for a prominence of small total optical depth, 
predicts for the emission a degree of polarization of 1°6 per cent at the 
photosphere, and 7-5 per cent at an altitude of 50,000 km., as compared with 
Lyot’s observed values of 0°8 to 1°2 per cent in general. For the D,; 
emission the theory gives 2°6 per cent at the photosphere and 8-4 per cent 
at an altitude of 50,000 km., as compared with the observed 0-9 to 1°5 per 
cent. The excess of the theoretical above the observational effect may in part 
be attributed to other processes of excitation like fluorescence and ionization 
and re-combination. The fact that the plane of vibration is not parallel to the 
solar limb, but shows deviations lying between o° and 30° for Ha and o° and 
40° for Ds, is not explained by the simple theory. 





1. Introduction.—The following is an attempt to explain the polarization of 
emission lines of Hx and D, in prominences, observed by Lyot.* The necessary 
theory was largely developed by Ohman+ making use of Heisenberg’s theory of 
polarization} for resonance or fluorescence under incident light. He also 
expected polarization in prominence emission lines, but only for the part of the 
prominence which is well above the photosphere, since he did not think of the 
effect of darkening towards the limb on the incident photospheric light. Actually 
this limb darkening results in a preponderance of light leaving the photosphere 
normally as compared with that for larger angles @ with the normal and this may 
result in a polarization of the light scattered by the lowest part of the prominence 
as well. So there are two effects intermixed: that of the limb darkening and that 
due to the elevation of the prominence above the photospheric surface. Effectively 
they both increase the intensity of that part of the radiation leaving the photosphere 
normally. If now for instance an atom of the prominence acts like a classical 
oscillator, it is easily seen that the scattered light would be partly polarized, there 
being a preponderance of electrical vibrations parallel to the limb. As Ohman 
pointed out and illustrated, Heisenberg’s theory also enables us to work out the 
polarization for other types of lines, as in our case the hydrogen line Ha and the 
line D,, A 5875-61, 28P—3°D of neutral helium. In Section 2 Heisenberg’s 

* B. Lyot, Comptes Rendus, 198, 250, 1934; Bull. Soc. Astr. de France, 51, 203, 1937. 


+ Y. Ohman, M.N., 89, 479, 1924. 
t W. Heisenberg, Zeit. f. Physik, 31, 617, 1925. 
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theory is applied to these two lines. The atom can partly be considered as an 
oscillator vibrating along the electric light vector and thus re-emitting polarized 
radiation, and partly as an isotropic scatterer re-emitting unpolarized radiation 
uniformly in all directions, and the fraction of oscillator scattering and isotropic 
scattering is computed. 

In Section 3 the polarization for prominence material close to the solar surface, 
which is due to limb darkening only, is considered, assuming a small total optical 
depth of the prominence in the spectral line. In Section 4 the general case of 
polarization due to the combined effect of limb darkening and elevation above the 
photosphere is treated, likewise for small optical depth, using Minnaert’s formulae, 
originally derived for scattering of photospheric kght by free electrons, but 
applying equally well to oscillator scattering of a spectral line. The result of the 
theory is contained in Table IV and this theoretical degree of polarization turns out 
to be stronger than the values observed by Lyot. Section 5 contains a remark 
regarding the derivation of Section 2. 

2. Representation of an atom by a classical oscillator plus an isotropic scatterer, 
worked out for Hx and D,.—Heisenberg considers a beam of plane polarized light 
incident on an atom or ion of any type in any state and capable of exciting this to a 
higher state by absorption of acertain line. In order to investigate the polarization, 
he introduces a magnetic field in the direction of the electric vector of the incident 
beam. ‘This magnetic field is not actually there, but merely serves as an artifice by 
means of which the various directions of vibration of the scattered light can be 
predicted. He then shows, from considerations akin to the correspondence 
principle, that this auxiliary magnetic field does not alter the vibrations of the atom, 
apart from the slight shifts in frequency due to the Zeeman effect, so that the 
degree of polarization remains the same. The excitation to the higher state, viz. 
to the Zeeman components of which it consists, takes place exclusively by the 
vibrations parallel to the electric vector, that is the 7 transitions for which AM =o, 
M being the magnetic quantum number. When however the atom falls back to 
its original or to any other state, he assumes, following Breit, or Gaviola and 
Pringsheim, that the emission is entirely spontaneous, and so contains the various 
7 and o components with the relative intensities predicted for the Zeeman effect. 
If now one lets the magnetic field approach to zero, the various 7 and o components 
tend to the same frequency, but their intensities and directions of vibration are 
not altered, so that one finally knows the degree of polarization if no magnetic field 
is present. 

As a preparation to Table I, Heisenberg’s procedure is illustrated in Fig. 1a 
for the line *S,—*P,, of a doublet. Upon application of the auxiliary magnetic 
field along the light vector E (Fig. 1B), the lower state J’ =} is split up into the 
components with magnetic quantum numbers M’ = —} and +4, and the upper 
state J” =1} into components with M” = —1}3, —4, +4 and +14. Only the z 
components AM=o are absorbed, so that only M”=—}4 and +4 are excited. 
Upon return to the lower state, the two components of total intensity 7, =2+2=4 
and the right-hand and left-hand o components, each of intensity 7,=1, are re- 
emitted, the latter being observed along the field so that they have their maximum 
strength.* 

* In general, the total intensity of all 7 components viewed perpendicular to the field will be 
termed 7, and of all right-hand circular components or all left-hand circular components viewed 
along the field 7,. Some authors define the intensities of the « components as those observed in 


the direction perpendicular to the field, that is the y or x direction in Fig. 1B. If this is done, 
the intensities in Fig. 1A should be halved. 
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For the purpose of linear polarization, one may replace the two right-hand and 
left-hand circularly polarized vibrations o, and og, each of strength i,=1, by two 
linearly polarized vibrations 7,, perpendicular to each other and in the plane of the 
formerly circular motion, as is shown in Fig. 1B (see Section 5). Consider first 
Heisenberg’s case of a plane polarized wave incident along the x-axis and vibrating 
in the z direction. When observed along the y-axis, the degree of polarization 
will be (i,—7,)/(¢,+7,)=%=60 percent. A space oscillator on the other hand 
would vibrate entirely along E and thus only have the 7 component, and a polariza- 
tion of 100 per cent seen along the y-axis. 


M" =1 Vs + / 





M' -/, +/2 


Fic. 1a.—Zeeman pattern for the transition J'=%4 to J°=14. M’ and M” are the magnetic 
quantum numbers of lower and higher states. The intensities of the o components (viewed along 
the field) and of the » components (in brackets) are indicated. 


X; 
! 


G 


4 
4 
7 


y* 


Fic. 18.—Intensities of the three vibrations. E is the vibrating electric vector. 





The total intensity of the three vibrations is 7,+27,, so that the fractional 
intensity of the 7 and o component is 


ty 


1, + 2tg and in t+ 2g (1) 
which yields 4/(4+2)=§ and 1/(4+2)=4 for the line under consideration.* 


* In the above example, the intensities of all 7 components are the same, viz. 2. When these 
7 components differ in intensity, as is the case for *P,,—*D,, and *P,,—*Dg, of Table I and all 
lines but the first of Table II, each absorbed 7 intensity is re-distributed amongst the 7 and o 
emissions from the state excited by it. The intensities of the Zeeman pattern then serve for the 
absorbed 7 intensities as well as for the ratios of the re-emitted 7 and oa intensities. 
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This is shown in the sixth row, columns 3 and 4, of Table I for the line 2*S, — 3?P,, 
of Ha. Exactly the same applies to the line of the next row 2?P, — 3*D,, which, 
though not belonging to the same doublet as 27S, —3?P,, coincides with it, as 
is indicated by the bracket. In the fine structure of Ha this composite line has the 
theoretical intensity 7-08 of the second column, so that finally z, = 4-72 andi, =1-18, 
given in the last two columns, in the intensity unit of column 2, which is common to 
all lines of the whole fine structure. 


TABLE I 
Polarization of Hx, 46562°79 


(The values of the intensities 7, and i, (Fig. 1B) are first obtained for each line of the 
structure and then added to get their value for the whole line Ha.) 


Lines of Int. f. str. Fractional int. Int. Int. 
fine str. . ig + 2t6 7 o js te 
2*Py 3°S, o'2 1/3 1/3 0°07 0°07 
2°P,,—3°D,, are) 41/75 17/75 0°55 0°23 
~ 3*D,, 13/25 6/25 4°6, 2'I¢ 
—3*P, aa 1/3 0°38 
3°S, I 3 
—3*P,, = 1/6 4°72 I'Ig 
¢ 3°D,, / 1/6 
Line Ha 18°42 10°49 4°09 


The other lines of Table I need but little elucidation, since the same procedure was 
followed. ‘The lines are all arranged in column I according to increasing frequency, 
the separation between the two strongest lines being 0:3285 cm.~! or 0-1415 A., 
the same as would be produced by a Doppler shift of 6-45 km./sec. The intensities 
of column 2 are the theoretical intensities of the fine structure. They were taken 
from a paper by Williams* who gives various literature references, and are 
naturally based on the occupation of the upper levels according to statistical 
weights, which should approximately be applicable to our case of excitation by 
photospheric light. For the Zeeman patterns and intensities needed, the reader 
may be referred to the current formulae and tables given in Candlert, Vol. II, 
pp. 108 to 120. 

All intensities of the second and the last two columns are now added and yield 
the desired i, andi, for the line Ha. The degree of polarization for Hx in Heisen- 
berg’s case of observations perpendicular to an incident plane polarized beam 
becomes (10°49 — 4°0,)/(10-49 + 4°02) =0-4,. We shall use the result in a different 
form. According to Fig. 1B, we may replace the intensities by three intensities 7, 
along the three axes and one intensity 7, —i, along the z-axis, viz. the direction of E. 
Thus the atom may be replaced by a space oscillator and an isotropic scatterer, the 
energies scattered being the fractions of the absorbed energy 


1, —1g Ste 
ta + 2tg ly +2tg’ 


(2) 





*Cf. R. C. Williams, Phys. Rev., 54, 558, 1928. For the most recent literature the reader 
may be referred to H. Kuhn and G. W. Series, Proc. Roy. Soc., 202, 127, 1950. 

+ Cf. A. C. Candler, Atomic Spectra, Cambridge, 1937. A diagram for the fine structure of 
Hz is found in Vol. I, p. 42, Fig. 3°17. 





No. 5, 1950 to explain the polarization in Hx and D, for prominences 495 


in general, or numerically for Ha : 35 per cent scattered by space oscillators and 
65 per cent scattered isotropically. The ratio of the absorption coefficient o' of 
oscillator scattering to o for isotropic scattering becomes, in general, 


o 6-1, 
or - 


or, in the case of Ha, _ = 0°53. (3a) 
Go 


’ 


For the polarization of D, of neutral helium, we have Table II which has been 
obtained in the same way as Table I for hydrogen. The case is somewhat simpler 
than for hydrogen lines, since the line represents only one multiplet, so that the 
intensities of each component line contained in the second row can be obtained 
from formulae or tables as given in Candler *, Vol, II, pp. 94 to ror. 


TABLE II 
A5875-618 : 2°P,, .—3°Das, 2,1 
A5875'960 : 2°P,— 3°D, 
Component Int. Fractional int. Int. Int. 
lines tn + 28 7 o in to 
2°P,— 3°D, 23°8 I ° 23°8 0'0 
2°P,—3°D, 17°9 1/2 1/4 8°95 4°47 
2°P,—3°D, 12 17/50 33/100 0°41 O*4o 
2°7P,—3°Dz 53°6 17/30 13/60 30°36 116, 
2°P,—3°D, 17°9 17/30 13/60 10°l, 3°85 
2°P,—3°D,; —100°0 259/525 133/525 49°4 25°3 
D, line 214°4 123°! 45°7 


Polarization of the Ds line of Hel 


The total intensities 7, and 7, for the whole D, line are again given at the bottom of 
the table. The degree of polarization for an incident plane polarized beam 
observed perpendicularly would become (123-1 — 45°7)/(123°1 + 45°7) =0°46. If 
instead of this we replace the atom by a space oscillator and an isotropic scatterer, 
formula (2) leads to 36 per cent and 64 per cent of the incident energy respectively 
scattered in these two ways. The ratio of the absorption coefficients o’ for 
oscillator scattering and o for isotropic scattering (3) becomes 


, 


~ =0:56 (36) 
oC 


for the D, line. 

Strictly speaking, one should take the Zeeman components of the hyperfine 
structure, and this results in a reduction of the theoretical polarization for Hz.+ 
However, for helium the nuclear spin is zero and consequently the D, line has no 
hyperfine structure, so that the computation for this line is quite correct. 

3. The effect of limb darkening: theoretical polarization of Ha and D, for low 
prominence material.—We now proceed to the theoretical excitation of a prominence 
by photospheric radiation, but assume the prominence material to be close to the 
solar surface, so that only the effect of limb darkening but not the effect of elevation 
comes in. 

For the continuous background of the photospheric spectrum at Ha A6562-82, 
the coefficient of darkening 8 = 1-20 follows from Abbot’s observations of intensity 

* Loc. cit. 
t A. G. C. Mitchell and M. W. Zemansky, Resonance Radiation and Excited Atoms, Cambridge, 


1934, p- 283. 
35 
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1(@) at cos @=1 and 0-312, or 8=0-g1 from the observations at cos 6 =1 and 0-661 
(cf. M.N., 101, 250, Fig. 2, 1941). For the excitation of Hz in prominences one 
should take the photospheric intensity /,(@) at the centre of Ha. There might bea 
displacement due to Doppler effect but, for instance, for 10 km./sec. this would be 
only 0-2 a., and presumably not carry one far away from the core. If r be the 
residual intensity at the centre of Ha, expressed in units of the continuous 
background, we have for /, and the corresponding coefficient of darkening £,, 


1,(8) = 1(8)r(8), (4) 
1,9) _1+B,cos@_ B, (1 —cos 6). (5) 
1,(0) 1+, 1+, 
For r(@) we use the values determined by Evans* in the second row of Table ITI 
for the cos @ of the first row. These are reduced to r(@)/r(0) in the third row. 
TABLE III 
Centre of Hx: Values of 1,(@)/I,(o) 
cos@=1:000 0°741 ~=—-0"364 
Evans r(0)=0'179 0°230 = 0224 
r(9)/r(o) = 1-00 1°29 1°25 
I(9)/I(o)=1-000 = 0875 0658 
Abbot 1@ilfcjm100 80:13 O82 


0-9 0-8 0-7 0-6 0-5 


)1(9)/I(o), x 1,(8)/T,(©) 


—> 


0-1 0-2 0-3 0-4 0-5 0-6 
> (l—cos @) 
Fic. 2.—Coefficient of darkening B, at the centre of Ha. 


Circles: Abbot’s observations, reduced to the crosses by means of Evans’ residual intensities. 
Straight line: Least squares solution for determination of B,. 


Erratum ; Ai further cross should have been plotted at cos 0= 1000 and 1,(@)/1,(0)= 1-00 


*D. S, Evans, M.N., 100, 156, 1940, Table XIII, 
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Abbot’s observations of 1(8)/J(0)* are plotted in Fig. 2 as circles, from which plot 
the values of the fourth row are obtained. By multiplying them by the previous row 
the desired values of /,(6)/J,(0) are obtained, given in the last row of Table III 
and represented by crosses in Fig. 2. Evans remarks that the experimental error 
in r is of the order of 3 percent. The scatter of the three crosses on the right is 
therefore largely experimental, while the maximum indicated by the second cross 
is probably real. However, we wish to use the linear relation (5), which can 
formally be done by introducing, for 6 =o0, the value 7(0) instead of the measured 
7(0), so that equation (5) assumes the form 


1,(8) _ r(o) B, 
7,(0) aL aa |, ) 


realizing that the difference r(0) — r(0) may only partly be an observational error. 
A least squares solution of (6) with the variables /,(@)/I,(0) and 1 — cos @ yields 
the straight line of Fig. 2 with 


and B, =0:67, (6a) 


for Ha. 

Thus £, for the line centre is smaller than the corresponding f of the adjoining 
continuous background. 

Now consider the theoretical effect to be expected for a prominence near the 
photosphere having a small optical depth in Ha. The radiation observed emerges 
at an angle 6=9go° with the normal to the photosphere. ‘This case of scattering 
material close to the photospheric surface was worked out by the writer.t We 
introduce o’, the coefficient of scattering per cm. by space oscillators, and a, the 
same quantity for the radiation scattered isotropically. 

Then the intensity of the oscillator radiation which is emitted by 1 cm.° of the 
prominence has the components ¢;, with the electric vector parallel to the solar 
surface and e, along the radius, given by 


ee ee _1,(0)_ 
ays (7+B,)” o'(4+3 35 P31), (7) 
I,(0) 


+B)” 


: 1,(0) 
=2 HO oh + Sh). (8) 
These are the formulae (4) and (5) of our former paper expressed in our present 
notation. The suffixes r and / mean that the plane of polarization is radial or parallel 
to the limb, but we now generally prefer to use the planes of vibration indicated by 
the capitals L and R. It should be remembered that a vibration parallel to the 
limb can also be indicated by its plane of polarization perpendicular to the limb. 
The intensity of isotropic radiation has the two equal components ¢€, and ex 
whose intensity is given by 
cr=€n=$ | h ~ - b. ws o(3 + 48;), (9) 
as follows from (5), dQ being the element of the solid angle. 


Epes Ka 2 


33 Pa — 338; sin? @) 


* A. Unséld, Physik der Sternatmosphdren, Berlin, 1939, p. 34, Table 3. 
1 H. Zanstra, M.N., 101, 250, 1941. 
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From (7), (8) and (9) follows the theoretical degree of polarization of radiation 
emitted towards the observer: 


, 


e 38 
_ (ex t+er)—(ert+€r) _ = 


p it 4a ¥4 * Topas ‘ 

(e, +€,) + (ent €p) = (14 “5 p,) +(1+48,) 

For Ha we had o’/o =0-53 (3a) and (6a) 8, =0-67, then (10) yields p=0-016 or 

1-6 percent. This result is entered in the first row of Table IV of the next section. 

Now consider the line D, of neutral helium. For A5876, Abbot’s observations 

at cos @=1 and 0-312 lead to a coefficient of darkening for the continuous photo- 

spheric background 8 = 1-51 and at cos @=1 and 0-661 to 8B =1-22, average 1-36. 

The Utrecht Photometric Atlas* shows at A 5875-62 only a faint line of depth 7 per 
cent, and at A 5875-96 no lines at all, and so we take 


B, =B =1-36. (66) 


Substituting the values (3 5) and (6 6) into (10), we obtain p = 0-026 or a degree 
of polarization of 2-6 per cent for the D, line of neutral helium. ‘The result is 
entered in the first row of Table IV of the next section. 

4. The combined effect of limb darkening and elevation. Comparison with 
observation.—The foregoing section dealt with prominence material just above the 
photosphere. The more general case of material at a point P of a prominence at an 
altitude h above the photospheric surface of radius Reo will now be considered. 
This can be done by replacing the formulae (7) and (8) for oscillator scattering by 
somewhat more complicated expressions originally derived by Minnaertt for 
scattering of photospheric light by free electrons, but which apply equally well to 
oscillator scattering of an emission line. If @ be the angle with the normal of the 
emerging radiation which is observed and ¢ the angle with the normal for the 
tangent from P to the solar surface, so that sind =1/(1+A4/Ro), the equations are 





(10) 


{ 


€,—€n = Ao’ sin? 6 cos ¢ sin? ¢ 


cos? 


+ $B, [3 sin? d—1I+ iad (1 +3 sin* ¢) log Ih (11) 


at 
c= Ao! {(3 —cos d— — *) 


2 
J 


cos? 

sind 
' a I 

with the abbreviation A= 0) (3. : 

1+f,\16 

These are Minnaert’s formulae (6) and (7) in our notation, except that the 


scattering coefficient per electron he used has been replaced by the scattering 


coefficient o’ per cm. for oscillator scattering in the line, so as to make it applicable 
to our Case. 


+48, [ 5 +sint 6 (5—sin? 4) log ==}. (12) 


cos 


(13) 


*M. Minnaert, G. F. W. Mulders and J. Houtgast, Photometric Atlas of the Solar Spectrum, 
Schnabel, Amsterdam, 1940. 


+ M. Minnaert, Zeit. f. Astrophys., 1, 209, 1930. 
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Further, for the isotropically scattered radiation the integration over dQ takes 
place for the solid angle within the tangent cone subtended from P by the solar 
surface, so that (9) for isotropic scattering is to be replaced by 

€, =€x= Ao{3(1 — cos $) + §B,(1 — cos* 4); (14) 
with A given by (13). 

For h=0, or 0=¢=7/2, these equations reduce to those of Section 3. 

Table IV contains the values of the degree of polarization from (10) com- 
puted from these general formulae (11) to (14) for the elevations 0, 25,000 km. and 
50,000 km., the latter being approximately the height of the average prominence. 
The quantity p/100 representing the percentage of polarization is given for the 
case 6=90", where point P is seen as far as possible from the limb and the case 
6 =¢, where P is seen at the limb. As might be expected, the polarization increases 
with increasing height: for H« it rises from 1-6 to 7-5 per cent and for the helium 
line from 2-6 to 8-4 per cent as the altitude becomes 50,000 km. If the prominence 
is seen at the limb, in front of or behind the solar disk, the polarization becomes 
only slightly less. 


‘TABLE IV 


Theoretical degree of polarization (per cent) for a prominence of small total optical 
depth at various heighis h: lines Hx and Dy, 





Ha: p/100 at Hel(D3;): p/100 at 
max. dist. limb max. dist. limb 





© 00000 ‘6 6 2°6 2°6 
25,000 0°0359 5°6 


50,000 0:'0718 


I 
5'2 6°5 6-0 
7°0 8-4 7:2 
t) Obs. (Lyot) 


t. 
bs. (Lyo 
o'8-1°2 o'9g-1'°5 

















In his 1934 publication Lyot’s conclusions from forty prominences are that 
prominences near the solar equator show a polarization in Ha of about 0-3 per cent 
with vibration approximately parallel to the Sun’s limb, or plane of polarization 
perpendicular to the limb. For latitudes between 36° and 44° the polarization 
becomes larger, between 0-5 and 1-2 per cent, but the plane of polarization is no 
longer perpendicular to the limb, but may deviate towards the north for prominences 
on the northern hemisphere, the deviations being between 0° and 22°. In three 
cases, contiguous prominences at high latitudes showed very different planes of 
polarization. Dense prominences generally showed uniform polarization, but 
in one case of a diffuse cloud the values varied much from point to point. The 
observations were made during a period of feeble solar activity and he judges the 
material to be too limited to allow general conclusions. 

In his 1937 note discussing observations of the year 1935 he states that polariza- 
tion of Ha generally is between 0-8 and 1-2 per cent and deviations of the plane of 
polarization occur up to 30°, in the same sense as above. ‘These values of p/100, 
which occur most frequently, are entered at the bottom of Table IV as “ observed”. 

For the D, line of HeJ, Lyot’s 1935 observations give polarizations between 
o-g and 1-5 per cent, but deviations of the plane of polarization up to 40°. 
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Prominences generally are observed at an altitude between 25,000 and 
50,000 km. The theoretical degree of polarization is therefore five to six times the 
maximum of the observed value. The discrepancy might partly be attributed to 
the fact that the theory has been worked out for the mechanism of resonance, 
whereas, in addition to this, the lines are excited appreciably by fluorescence and 
ionization and re-combination, while also electron impact might perhaps play a 
part. One might also think of other depolarizing agencies like phase-disturbing 
collisions *, while in the case of Ha the assumption of low total optical depths 
probably is not justified, and the hyperfine structure was not taken into account 
thus far. However this may be, it is gratifying that the theory is more than 
sufficient to explain the polarization observed by Lyot, although, in the simple 
form it is given in the present paper, the increase of polarization towards higher 
latitudes and the deviation of the plane of polarization that he observes are as yet 
unexplained. 

5. Remarks.—Considering the passage of circular into linear polarization 
represented in Fig. 1B, the following may be remarked. We take the average E* 
as a measure of intensity, E being the value of the vibrating electric field. Let 
E, =const. be the electric vector for each of the circular vibrations o, (clockwise) 
and og (counterclockwise), then the energy of each is E5=7,. When the auxiliary 
magnetic field is removed, the angular frequencies w of these two circular motions 
become the same, but their phase difference remains arbitrary. ‘Their resultant 
then is a linear vibration in the xy-plane of arbitrary direction, amplitude 2E, and 
energy }.4E;=2E;. Statistically the average energy of vibration along the 
y-axis is then equal to that along the x-axis, or Ej, viz. i,. The assumption made 
in the text is therefore justified. Also, the third component 7, has an arbitrary 
phase difference with the components 7,, so that statistically speaking the energies 
of the three components may be added to 2, + 27,, as was assumed in the text. 


‘The writer is indebted to Professor M. Minnaert of Utrecht and to Mr P. F. A. 
Klinkenberg of the Zeeman Laboratory in Amsterdam for valuable information. 


Astronomical Institute, 
University of Amsterdam : 
Fuly 1950. 


* Cf. H. Zanstra, M.N., 101, 266, 1941; 106, 225, 1946. 
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The following results are in continuation of those given in M.N., 109, 481, 
1949, and are derived from the measurement at Greenwich of photographs taken 
at the Royal Observatories of Greenwich and the Cape and at the Kodaikanal 
Observatory, India. ‘There were six gaps in the series. For these days, copies 
of original solar negatives were together supplied by the Mount Wilson 
Observatory, California, and the U.S. Naval Observatory, Washington, D.C. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the year 1945; the means for each 
year as a whole are included in Table II, which summarizes the yearly values 
from 1933 to 1945. 

TABLE I 


Mean Daily Areas 


Projected* Corrected for Foreshorteningt 
No. Rotation Days 
of Commenced Photo- 
Rotation Ried « 


1944-45 








~ 


Whole 


Umb Whole 
graphed ~™°T® — Spots 


Spots Faculae 


Faculae Umbrae 


1221 
1222 
1223 
1224 
1225 
1226 
1227 
1228 
1229 
1230 
1231 
1232 
1233 
1234 


Dec. 
Jan. 

Feb. 
Mar. 
Apr. 
May 
May 
June 
July 

Aug. 


Sept. 


Oct. 
Nov. 
Dec. 


18°74 
15°07 
11°42 
10°75 

7°95 

4°30 
31°52 
27°72 
24°92 
21°15 
17°40 
14°69 
10°99 

8-30 


66 
36 
62 
161 
123 
51 
103 
140 
go 
29 
237 
163 
99 
73 


358 
184 
323 
928 
669 
246 
541 
753 
480 
150 
1365 
935 
554 
374 


704 
486 
516 
631 
697 
729 
955 
866 
779 
556 
739 
IOI9g 
1139 
979 


56 309 
29 144 
47 241 
117 674 
93 499 
44 220 
82 423 
102 554 
66 349 
24 126 
179 1012 
135 779 
75 422 
56 294 


867 
592 
639 
744 
835 
895 
1135 


: 1088 


936 
690 
929 
1233 
1364 
1168 


* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 


Table III gives for each rotation in the year 1945 the mean daily area of the 
whole spots (corrected for foreshortening) and the mean heliographic latitude of 
the spotted areas for both northern and southern hemispheres. 

The mean heliographic latitude of the entire spotted area and the mean 


distance from the equator of all spots are also tabulated. 
the year 1945 are included in Table IV. 

Tables II and IV are in continuation of similar tables given in Monthly 
Notices; for the years 1874 to 1888—4g, 381-382; 1889 to 1902—63, 465-466; 


The mean values for 
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1901 to 1914—76, 402-403; 1913 to 1924—85, 1007-1008; 1923 to 1933—94, 
870-871, 

The rotations in Tables I and III are numbered in continuation of Carrington’s 
series. 

The chief features of the record are as follows :— 


(1) The steep rise in sunspot activity that began after the rather shallow and 
brief minimum at the epoch, 1944.2 (Tables I and II). The vigorous rise is 
particularly well shown by the mean areas of faculae for each solar rotation given 
in Table I. 

Taste II 
Mean Areas and Heliographic Latitudes 
Mean Daily Areas 


No. of Days Projected* Corrected for Foreshorteningt 


Photo- Without ,, _. Whole a . _ Whole ss 
graphed Spots Umbrae Spots Faculae | Umbrae Spots Faculae 


24 120 225 18 88 267 
27 142 297 22 354 
152 810 gIo 117 624 1100 
288 1529 2106 217 2545 
509 2795 2856 375 3505 
501 2706 2578 371 3205 
404 2163 1907 295 2349 
366 279 1434 1224 200 1522 
365 167 gI5 1041 120 1282 
365 16 112 591 666 79 809 
365 48 73 410 457 51 561 
366 157 30 160 281 23 338 
365 14 102 560 773 78 939 
* Expressed in millionths of the Sun’s disk. 
+ Expressed in millionths of the Sun’s hemisphere. 


365 3 
365 4 
365 2 
366 
365 
365 
365 


4 
5 
I 
° 
° 
° 
° 
° 
8 


(2) On only 14 days were no sunspots recorded, and even on these occasions 
faculae near the Sun’s limb were present. 

(3) From the onset of the new cycle up to the end of 1945, the southern 
hemisphere was more active in sunspots than the northern (Tables III and IV). 
In the latter phase of the old cycle, the northern hemisphere predominated in 
this respect (M.N., 109, 482, 1949). 

(4) The number and distribution of spot groups of duration (a) two days or 
longer: (5) one day, are as follows: 


(a) (d) 
North spots 56 13 
South spots 101 35 
‘Total 157 48 


The highest latitudes in which the longer-lived spots appeared were 35° north 
and 383° south. Corresponding latitudes for the 1-day spots were 29° north and 
43 south respectively, 
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TAaBLe III 


Spots North Spots South 
of the Equator of the Equator Mean 
: Mean _ Distance 
Mean a Mean Latitude from 
Helio- .1,,  Helio- of Entire Equator 
graphic graphic Spotted of all 
Latitude ‘ Latitude Area Spots 


No. Rotation Mean 
of Commenced Daily 
Rotation U.T. Area 


1944-45 ‘ ‘ 
Dec. 18°74 177 23°74 26°75 2°15 25°03 
Jan. 15°07 70 22°67 75 10°95 + 5°26 16°60 
Feb. 11°42 46 22°50 15°09 — 7°94 16°50 
Mar. 10°75 19 30°80 24°88 —23°28 25°05 
Apr. 7°05 4 27°35 12°02 —11°75 12°13 
May 4°30 20 28°14 19°90 -15°54 20°65 
May 31°52 23 28°59 24°21 —21°29 24°45 
June 27°72 148 22°29 20°27 — 8-87 20°81 
July , 24-92 332 23°32 20°2 + 21°27 23°18 
Aug. 21°15 27 25°24 17°15 7°96 18-90 
Sept. 17°40 388 8-48 20°86 — 962 16°12 
Oct. 14°69 242 24°37 23°19 ~ 8°43 23°56 
Nov. 10°99 146 24°69 77 19°55 - 4°30 21°33 
Dec. 8-30 125 23°34 I° 17°06 + 0:87 19°85 


TABLE [V 


Spots North Spots South 
of the Equator of the Equator Mean 
Mean Distance 
Mean Mean Latitude from 
Mean Helio- Mean Helio- of Entire Equator 
Days Daily graphic Daily graphic Spotted of all 
Photographed Area _ Latitude Area _ Latitude Area Spots 


365 86 
365 44 


3654 No Cycle 19 


_~ 


New Cycle 26 

365 205 
366 463 
365 1317 
365 899 
365 649 
366 499 
365 427 
365 252 
365 249 
6 t= Cycle 249 
3 New Cycle ° 
366 42 
Old Cycle 7 

3664 Xo Cycle 35 
365 121 


New 


wn 


yseeu est 


a 


~~ om 
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(5) The mean weighted latitude of all spots was 20°-22 (Table IV). The 
characteristic overlap in latitude of old and new cycle spots at the minimum epoch 
practically ceased in 1945. Only eight spots of the old cycle were indicated by 
the Mount Wilson observations of magnetic polarities. These few low latitude 
sunspots (mean latitude 7°) do not materially affect the mean latitude for the year. 

Other notes on the sunspots of 1945 will be found in M.N., 106, 76, 1946. 

An appended table gives the mean daily areas of whole spots and of faculae 
(corrected for foreshortening and expressed in millionths of the Sun’s hemisphere) 
for each calendar month of 1945. 


Monthly Mean Daily Areas of Sunspots and Faculae 


1945 
Month Spots Faculae Month Spots Faculae 





Jan. 182 2 July 06 1068 
Feb. 92 2 Aug. 816 
Mar. 515 2 Sept. 826 
Apr. 585 807 Oct. 1148 
May 283 963 Nov. 1365 
June 399 1083 Dec. 1195 





Royal Greenwich Observatory 
Herstmonceux Castle, Sussex: 
1950 June 9 
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